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HYDRAULIC DESIGN OF HOLLOW- JET VALVE STILLING BASINS» a J 

nd A. J. Peterka,2 F, ASCE 


fal By G. L. Beichley,! M. ASCE an 


SYNOPSIS 


Hydraulic model and prototype test s made to generalize and prove the hy- © 
draulic design of a new type of stilling basin which utilizes the hollow- jet 
valve for discharge control are described, Dimensionless curves are derived — 
7 from model data and are used to define the important dimensions of the beste 
4 for the usual combinations of valve size, operating head, and discharge. Sam- 
_ ple problems are presented to illustrate the use of the design curves and the 
_ general hydraulic design procedures. . Prototype tests on the Boysen and Fal- 

_ con Dam stilling basins are described and analyzed to help establish the reli- -_ 
_ ability of the recommended basins. Basin dimensions obtained from individual — a 
model tests on six stilling basins are shown to compare favorably with the _ 7 

; mensions obtained from the dimensionless curves and “methods given in this — 


and is used to ‘Glesigate hydraulic energy at the downstream end of an outlet 
works control structure. To reduce cost and save space, the stilling basin is 
usually constructed within or adjacent to the powerhouse ‘structure as shown — 


in Figs. 1 and 2. The hollow- jet valve, Fig. 3 3, controls and regulates the ‘flow. 


- Note.—Discussion open until February 1, 1962. To extend the closingdate one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part _ 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- So- 
ciety of Civil Engineers, Vol. 87, No. HY 5, September, 1961. 
_ 1 Hydraulic Engineer, Bureau of Reclamation, Denver, Colorado. . 


2 Engineer, Bureau of Reclamation, Denver, Colorado. 
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HORIZONTAL SECTION THROUGH 
Reservoir normal ws. El. 4725.00 
Outlet Works discharge 1320 cfs 
valves 100 percent open. 
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POWERHOUSE SECTION-THROUGH OUTLET STILLING BASIN 
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q FIG. 1. —BOYSEN DAM OUTLET WORKS STILLING BASIN: AND POWE 
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FIG. 2- LLOWTAIL DAM PROPOSED OUTLET WORKS STILLING BASIN AND 
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STILLING BASINS 


Regardless of the valve opening or head, the outflow has the ss same pattern, an 
annular or hollow jet of water of practically uniform diameter throughout its — 
length, Fig. 4. The stilling basin is designed to take advantage of the hollow-. 

shape; solid jets cannot be used inthis basing 
_ The hollow-jet valve was developed by the Bureau of Reclamation in the 
early 1940’s to filla need for a dependable regulating valve. _ The design was 
accomplished with the aid of a complete 6-in.-diameter hydraulic model and a 
sectional 12-in.-diameter air model. . These models were tested in the foc f 
of Reclamation Hydraulic Laboratory. To evaluate the valve characteristics 
at greater than scale heads, a 24-in.-diameter valve was tested at Hoover Dam 
under heads rar ging from 197 ft to 349 ft. 


(a) V fully 


(b) Valve 50 percent open 


4. 4. —SIK- INCH HOLLOW-< JET ‘VALVE DISCHARGING 


Piezometer pressure measurements, thrust on the valve 
- needle, and rates of discharge were studied in both field and laboratory tests. 
It was found that the hydraulic characteristics of the larger valves could be > 


_ predicted from the performance of the smaller model valves. From these 
tests and investigations of prototype valves up to 96 in. in diameter, the valve 
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‘te Cavitation damage, found on a few of the ‘many ones valves in use, was. 
_ minor in nature and was caused by local irregularities in the body casting and 
by y misalinement of the valve with the pipe. These difficulties have been elim- 
_ inated by careful foundry and installation practices. On one installation, dam- 
age that occurred on the cast iron valve support vanes may have been caused 
abrasive sediment in the water. . The design itself is cavitation free. 
“ea _ Because a large valve operating at high heads can discharge flows having 
an energy content of up to 150,000 hp,a stilling basin is usually required down- — 
stream from the valve. In early designs, the valve was discharged horizon-— 
_ tally onto a trajectory curved f. floor which was sufficiently long to provide a uni- 
_ formly distributed jet entering the hydraulic jump stilling pool. This resulted | 
in an extremely long structure, twice or more the length of the basin recom- 
mended herein. When two valves were used side by side, a long, costly divid- 
ing wall was also required. Hydraulic 1 model tests showed that the basin | length» 
could be reduced more than 50% by turning the hollow- jet valves downward and 
using a different energy dissipating principle in the stilling basin. The first 
stilling basin of this type was developed for use at Boysen Dam, a relatively 
low-head structure . Basins for larger discharges and higher heads were later 
_ developed from individual hydraulic models of the outlet works at Falcon, Yel- > 
lowtail, Trinity, and Navajo Dams. It became ‘apparent at this time that gener a 
alized design curves could be determined to cover a wide range of operating 
heads and discharges. Therefore, a testing program was initiated to provide 
the necessary data. A brief description of the individual model tests made to — 
the basin 1 type is given in the following section. Table 1 givesa sum- 
mary of basin dimensions, valve sizes, test heads, and discharges for. these 


Dam.—In the Boysen Dam studies, ‘a series of basic tests” 
- was made to determine the optimum angle of entry of a hollow-jet into the tail 
water. For flat angles of entry, the jet did not penetrate the pool but skipped 7 
is along the tail water surface. For steep angles, the jet penetrated the pool but 
_ rose almost vertically to form an objectionable boil on the water surface. When 
the valves were depressed 24° from the horizontal, Fig. 1, anda 30° sloping 
floor was placed downstream from the valve to protect the underside of the jet 
from turbulent eddies, optimum performance resulted. i The ‘submerged path 
of the: valve jet was then sufficiently long that only a minitum boil rose to the 
surface. The size and intensity of the boil were further reduced when converg- 
ing walls were placed on the 30° sloping floor to protect the sides of the jet 
until it was s fully submerged. ‘The converging walls have another function , how- 
ever; they cc compress the hollow-jet between them to give the resulting thin jet 
greater ability to penetrate the tail water pool. Sudden expansion of the jet as 
it leaves the converging walls plus the creation of fine grain turbulence in the 
_ basin account for ‘most of the. energy losses. in the flow. Thorough breaking-u -up 
of the valve jet within the basin and good velocity distribution over the entire 
flow cross section account for the low velocities leaving the basin. © Fig. § 
shows the performance of a hollow-jet basin both with and without the converg- 
_ Pressures on the inside face and downstream end of the converging walls 
were measured to determine whether low prensures which might induce cavi- 
a were present. The lowest pressure, -measured on the end of end of the wall, I, 
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was” 3 ft of water 2 above atmospheric; therefore, cavitation should not occur. 


- Pressures measured. on the sloping floor, and under and near the impinging 
jet, were all above atmospheric. Maximum did not exceed one- 


fourth of the total at the valve. baie 
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in each box were developed from individual model studies, 
BD Valve tilt 24°; inclined floor 30° in all cases. Te; on 
See Figs, 1, 2, 6,7,8,9,and1l, 
d Special case, for structural reasons. 
7 € Special case, for diversion flow requirements (dentated sill used and basin size in- 


Scour. downstream from the sill was wave heights 


=, were only 0.5 ft in the river channel. A vertical traverse taken near the end 
showed velocities to be about 5 5 fps, ‘decreasing to 


; 7 2 Upper values in each box were calculated from Figs. ill through 1 15 ; 
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— Falcon Dam. =In the Falcon Dam tests, two separate basins were developed, — 
one for the United States outlet works and one for the Mexican outlet works, 
Figs. 6and 7. In these tests, the basic concepts of the Boysen design were © 
proves to ‘Satisfactory for greater In addition, it was confirmed 


on with short converging wells 


action with recommended convergi 
= 


HOLLOW-JET VALVE STILLING BASIN WITH AND WITHOUT 
‘ that dentils on the end s sill wer were not necessary and that the center dividing wall 
: _ need not extend the full length of the > basin. JA low 2:1 sloping end sill was suf- 
ficient to provide minimum scour and wave heights. Maximum on 
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STILLING BASINS 


the floor beneath the impinging jet ‘were found to ‘be about one- third of the he oat 
head at the valve, somewhat greater than found in the ‘Boysen tests, but still 
Yellowtail Dam.—In the Yellowtail Dam model studies, the head and 
_ charge were both considerably higher than in the Boysen and Falcon tests. Be- 


‘- 7 cause of the high velocity flow from the valves , it was found necessary to ex- 


a and to ‘reduce the wall gap | to about one- .-quarter of the basin width. "These re- 
 fnemente improved the stilling action within the basin, Fig. 5 (c), and made it 

_ possible to further reduce the basin n length. Scour was not excessive, and the - 
water surface in the downstream channel was relatively smooth. Pressures on 
the converging walls and other critical areas in the basin were found to be 

Trinity Dam. —The Trinity Dam outlet works utilized a head almost 4 times 
greater rong a discharge 5 times greater than at Boysen Dam, In the develop- 
ment tests, it was found that the performance of this type of basin would be 

satisfactory for extremely high heads and discharges. Although several a 

- ations in the basin arrangement | were investigated, no new features were in-— 

_corporated in the design. Fig. 8 shows the developed design. 

Navajo Dam.—The experimental work on the Navajo outlet works was com- 


ce by the fact that the hollow - jet valve basin, Fig. 9, had to first serve 


the pr oper appurtenances in the temporary basin ‘convert it toa per- 
_ manent outlet works basin - The development tests indicated that a larger than 
at necessary basin does_ not in itself guarantee satisfactory performance of the 
“hollow- -jet valve basin. fb - Best outlet works performance was obtained when the _ 
«EEE basin was reduced in size to conform to the optimum size oil was 


for the permanent structure. © ‘Since the Navajo Dam outlet works model was 
viilable both during and after the generalization tests, the model was used 


bot! toaid in obtaining the generalized data and to that the design: curves) 


‘Because ‘development basins had reached a point where 

the general arrangement of the basin features was consistent, and because the 

basin had been proved satisfactory for a wide range of operating conditions, a 
- testing program was inaugurated to provide data for use in generalizing the : 

basin design. The purpose of these t tests was to provide basin dimensions and 
hydraulic design procedures for any usual combinations of valve size, dis- 
charge, and Operating head. The main purpose of this paper is to describe 
_ these tests, explain | the dimensionless curves which are ‘derived from the 
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pr os -jet valve basins, that fit the dimensionless ye in the gen- _ 
eral study, will perform as well in the field as predicted from the model tests. _ 
Equipment.—The outlet works ‘stilling basin model shown in 


was used for the generalization tests . The glass-walled testing flume con- 
tained two stilling basins separated by. a dividing wall. is he right-hand basin 
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length, width, a and depth requiremetns; both basins were used to study | the per- 
formance with and without flow in an adjacent basin, 

_ The glass panel permitted observation of the stilling action and the flow _ 
currents within and downstream from the basin. The length, , width, and depth — 
of the basin were varied by inserting false walls or byr moving the basin with- — 
in the test box. The tail box contained an erodible sand bed to > represent the 


__ The test valves were exact models of a prototype valve in that the flow sur-_ 
faces were exactly reproduced, and could be opened and closed to any partial — 


opening. The models were 3-in. valves machined from bronze castings. ty 


FIG.  10.-HOLLOW- JET VALVE ‘STILLING BASIN MODE L USED 


_ GENERALIZATION TES TESTS 


Discharges were / measured using ‘calibrated venturi meters permanently in in- 


controlled with a hinged tailgate in the tail box. Tail water elevations were | 
determined visually from a staff gage on the tail box wall located approximate- _ : 

bd 62 valve diameters downstream from the valves. eh ain 
_ Preliminary Procedures. —The was begun by » the 
Falcon, Yellowtail, and ‘Trinity outlet 


‘fan on these dimensions, a model was constructed as shown in Fig. 11, us- 
ing the 3-in, valve dimension to establish the absolute model size. More weight a 


was given to the Yellowtail and Trinity basins because they | were developed for 
Pei heads and contained refinements in the converging wall design which 
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= provide pra practical ciliate Limits : for the tests, ate 3-in. model was as 
‘sumed to represent an 84-in. prototype valve | making the model scale 1: 28. = _ 
_ Discharges of 2,000 sec-ft to 4,000 sec-ft with one valve open 100% were con- | 4 
sidered to be the usual design discharges for a valve of this size. To produce © 
these discharges, heads of 100 ft to 345 ft of water at the valve would be re-— 


Initial tests were made with the stilling basin apron longer than necessary .? 


“és with no end sill in 1 place. For a given discharge, the ideal depth of tail wa- 

a ter was determined from visual inspection of the stilling action as it occurred B 

over a range of tail water elevations. For each ideal tail water determination, © ; 

a the ‘minimum length of concrete apron was estimated after an inspection ofthe | 

- flow currents in the model had indicated where an end sill should be placed in 7 
the prototype. Confirming tests were then conducted successively on a repre- 8 
sentative group of basins having the apron lengths previously determined and : 

. having an end sill at the end of the apron. . Adjustments were then made as ne- 
 cessary to the preliminary values to obtain final ideal tail water depths and © 
apron lengths. In the latter tests, the height of the valve above the maximum 
water elevation was adjusted to simulate a typical prototype installation. 
Similar tests were then made with the valve open 75% and 50%. Finally, a 
_ series of tests was made to determine the ideal width of stilling basin and the 
4 _ range of widths over which satisfactory performance could be expected. a ad 
be a _ Preliminary Tests.—In a typical test, the desired discharge was set by 
: means of the laboratory \ venturi meters and passed through the hollow- jet valve 
or valves opened 100%. The tail water elevation was adjusted to provide the 
best energy dissipating action in the basin. . The optimum value, tail water | 
depth D in Fig. 11, was judged by the appearance and quality of the sieeeees ac- 

tion in the basin and on the smoothness of the tail water surface. 
_ Fordischarges of 2,000 sec-ft to 4,000 sec-ft, it was found that ‘the tail wa- es, _ 
ter could be raised or lowered about 3 ft (0.1 ft in model) from the ideal tail — 
water elevation without adversely affecting the basin performance. — 


the tail water depth beyond this margin reduced the efficiency of the stilling 
action and allowed the jet to flow along the bottom of the basin for a greater 
distance before being dissipated. This aiso produced surges in the basin and 
7 _ increased the wave heights in the discharge channel. Decreasing the tail wa- 
: - ter depth below the 3-ft margin moved the stilling action downstream in the _ 
dal basin and uncovered the valve jets at the end of the converging walls. This in- 
, _ creased the flow velocity entering the discharge channel and increased the 
_ tendency to produce bed scour. Uncovering of the stilling action also produced 
“objectionable sp splashing a at the upstream end of the basin. If the tail water depth | - 
was decreased further, the flow swept through the basin with no stilling action © 
_ having oceurred. The I latter tail water depth was measured and recorded a 
: the sweep-out depth Dg. _ These tests were made with the dividing wall extend-. 
Ms ies to the end of the basin, since this provided the least factor of safety against 
sweep out. Witha shorter wall, out occurs at a tail wa- 
_ With the ideal tail water depth ost for a desired ~— the action in | the _ 
was examined to determine the ideal length, L, of the basin apron, Fig. 11. 
The apron length was taken to the point where the bottom flow currents began 
to rise from the basin floor of their own accord, without assistance from an 7 
end sill, Fig. 5 (c). The water surface directly — and downstream from | 
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‘this point was fairly smooth, indicating that the stint action had been com- 
pleted and that the paved apron and training walls need not extend farther. _ = 
the preceding individual model studies, it had been found that when the basin was 
appreciably longer than ideal, the ground roller at the end sill carried bed mate-— 
_ rial from the discharge channel over the end sill and into the basin. If this ac- 
_ tionoccurred ina prototype structure the deposited material would swirl around : 
in the ene: end of the basin and cause abrasive damage to the concrete 
and end sill. It had also been found that scour tendencies in the discharge 
> channel were materially increased if the basin was appreciably ‘shorter than 
ideal. Therefore, the point at which the currents turned upward from the — 
: plus the additional length required for: an end sill, was s determined to be the opti- ei 
" any scouring tendencies would be reduced by the sloping end sill to be added later. 4 
Practical difficulties were determining the exact length of 


a ' to s D in a period of 15 sec to 20 secin the model. An average apron stl 


_ was therefore selected in the preliminary tests. _ For this reason, too, the end 
sill would help to neutralize the scouring tendencies which increased as the 
_ The depth D, sweep-out depth Dg, and length Lwere then determined for the - 
range of discharges possible with the hollow-jet valve open 75%, and finally 
a 50%, using the testing methods described in the preceding paragraphs. ~Par-— 
tial openings were investigated because the valve size is often determined for 
the minimum operating head and maximum design ‘discharge. the same 


R quantity - is discharged at higher heads, the valve opening must be reduced. Tt a 


may be necessary, therefore, to design the basin for maximum discharge with — | 
the valves opened less than 100%. When the relation between head and velocity 
in the valve is changed materially, the minimum required basin dimensions *] 


will be affected. _ The data for the partially opened valves are also useful in 


‘joan of the energy in 1 the flow at the valve. a energy may ybe represented 
by the total head, H, at the valve, Fig. 11. Therefore, to provide dimensionless 
_ data which may be used to design a basin for any size hollow-jet vaive, D, Ds, : 
and L values from the preliminary tests were divided bythe valve diameter d, 
7 te each variable was plotted against H/d. The resulting curves, similar to — 


— > 


those in Figs. 12, 13, and 14, were used tv obtain dimensions for a group of | 


model basins which were tested with the end sill at the end of the apron and -> 
with the valves placed to give the proper vertical distance between the valve a 
and the tail water. | For each model basin, a 3:1 upward sloping erodible bed, — q 
composed of fine sand, was installed cownstream from the end sill. The bed oo 
was kept sufficiently low that it did not interfere with tail water manipulation, — 
even when the tail water was lowered for the sweep- -out tests. Test procedure — 
_was essentially as described for the preliminary tests. 
 Ieste Depth and Length. —The preliminary depth curves for both ideal tail 
vel depth and sweep-out tail water depth needed but little adjustment. The 


preliminary basin lengths were found to be too — — es the 1e high he heads and) too 
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design flow conditions, 
__ Final Tests and Procedures .—The final tests were made to correct orveri- — . 
—— fy the dimensions obtained in the preliminary tests and to investigate the effect 
_of varying the basin width. Scour tendencies were also observed to help eval-___ 
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short for the lower heads, although both adjustments were relatively ‘minor. _ 
The adjusted and final curves are shown in Figs. 12, 13, and 14. 
_ It was observed that a longer apron than indicated by Fig. 14 was necessary 
when the tail water depth exceeded the tail water depth limit in Fig. 12. As the 
stilling action became drowned, the action in the basin changed from fine- 


grain turbulence to larger and slower m moving vertical eddies. The bottom flow 


the apron for a greater distance, thereby increasing the necessary length of 
basin. The action is similar to that observed in hydraulic jumps which are 
drowned by excessive tail water depths. _A moderate amount of drowning is 
tolerable, but it is important that the ideal tail water depth be maintained with- 
in stated limits if the best performance is desired. The tail water depth limits, | 
0.1 ft above and below the ideal depth, expressed in dimensionless form is 
0.4 d. this limit is exceeded, a model is 


greoter or less than 0. -_ 
H,ond d and ore defined in " 


Represent dota points shown in Figures IG and 


lad es Basin Width.. —To determine the effect of basin v width, tests on several ba- 
gins were made in which only the basin width was varied. It was found that the 7 
width could be increased to 3.0 times the valve diameter before the action be- 
came unstable, | The width could be decreased to 2.5 times the valve diameter 
before the stilling action extended beyond the ideal length of basin. However, 
the H/d ratio and the valve opening were found ‘to affect the required basin 
width as shown for 100%, 75%, and 50% valve openings 
4 Basin width is not a critical dimension but certain precautions should be ~ 
taken when selecting a minimum value. If the tail water is never to be lower — 
than ideal , as shown by the curves in Fig. 12, the basin width may be reduced 
3 2.5 d. If the tail water elevation is to be below ideal, however, the — 
_ values for width in Fig, 15 should be used. In other words, the lower limits © 
‘ for both tail water and basin width should not be used in the same structure. 7 
F The combined minimums tend to reduce the safety factor against jump sweep- | 
out and poor overall performance results. The basin width should hot be in- 


aetna above 3.0 dto ‘substitute for. some of the required length « or depth of 
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the basin, If unusual combinations of width, depth, , and length are needed to fit 
a particular space requirement, a model study is recommended. | 
: _ Basin Performance.—The six model basins shown operating in Figs. 16 and 
17 illustrate the performance to be expected from the recommended struc- - 


tures, The in 16 and 17 correspond to points shown 


‘the model appearance ‘the prototype basin and may help to pro 
_ Vide a visual appraisal of the prototype structure. Wave heights, boi! heights, — 
_ or other visible dimensions may be scaled from the photographs (using the 

scale shown in the photographs) and converted to prototype dimensions by 
"multiplying the scaled distances by the model scale. To determine the model 
_ scale, the prototype valve diameter in inches should be divided by 3 (the model 
; valve diameter). To determine which of the six photographs represents the > 
prototype in question, -H/dr atio should be. used to select the photograph — 

which most nearly represents the design problem. Itis permissible to inter - 

‘polate between photographs when 


Po. 

MOTE: Ds is the depth of toilwoter above the apron the 
_ flow from the valve first begins to sweep out of the basin 


0nd ore defined im Figure 


Center Dividing Walt —Prototype ‘stilling usually have two valves” 

placed a minimum distance apart, and alined to discharge parallel jets. It is” 
necessary, without exception, to provide dividing walls between the valves for 

P _ satisfactory hydraulic performance. When both valves are discharging with- 

out a dividing wall, the flowi in the double basin sways from side to side to pro-— 

duce longitudinal surges in the tail water pool. . This action occurs because 
the surging downstream from each valve does nct have a fixed period, and the 
resulting harmonic motion at times becomes intense. When only one valve is — 


_ discharging, conditions are worse, The depressed water surface downstream E x a 


from the operating valve induces flow from the higher water level on the non- ; 
operating side. Violent eddies carry bed material from the discharge channel — 
into the basin and swirl it around. This action in the prototype would damage > 7 
thet basin as well as the discharge channel. In addition, the stilling actionon 7m 
_ To provide acceptable operation with one valve operating, the dividing wall - 
should extend to three-fourths of the basin length or more. However, if the : ; 


two adjacent valves discharge equal quantities of flow | at all times, the length 
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W,H,and d are defined in Figure 


Represent dota points shown in 
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performance is for widths 


_ Good performance occurs over range o 


Best hydraulic 
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he of the center r dividing wall may be reduced to one- half of the basin length. The 
| margin against sweep out is increased, but the stability of the flow pattern is 
_ decreased as the dividing wall is shortened. In some installations, a full-— 
~ length wall may be desirable to help support the upper levels of a powerplant, 7 
“Fig. 1. If other arrangements of the center = are required a model study 


of the possibility of cavitation occurring within the valve. However, the valve 


Valve Placement.—A hollow - jet valve should not operate submerged because ause 


may be set with the valve top at maximum tail water elevation, and the valve 
not be ‘underwater at maximum discharge. The valve jet sweeps the tail 
_water away from the downstream face of the valve sufficiently to allow usual 
ventilation of the valve. | Z However, as a general rule, it is recommended that 
the valve be placed with its center (downstream end) no low lower than tail water — 


ie _Riprap Size.—A prototype basin is usually designed for maximum discharge, — 
= but will often be used for lesser flows at partial and full valve openings. For 


these lesser discharges, the basin will be larger than pers and in most 


‘design discharge, particularly those close to the design discharge, the ground 


a roller will tend to carry some bed material upstream and over the end sill in- a" 


to the basin. The intensity of this action is relatively mild over most of the 

i discharge range, and movement of material may be prevented by placing rip- 

rap downstream from the end sill. - Riprap, having 50% or more of the indivi- 
dual stones 24 in. to 30 in. or larger in diameter, should provide a stable chan- 
nel downstream from the end sill. The riprap should extend a distance D, or 
_ more, from the end sill. . If the channel is excavated and slopes upward to the 


a i natural river channel, the riprap should extend from the end sill to the top of - 

slope, or more. riprap should not be terminated on the slope. 

a The justification for choosing riprap as described is as follows: "Because 
of the fixed relationships between depth and width of basin , the crea veal 


, ‘eee velocities will therefore seldom exceed 7 fps to 8 fps and bottom velo- . 
cities 3 fps to 4 fps. To protect against these velocities, stones 10 in. to 12 

aa in, in diameter r would b be ample. However, the critical velocity for riprap sta- 

; - bility’ is the upstream velocity of the “ground roller which has a curved path 
and tends to lift the stones out of place, — Model tests showed that graded rip- 
_ ‘rap up t to 24 in, to 30 in. in diameter was sufficient | to provide bed stability. 

APPLICATION OF RESULTS 


a stilling basin for 1 hollow- jet valve ve discharging 


= 


1,300 cfs, ane (b) a double basin for 2 valves discharging 650 cfs each. In eal 


Stilling Basin Design.- 
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in is mie is Of Uischarpe, A IS the 
ae a inlet area to the valve, g is the acceleration of gravity, and H is the usable or : > 
_ total head at the valve with the valve center placed at maximum tail water el- 
oon s evation. In this example, the usable head at the valve is estimated to be80% of = = 


STILLING BASINS is 


From Fig. 3, for 100% valve opening: 


“ad 
= 25 sq ft 


in which d is the inlet diameter of the valve and also the nominal valve size, - 


Since nominal valve sizes are owe graduated in 6-in, innreanenmn, 


” Having determined the valve size and therefore the diameter of the supply | 
conduit, the probable head losses in the | system ‘from reservou to io valve may 
be computed. In this example, the computed losses are assumed to be 20 — 
which leaves 88 ft of head at the valve. Using Eq. 1,C is computed to be 0.61; 
_ from Fig. 3, , the valve o opening necessary to ugha the design discharge at the 
=e The basin depth, length, and | width — be determined from Figs. 12, 13, 14, 


14, 


how e depth rato 

therefore, th the apron. is 20. . below th the maximum tail water elevation, 


d 
atin 
length ofthe basinis 
For 89% valve Fig. 15 shows the width ratio 
width « of. the basin is 
y= 15 ft 
The of other components of the basin may be determined from 
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September, 
‘The tail water at which “the flow will s weep p from the basin may 


determined from Fig. 13. For 83% valve opening, the sweep-out ratio. 


be placed 2.4 ft higher to provide only 18 ft of depth — a ie 


a ‘If the tail water depth from Fig. 12 is adopted, the water surface aiiial will 
be similar to that shown in Fig. 16 (a), ‘since the H/d value of 16 in Fig. 16 ( (a) > 
is comparable to 14.67 in this example. If tail water depth 2 ft greater or less 
' _ than the ideal is adopted for the prototype, the water surface profile will be 
moved up or down accordingly. Water surfaces may be estimated by multiply- 
ing the variations shown in Fig. 16 (a) by the quotient obtained by dividing the 
- _ prototype valve diameter of 72 in. by the mdoel valve diameter of 3 in. Wave 
heights in the downstream channel will be considerably less as in 
_ _Two-valve Stilling Basin Design.—If two valves are to be used to discharge 
_ the design flow of 1,300 sec-ft,a double basin with a dividing wall is required. 
a _ The discharge per valve is 650 cfs, , and at 100% valve opening the vaive coef- 
ficient is 0.7, Fig. 3. _ The head on the valve is estimated to be 86 ft as in the 
- first emusic, fren Eq. 1, the inlet area of the valve is found to be 12.48 sq 
‘ft. 48-in. valve provides practically the exact area required. 
For this example, it is assumed the computations to determine 
losses have been made and that the estimated head of 86 ft at the valves is 
correct, Therefore, 100% valve opening will be to ‘pass the the design 


‘Using the methods { wei in detail in the first example: pee 


—_ 


4.06, from from Fig. 120 


| 


— 
— 
‘| i Since 20.4 ft of depth is provided, the basin has a safety factor against sweep- ’ 7 are 
ae €=—soutt of 4.2 ft of tail water depth. In most installations this is sufficient, but if _ 
a of desired, the apron elevation may be lowered 
ne 
al 
— 
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STILLING BASINS 

The tail water depth for sweep out is therefore 3.0 ft belowthe ideal tail water a) 
_ depth. _If more or less insurance against the possibility of sweep out is de- 
sired, the may be set lower aad amount 


| 


other discharges on the tail water range may need to be —m, - 


14.4, from Fig. 14 


the thickness of the center ‘aridine wall. The length of the center dividing wall 

_ should be three-fourths of the apron length or 43.5 ft long, Fig. 11. If it is 
certain that both valves will always discharge equally, the wall need be only 
one-half the apron length or 29 ft long. The hydraulic design « of the basin may 

ton _ If the tail water depth determined from Fig. 12 12 is adopted, the water sur-_ 
_ face profile for determining wall heights may be estimated by interpolating be- 

_ tween Fig. 16 (a) and (b). Water surface variations may be predicted by mul- - 
valees trom the by the ratio aay 3 


PERFORMANCE 
and Falcon Dam ‘outlet works stilling 1, 
and 7, fit the design curves derived from the generalized study quite well, and 
have been field tested and found a in an excellent manner. Table 1 


those obtained from the individual model tests. — - 

Boysen Dam.—The outlet works basin at Boysen is for 1 320 

“@  efs from two 48-in. hollow- jet valves 100% open at reservoir elevation 4725. 00. 

ce “ _ Design tail water elevation at the basin is 4616, 00. |The model performance of j 

this basin is’ shown in Figs. 

‘The prototype tests, Figs. 20, 21, and 22, were conducted with the reser- 

= at elevation 4723. 5 and with ‘the powerplant both operating and shut down. 
_ ‘spillway was not operating. The outlet works discharge was measured 

- at a temporary gaging station located about 1/2 mie downstream from the 

dam using a current meter to determine the discharge. Tail water elevations 
were read on the gage located in the powerhouse. 

_ _The prototype performed as well as predicted by the model - was con- 

: sidered satisfactory in all respects. However, the field structure entrained 
_ more air within the flow than did the model. * This caused the prototype flow to_ 
appear more bulky, cy, and “white water” extended farther into the “downstream — 
channel than was indicated in the model. _ A comparison of the model and pro- 

totype Figs. 19 19 and 22, illustrates this difference. _ Greater air 
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«FIG. (18. —BOYSEN DAM, LEFT VALVE OUTLET WORKS BASIN, DISCHARGING ING 
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Both valves fully open and both turbines oper 
elevation 4723.5. Tailwater elevation 4617. 
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4,570 cfs - res. elev. 300 approx. - T. W. elev. 181.2 
1:30 scale model +) 


_ FIG. 23.-MEXICAN OUTLET WORKS - FALCON DAM 
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comparisons, particularly when the ‘difference between model and | 


i: velocities is appreciable. In other respects, however, the prototype basin was _ | 
as good or better than predicted from the model tests. 
the initial prototype test, only the left outlet valve was 
powerhouse w was not operating. At the gaging station, the discharge was meas- 
ured to be 732 cfs after the tail water stabilized at elevation 4614. 5. (This is 
a greater discharge than can be accounted for by. calculations. It is presumed 
tha valve overtravel caused the valve opening to exceed 100% even though the — _ 
that val showed 100% open.) It was possible to descend the steel ladder, Fig. [i 
a to closely observe ai and photograph the flow in the stilling basin, Figs. 20 and 
21. The basin was remarkably free of surges and spray; the energy dissipating — 
action was excellent, There was no noticeable vibration at the valves or in the 
basin, The flow leaving the structure caused only slightly more disturbance in 
ts tailrace than the flow from the draft tubes when the turbines were operat- -y 


e _ Operation of the prototype provided an opportunity to check the air require- 
‘ments of the structure, which could not be done on the model. With the inspec- 
“tion cover removed, Fig. 1, the basin was open to the rooms above, / - move-_ 
_ ments through the inspection opening and in the powerplant structure were neg- 
ligible, which indicated that ample air could circulate from the partially Open = 
end of the stilling basin, Fig. 21.00 
= When both valves were discharging fully open, open, the tail water stabilized at. 
- elevation 4615. A discharge measurement at the gaging station disclosed that 
both val valves were discharging 1,344 cfs. * Since the left valve had been found to 
‘discharge 732 cfs, the right valve wa was discharging 612 cfs. 
_ The reason for the difference in discharge is that the 57- inch- inside- 
diameter outlet pipe to the left valve is short and is connected to the 15-foot- 
- diameter header which supplies water to the turbines, Fig. 1. ‘The right valve : 
_ is supplied by a separate 66-inch-diameter pipe extending to the reservoir. : 
a - Therefore, greater hydraulic head losses occur in the right valve supply line, 
_which accounts for the lesser discharge through the right valve, _ Although it 
was apparent by visual observation that the left valve was ‘discharging more 
than the right valve, Fig. 22, no adverse effect on the performance of the out- 
__ let works stilling basin or on flow conditions in the powerhouse t tailrace could _ 
The outlet works basin - performance w: was also observed with the turbines o_ 
ae operating and the tail water at about elevation 4617. No adverse effects of the a 
— outlet works discharge c on powerplant performance could be detected, Flow 
a conditions in the tailrace area were entirely satisfactory, Fig. 22. ‘Since the 
_ tests were made at normal reservoir level and maximum discharge, the still- _ 
q basin was subjected to a severe test. 
. i Falcon Dam. .—The outlet works basin inte Mexico side at Falcon Dam is 
. - designed to accommodate 4,570 cfs from two 90-in. valves or 2,400 cfs from | 
- one valve, with the valves 100% open and the reservoir at elevation 300. . The _ 
tail water elevation is 181.2 when the powerplant is discharging 5,400 cfs ) in 
- conjunction with both valves. The model performance of this basin is shown 
_ The outlet works basin on the ‘United | States side at Falcon Dam is designed 
to discharge 2,920 cfs from two 72-in. valves, or 1,600 cfs from one valve, 


with the valves 100% open and the reservoir at elevation : 310. Tail water isat 
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90-inch left ‘valve di ay 100 percent open 
2, 400 cfs - res. elev. 300 approx. - T. W. elev. 181.2 
1:30 scale model 


72-inch valves open 100 percent 
% 2920 cis - res. elev. 310 approx. - T. W. elev. 180. 8 a 


1:24 scale mo e model 
25.- _UNITED STATES OUTLET WORKS - FALCON DAM 
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72-inch right valve discharging 100 percent open 
1600 cfs - res. elev. 310 approx. - W. elev. 180. 
24 scale model 


72- “inch valves open 100 | percent 
2920 cfs - res. elev. 310 approx. - T. W. elev. 180.8 
scale model 
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inch outlet works open 100 percent discharging 
4,500 cfs approx. - T. W. elev. 183.6. 
gates 72 ‘percent open - load. 


72- inch outlet works valves open 100 percent discharging ee i 
3,000 cfs approx. - T. W. elev. 184.1. 

Turbine gates 72 percent open - 100 percent load. 


FIG. 30. —FALCON DAM MEXICAN & UNITED STATES POWER- 5 
PLANTS & OUTLET WORKS DISCHARGING AT RESER- 
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erating. ‘The model performance of this basin is s shown ; in n Figs. 25, 26, and 27, 
The prototype tests at Falcon, Figs. 28, 29, and 30, were conducted at near 
muons conditions ; the reservoir was at elevation 301.83, and the valves 


were open. 


The prototype valve were determined discharge curves based 
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31.—DEVELOPED BASIN 

Here, too, more white water was evident in the prototype than in the model. 

The greater ‘amount of air entrainment in the prototype, evident in the photo - 3 

graphs, caused bulking of the flow at the end of the stilling basin and a higher — 
water surface than was observed in the model, However , the prototype tail wa- 
ter is 3 ft to 4 ft higher than shown in the model photograph, and this probably ~ 


_ helps to produce a higher water eurtace boil at the end of the ba- 


September, 1961 HYS5 
th -va operation represents an emergency condition and , 
s the stilling basin tothe test, Figs. 28 and 29. All turbines 
owerplanics we) and 100% load during all tests. 
at 
ing 


by reducing the efficiency « of the stilling action, other respects, the 


The schematic drawing, Fig. 31, shows the developed basin and the neal: a 


- 1. Using the design discharge Q, the total head at the valve H, and the hollow- 
_ jet valve discharge coefficient C from Fig. 3, solve the equation Q=CAV2¢gH H 
_ for the valve inlet area A and compute the corresponding diameter d which is 
alsothe nominal valve size, 0 
ca ° Use H/d in Fig. 12 to find D/d and thus D, the ideal depth of tail water = 
‘in the basin. is Determine the elevation of the basin floor, tail water « elevation : Snes 
- minus D. It is permissible to increase or decrease D by as much as 0.4 (d). = 
Use Use H/d in ‘Fig. 14 to find and thus L, the length of of the horizontal 


*. Use H/d in Fig. 15 to find w/a and thus W, the width of the basin fo or { 


 %, Use the H/d ratio to select the proper elas in Figs. 16 and 17 to 
see the model and helpvisualize the prototype performance of the design. The © 
- water surface profile may be scaled from the photograph using the scale on 
the photograph. To convert to prototype dimensions, multiply the scaled edvalues 


a.  Stilling basin dimensions calculated as indicated above are in close agree- 
ment with the dimensions obtained from individual model tests of the basins — 
for Boysen, Falcon, Yellowtail, Trinity, and Navajo Dams, Table 1. Since the 
Boysen and Falcon basins performed satisfactorily during prototype tests, 7 
7 believed that isfactory fu be hydr uli all desi ed 
is belie mat sat ctory future projects aulically 
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= The prototype tests at Boysen | Dam were made » with the - cooperation on of the 
Bureau’s Region 6 office in Billings, Montana, _ and the Yellowstone- Bighorn 
Projects Office, Cody, Wyoming. Bureau personnel at Boysen Dam “operated 
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_ cal Survey personnel at Riverton, Wyoming, made the river gagings, and State — 
of Wyoming made downstream river adjustments to permit above- 


Commission, El Paso, Texas. 
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PREDICTION OF CAVITATI 
A non. number called cavitation on damage number has been 
‘ shock - waves due to the collapse of the cavitation bubbles. | This number gives 
the ratio of the energy absorbed by the material in deformation to the energy of 
the collapse of the bubble. The cavitation damage number has been correlated _ 
with the number of cavities collapsing in a given period. The present concept | 
adopted explains most of the experimental observations made by various re- | 
search workers in this field. (A few of them being the threshold velocity for 
“the inception of damage, incubation period, effect of time, effect of corrosion, = 7 
and so on. i” Experiments show that it is possible to predict the eroded length, — 
breadth, and area if the ¢ geometry | of the cavitating body and the value of the 


cavitation parameter are known. | These results indicate the possibilities . 
an possible cavitation © damage and thus aid in the proper choice of the 
‘material for fabricating any hydraulic machine or a structure. | 


mee 


oft the Journal of the Division, Proceedings of the American So- 
of Civil Engineers, Vol. 87, No. HY 5, September, 1961. 
_ I prof. of Civ. and Hydr. Engrg., Indian Inst. of Science, Bangalore 12, a <4 


2 Senior Research Fellow, Civ. and Hydr. > Sect., Indian Inst. of Science, 
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Although Osbourne Reynolds began experiments on cavitation as earlyas 

* 1894, the method of prediction of cavitation damage has still not been ration- § 7 

i Note.—Discussion open until February 1, 1962. To extend the closingdate one month, 
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‘ “alized. . Because of the se ‘serious interest created by. the severe destruction of _ 


ios, propellers, work was first initiated by the propeller sub-committee in | , 
1915 and t their report | was in In ‘Spite of an impressive ac- 


predicting the of damage, the cumulative damage and the expected 
7 life of a material in a given flow system have not been solved. The work re- _ 
ported in this paper is a start in providing data on which the future progress - 

‘in this field of study m may be made. 
_ Types of Cavitation.—Though there are various classifications available. in 
_ literature, cavitation in hydraulic flow can be grouped broadly into two types: 
(a) cavitation that has its inception in the boundary layer (2) cavitation that has 
_ its inception in the free turbulent shear layer. For example the boundary lay- _ 
er cavitation takes place over a cylindrical body with hemispherical i 
7 [Fig. 1 (a)] , and the free turbulent shear layer cavitation takes place in the 


wake of a circular cylinder (Fig. The results ‘presented pertain to the 


second group; namely, the cavitation occurring in the wake of anobstruction or 

_ a discontinuity called “cavitating body” [Fig. 1 (b)]. er ae) 
a gis Existing Theories of Damage. _—The various theories of damage previously 
* propounded are given in Table 1 in chronological order. ‘9 These theories can 
+: be classified into four main groups: (a) Mechanical action theory, (b) Electro- 
chemical corrosion on theory, (c) and (d) Instan- 


: ory differed in the actual mechanism or proc op by which the material is re- 
arsons an 6 and deHaller7, 8 came to the 
conclusion that the water. hammer pressures were the main reason for voll 
damage. concept is practically the origin of all the other theories. Boet-— 
- _cher9 introduced the idea of surface fatigue and found that the pieces were re- 
moved from surface by fatigue cracks. This idea was followed by the indivi- 
dent formation Thane dents were photographed as early as 1935 
Cc. _Hunsaker. 10 Poulter!1 imagined that liquid particles penetrated 
: _ through the crevices in the material and shattered the material due to high 
‘pressures of collapse of the bubbles, 
” Blectro- Chemical Corrosion Theory.—The electro- chemical corrosion the- 
ag has gained ‘support because of the observed decrease in damage due to. 
~ eathodic protection reported by Petracchi ,12 Nechleba!3 and Wheeler. Put the 
idea of corrosion fatigue explains all these observations. = 
_ Thermodynamic. Action Theory. —The thermodynamic action theory was sug- 7 
pea by ' Wislicenus. g.14 But there is no experimental evidence that the adia- 


Items 3,4 Table 
Items 3,4 Table 1. 
7 Items 4,6 Table 1 

8 Items 4,6 Table 1. 
‘ 9 


“a 11 Table 
12 Item 13,14 Table 1, 
js 13 Items 13, 14 Table 1, 


14 2 
item 12 Table 1, 


— 
| 
Parsons and S. S. Cook, Engineering, Vol. 107, 1919, 
— 


Details of the Mechanism Proposed 


Parsons and Cook 

Rayleigh 
_Ackeret, 


Ackeret and de 
Haller 


Vater, 
Mueller, 


(1948) 
(1947) 


Knapp, R. T 


Poulter 


_ Wislicenus 


(1949) | 


(1955) | Nechleba, 


Plesset and Ellis 


Irving Taylor 


‘ _ Conducted nose cone experiment demonstrating wa- 
hammer pressures, First to propose the me- 
chanical action theory | 


Calculation of the energy of collapse of the bubble 


Showed the possibility of high intensity pressures — 
_ ‘First to suggest the effect of flow velocity on the 
intensity of erosion. Based on this fact he dis- _ 
carded the corrosion theory and accepted the *- 
chanism of cavity collapse as the case is of 
He concluded that in several tests the pressures 
_ developed by the collapse of voids, both calculated | 
and observed, are of too small an order of ade- 
quately to account for cavitation damage 
‘Surface fatigue’ or ‘cavitation fatigue’. First to 
find a correlation between damage and hardness, 
remedy onl 
of used a percussion wave to » produce superfi- 


cial damage on various metals and demonstra- 


| tion of mechanical damage 
Suggested ‘ Periodic strains’ 


First to suggest the concept of threshold velocity 
tor each material below which no damage can be | r 
7 caused to occur, Compares threshold velocity to 

fatigue limit. Suggested that the mechanical _ 
action depends on on the corrosive influence of the 


“nt surface is disrupted local stresses re- 
sulting from the collapse cf vapour pockets, Ac- _ 
cepts an influence of corrosion similar to its in-_ 
= fluence on fatigue limit 
_ Studied the life history of the bubble and found — 
agreement with Rayleigh’s calculations 
= liquid penetration through crevices and 
_ pushing the material out by the liquid mo momentarily 
under high pressure 
Increase in temperature due to adiabatic compres- 
j sion of the bubble and local melting of the mater-_ 
+ Influence of cathodic protection. Electric currents 
due to local mechanical stresses 
of corrosive currents to a thermo- 
_ couple action resulting from local heating of the 
‘metal by the temperature rise of the collapsing _ 
cavity. Supports Poulter’s view 
First to prove plastic deformation by X-ray diffrac - 
“Suggested that at the instant of collapse, minute _ 
i” regions of high intensity ‘water hammer’ pres- — 
de sures are produced, The stored energy ee 
an from the collapse centre as a spherical shock _ 
-wave whose intensity decreases rapidly as its 


_ First to photograph the shock wave radiated from 
_ the bubble by means of Schlieren’s Photographic 


about instentensous chemical | action 


4 1,—THEORIES OF DAMAGE 
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Year nl 

— 

— (1985) | Boetcher, H. M. 

886) 

— (1937) — 
— (1938) 
— 

(1941) | Beeching = 
| 
Petracchi 
(1955) 

(1955) 


1961 
batic c compression of the tthe ts can n produce that order of ncemnenatueve So as 
to melt the materials. Materials such as stone, concrete, and brick get dam- 
Le aged at much lower velocities, in which case ‘the temperature produced | would — 
be very low when compared with the temperature required to melt these mate- 
instantaneous chemical reaction theory proposed by 
“a Theory Adopted.—Theory adopted in the present measnadnene is based on 
the mechanical action theory. _ During the positive pressure phase of the life | 7 
cycle of the cavitation bubble, it collapses due to the condensation of vapor 7 
: ae the bubble as a result of the increase in pressure of the surrounding 
liquid that then flows radially inwards, Near the center, the motion is arrested 
— and lead kinetic energy of the violent ‘hydrodynamic mass is largely 


— — Flowing liquid 


Plastic deform “ton 


“3 
Compression 


4 


(b) FREE TURBULENT ©) STRESS S DIS 
LAYER CAVITATION 


_ 


it is this high pr: pressure 
a shock wave that is mostly responsible for the deformation of the material. The 
dents ‘so so formed were photographed by R. T. . Knapp!6 b by using soft annealed 


od, 17 "The writers have produced similar dents in several materials 
brass, copper, stainless steel, andaluminum. 


al If a bubble produces a shock wave whose intensity is more than the dynamic — 
yield strength of the material, then the material is plastically deformed. It is . 


(16 “Recent Investigations of Cavitation and Cavitation Damage,” by R Knapp, 
~ Transactions, ASME, Vol. 77,1955, p.1050, 
17 "The Formation of Pressure Waves by Cavitation,” by W. Guth, Proceedings, 
ymposium on Cavitation in N.P.L. 1, 1955. 
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CAVITATION DAMAGE 
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a (b) ON ALU 
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| (c) ¢ 


FIG. 3.—DAMAGE BEHIND CIRCULAR CYLINDER 
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displaced sideways towards the free surface, undergoing relatively severe dis- =e 
tortion to produce the raised of indentation 2). a 


ad 
-BULGING OF THE A SHOT- PEENED PLATE 


OF DAMAGE 
a 


@) FORMATION oF DEEP CRATERS 


a 
bubbles are pits are formed and the sur- 
face becomes spongy [Figs. 3 (a), 3 (b) and 3 ( (c)]. The test bulges outas 


‘in the of shot- ~peening (Fig. 
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‘CAVITATION DAMAGE 


9 By m micro- soneiiieds it is ; possible. to detect the distorted metal and oe 
cracks, Figs. 5 (a), 5 (b), 5 vin and 5 (d ) show the oa view of the test piece > 


UNDAMAGED cross SECTION 


ELOPMENT OF SLIP LINES 
TION OF CRACKS 


 (d) CHIPPING OFF OF METAL 


6.—CROSS SECTION OF DAMAGED 


_ shows the undamaged surface, and Figs. 5 (b), 5 (c) and 5 (d) show how the 
lastic deformation develops into big craters. _ Because the shock dnnesaneeglll 


we 


if 
— 
4 
— 
— 
= and 6 (e) show the cross- 
(b), 6 (c), 6 (d) and 6 (e) s 
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tion of the test 84 tir times at stages 

_ Of the deformation, Figs. 6 (b) and 6 (c) show how cracks are formed due to 
repeated impacts. — Figs. ¢ 6 (d) and 6 (e) show how material gets chipped off __ 
_from— the surface finally. _ Fig. 7 shows the spherical dents produced by the a 

shock waves going through towards the backside the The has 


been taken from the backside of thet 
i 


FIG. — BACKSIDE OF PLATE SHOWING DENTS 
- Notation. —The letter “symbols it for use in this paper are defined 


P her they first appear, in the illustrations or in the text, and are — 


for convenience of reference, in the Appendix. > 
CAVITATION DAMAGE NUMBER 
‘i one imagines the shock waves as solid balls’ dynamically impinging over. 
= the solid surface, just as in the case of the dynamic hardness tests, the ee 
deformed is proportional to the energy creating the indentation. "The energy 
creating the indentation is proportional to the energy of collapse of the bubble. 4 
= Hence the volume deformed, the dynamic c yield s stress of the material, , and t the 


energy of collapse are related by 

in which Y is the dynamic yield stress “ the e material, Av denotes the volume 

_ deformed per collapse of the bubble and E, refers to the energy of collapse of 

| 
by 


— 
— ff. 
= 
— 
— 
al 
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in Po is the of the liquid the the 
moment the bubble begins collapsing and R, is the maximum possible radius" 7 

of the bubble at the beginning of the collapse. The value of R, can be esti- 
mated from the equation 


in which Pmin is the minimum pressure the liquid can withstand (that is, the ¥. 
_ tensile strength of the liquid), gis the surface tension of ms liquid and Py , is — 
the of the liquid, 


corrected for the partial pressures of gases in the liquid. Seen sore 
7 Now it is necessary to estimate the number of bubbles per cavity. Each 


= would occupy an area equal to 1 sol _ If the area of the cavny is A, then 


hee will be a function of the gas content of the liquid and Py will have to be a 


volume deformed per bubble is ‘Av, 1 and 2, 
Po Ra 


Therefore, ‘the total volume V deformed | per cavity is 


: 


. .(8) ‘aa 
v/Ac ‘is the average depth of erosion if the area of ereaten Ae is 
equal to the. area of the — A. (Fig. 18); thatis — 


A. 


that is, vis the v volume deformed per per cavity, Re is s the : area ad erosion ares : 
is equal to the area of the cavity), and iis the average depth of erosion. 


Eq. 8 yields the ratio of the energy y absorbed by ‘the ‘material to the energy 
of ‘collapse of the bubble for one cycle. If there are N cycles, the cumulative - 

energy a absorbed by the material will be equal to Y » i. . The total energy of 


|. 
— 
q 
“4 
2S 
is iim 
il 
— 


will to N Po Ry Hence, ‘the: ratio - 
»¥ sorbed by the material for the deformation to the total energy te by a 
- cavities in doing w work after N — is given by 


Eq. 10° gives : a non-dimensional number known as the “cavitation n damage num- 


“uly developed damage zone 


= 


aw 


(1) Threshold zone Re 


ese exact relationship between the cumulative absorption of energy by: t 
‘material and the number of cycles of dynamic indentation is an unsolved prob- 
lem in plasticity. If the mechanism of repeated dynamic indentation is as- 


a 
7 
lz: 
— FIG. 8.—E-N DIAGRAM 
m 


~@ 


_ sumed to be similar to the Charpy’s impact fatigue, +° 18 then the E - N diagram | 
in which E is the energy of indentation and N is the number of a, 
(Fig. 8) can be divided into three zones. 
_ Zone 1.- —When the energy of indentation E is below the fatigue limit of i 
_ dentation (Ef) which will produce no plastic deformation for infinite number of | 
_ cycles, then there is no ® Gamage at all. This is called the threshold zone 
>i - Zone 2.—When the energy of indentation is above the fatigue limit (Es) but 
below the yield limit (Ey) (the yield limit being the minimum energy of inden- 
- tation for which the plastic deformation takes place at the first cycle itself) _ 
_ then the damage will start after a given number of cycles depending on the en- 
level of the indentation. This zone is called the zone, The 


\VITATION DAMAGE 7 


byld 
in which m is a constant depending on the material. =e 
_ Zone 3,—When the energy of indentation is more than the j yield e energy, y, then =, 
the damage starts from the first cycle. This is called the fully developed dam- 
age zone. Inthe present paper only the third zone _ be hated in | detail. 


ia the number of emiin ‘through which the indentation has undergone is 


equal to the number of cavities produced by the flow. In 1935, Hunsaker!0 
= - found that the Strouhal number was a constant if the length of the « of the cavity ! is taken 


cavity, and Vo is the velocity of flow. In 1955, K. K Shalnav20 reported 
for a given length of the cavity the frequency was directly proportional to the 


: in which f is the frequency « of shedding of the cavities, lis the length of the 
_ velocity. The writers, “using Knapp’s data21 found that Hunsaker’s finding 


lationship is true irrespective of the geometry, and the constant of pro- 
_ portionality will change with geometry. 


“Experimental of the Intensity of Erosion Du Due to. Cavitation,” by K.K. ‘Shal- 
nav, Proceedings, Symposium on Cavitati 1 in Hydrodynamics, N.P.L. London, 1955, — 


was true even for axis- symmetric cavities, Hence, it is assumed that this re- _ 


pat 21 “Further Studies of the Mechanics and Damage Potential of Fixed Type of Cavi- 
‘ties, by R. T. Knapp, on in Hydrodynamics, N.P.L. London, 


1955, p. 19-1. 
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Eq. 12 can 


and l/dis s the ‘relative length of the. 


‘4 


_ Experimental investigations we were conducted to determine ome elie be- 

tween the re ene number r Cp and and the number of cycles of indenta- 

EXPERIMENTAL SET- UP AND P PROCEDURE URE 


speeds by means of a variable speed motor (Fig. 10). The test liquid was cir- = 

culated through the chamber by means of a pump. The pressure and the quan- — 
i tity of the liquid in circulation could be controlled by means of two valves ve 

upstream and one of the chamber. By controlling the quantity of 


which | a eiroular ‘metallic disc was rotated at 


disc was made from two 1/16 in. metal plates fixed to “each other. ¥ Several 

»* holes of various diameters were made through the disc and the test specimens 

_ ‘were fixed behind the circular and rectangular holes of various sizes. = — 

Two-dimensional Open Circuit Water Tunnel.—The two-dimensional open 

water tunnel consisted of a 100 turbine pump, pressure ‘regulating 


ia 22 “Some Experiments on Cavitation Erosion in Water Mixed with Air,” by R. ~ > a 
_ Rasmussen, Symposium on Cavitation in Hydrodynamics, N.P.L. London, 1955, p. 20-1. _ 
28! “Study of Corrosion and Cavitation Erosion Damage,” by J. Z. Lichtman, D. — 
‘Kallas, C. K. Chattan, and E. P. ‘ASME, Vol. 80, 1958, p. 1325. 


(4) 
This relationship is shown in Fig. 9 plotted from Knapp’s data.220 00 
24 ___Now the number of cycles of indentation will be equal to the number of cay sy i 
Fee eee ities that have passed in a given time t. Hence, _ i if 
| 
— 
— 
| rotation of the liquid within the chamber aiong with the rotating disc was re- 
ae duced by providing baffle plates on both sides of the rotating disc. A total - 
ar clearance of 5/8 in wac lowed the diectnon rotat hetween the ffl e Tho | 
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O Stroboscope 


K Pressure gage 
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Baffles” 


E Rotating disc 
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001. 
7 oma manifold, , bypass needle valve, ‘contraction cone, test se section, ,and diffuser. The 
fa test section was 4 in, by 3/4 in. by 13 in. The maximum velocity attainable at *: 
rad the test section was 100 fps: with a static pressure of 200 ft . Fig. 14 shows the 
_ layout. A circular cylinder was fixed in the test eet and the damage | shown 


Ye ~ Contraction 


meter 


To Pelton 


.- Figs. 3 (a), 3 (b), was reproduced in the \ wake of the ner sac de- 
tails regarding this tunnel are available. 4a 


Prediction of Damaged Area. -Firstly it is necessary to predict the length, 

the breadth, and the area of damage so that that particular area can be pro- 

tected bya stronger material. Further the area of damage is required to com- 

_ pute the total volume loss from the cavitation damage number, Therefore, ex- 
periments were conducted with a circular cylinder as the cavitating body that 
produced cavitation when it was kept in a . stream of flow. _ These experiments 

_ were conducted in the two-dimensional open-circuit water-tunnel described = 

previously. The geometry of the flow behind the circular cylinder can be con- a 

as being of three regions at high Reynolds numbers. a thin layer 


‘lig 


« “Cavitation and Cavitation Damage, ” by A. ‘Thiruvengadam, thesis to 


“ the Indian Inst. of Science, in Bangalore, India, in 1959, 1358, in partial fulfillment « the re- 


quirements for the degree of Master of Science. 
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at the front and sides of the body the flow is irrotational; behind the body | 


there will be the third region of free turbulent shear flow or wake flow. The 
Pex region consists of two zones divided by the free streamlines (Figs. 12 and 


13). The first zone is a region 1 of high shear due tc to resing nananeioes and it is it is 


- Free stream line 


extremely t turbulent, The second : sone is a region of low dis- 


turbance. The s shear —_ springing from both the sides of the cylinder 


— 
— 
— 
— 4 
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body . In the case of high velocity flow of a liquid, the second z zone becomes © ie, 
gaseous filled with the vapor of the liquid. The inner boundary of the shear 
layer becomes the free-streamline of the classical theory (the condition of 
free- -streamline being that the velocity and, should be co constant a 
the method of transformation coupled with the free-streamline 
theory, Riabouchinsky, Birkhoff2° and various other workers have tried toeval- 
uate the boundary of the wake cavity. Such a procedure deals only with an ideal 
hj — liquid. They cannot be applied to a real liquid that is in a highly turbulent two 7 
: phase flow. The classical theories deal merely with the overall forces or ; 
changes in pressure rather than the local details of the flow pattern. i, In fact, 
the real conditions of flow are quite contradictory to the free-streamline the- — 
eee at the junction of the boundary of the wake and the irrotational flow cover- | 
_ ing it. Because | of the high shear along the unguided surfaces of discontinuity vee 
2 at the boundary, the wake, which is a large cavity itself becomes the seat of 
intense turbulence that diffuses laterally as it is carried downstream (espe- 4 
at f flows of high Reynolds numbers). This can seen through a stro- 
7 boscope. The e cavitation ‘damage takes place only at the boundary of the cavity. b> 
_ Each vortex in the shear layer is entrained with a cloud of vapor bubbles, and \ : 
the diffusion of these -filled vortices into the causes: the 


Fig. 13). ‘Fig. 12 is an actual plot of this shear omen in relation to the circular 
os cylinder (full scale), Fig. 14 shows the variation of the width of ; erosion with © 
_ the longitudinal distance. From this it can be concluded that the width of ero- 4 5 
sion follows the pattern « of vortex diffusion in the turbulent shear layer. sn 7 
as Experiments were conducted to study the manner : in which the eroded di- id 
mensions vary with cavitation parameter, The rotating disc apparatus (Fig. 7 
10) was used for this purpose. The velocity, pressure and size of the hole 
were varied independently in the rc rotating disc apparatus. . Figs. 15 (a), 16 (a) _ 
_and 17 (a) show how the relative length of aenenne, relative breadth of erosion, — 
_ and the area of erosion vary with the cavitation parameter for circular holes 
1.27 cm, 0. 95 cm, 0.635 cm. . Figs. (b ), 16 (b) and 17 (b) show the 
same relationship for square holes of 1.27 cm, 0.95 cm, and 0.625 cm. These 
- results indicate that the eroded area and the relative length of the cavity de- 
pend only on the cavitation ‘parameter : and the geometry | of the flow system. 


4 


The area of erosion is a function of time because of the cumulative nature 
8 of the damage. ‘The area increases with time until it becomes equal to thearea 
of the cavity itself. This relationship is plotted in Fig. 18 for two values of K. 


= Wakes, Jets and Cavities, by G. Birkoff and E. AL. Zarantonello, , Aca ademic Press, 
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nik: 7 merge with each other after some diameters from the cylinder, The velocity | — 
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OF DAMAGED AREA WITH LONGITUDINAL | 


Diameter 

0.635cm 
© 0.950cm 
e 1.270c¢ 
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ine the relationship bewteen the cavitation damage number Cp and the n num- 
er of cycles of stress ee. _ The cumulative increase in the westadl 


(b) SQUARE HOLE 


lo -0.950cm 
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FIG. 17. LATIVE AREA VERSUS 
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a By depth of erosion > i was computed by finding out the weight loss and dij 


= 


DA MAGE. 


in which y aw Aw is the weight loss, py, denotes the density of the material, and 


Ae refers to the final area of erosion. The number « of stress applications was 4 
taken to be equal to the number of cavities given by Eq. 15. | The pressure at 
which the bubbles start collapsing was assumed to be proportional to the free € 
stream pressure maximum radius Ro of the bubbles was assumed 
_be constant for a given temperature, ~The value of Cy was calculated, based on 
these assumptions. As the material used was 99% pure aluminum, the value of 

the yield stress as given by the manufacturer was taken to be equal to 786 | kg 

per sq cm, _ The pressure was © varied from 0.351 kg per sq cm to 1.05 kg per 

sq cm, and t the velocity was varied from 34.6 m per sec to 46.6 m per sec, 


‘FAG. 18. EFFECT OF TIME AREA OF DAMAGE 


intervals. For each pressure and duration of tests varied 
; from 2 hr to 3 hr, Fig. 19 (a) oo the relation between | Cp and N for r square — 
a. holes. The variables were Vo Po, Ant and d. There seems to be an excellent 
. correlation between Cp and N. . Fig. 19 (b) shows Cp as a function N for cir- _ 
cular holes. Fig. 19 (c) shows the results for a circular cylinder kept in the 
7 two- -dimensional open circuit water tunnel. . In this case time was the only 
- variable. These experimental results support the theory of the / shock wave in- 
7 dentation and the concept of the cavitation damage number. In addition to this, 
this correlation is useful for the prediction of the cumulative damage « caused 
ms The phenomenon of cavitation damage involves three major groups of var- 
: iables: (1) the properties of the materiel; (2) the properties of the liquid; and 
Experiments were conducted in the two-dimensional open circuit tunnel to 
7 . investigate the effect of the mechanical properties of the material. Five sheets 


of pure ( aluminum with different ‘mechanical proper 
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were procured from the Indian. Aluminum Company and tested. . The results - 
are shown in Figs, 20 (a), 20 (b), 20 (c) and 20 (d). It can be seen that )'i is 
_ inversely proportional to all the mechanical properties because each mechan- _—E_s| oo 
ical property is proportional to the other. The same chemical composition — a ae 
was used to ensure that there was no corrosion otherwise. Another important e. a 
point is that none of the five sheets underwent any static corrosion in the pa 
months they were used for tests inthe laboratory water, = 


Ph... the case of brittle materials like | concrete, brick, and stone the damage ' 
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“Fatigue str = 


It loss, in oz/ hr 


oz 


ht loss, i 


veig 
Nm 


Camere stress, in 
thousands of Ib /sq in. 5 


corrosion the ‘the m property may be used to represent 
property of the ‘material. This" explains the inclusion of yield stress” in in the 
_ The properties of the liquid had not been varied in this investigation. | The 
_ liquid tested was laboratory water at 32°C. The temperature was that of the — 
flowing water. But the concept of cavitation damage number still takes into 
account the physical properties of the liquid such as surface tension, vapor ’ 
pressure, and the tensile strength of the liquid that, again, depends on the gas a 
content of the liquid. — ‘It is assumed that Ro is is constant as long as these pro- a 
perties a are the same. The validity of this assumption is yet to be established. © 
 -It is now necessary to investigate the manner in which the ——— 


of the flow such as the geometry of the flow, pressure, , velocity, and time af- 
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is a function of the tt rh of flow, time, the characteristic linear dimension, 7 
a the relative length of the cavity. had The relative length of the cavity again 
_ depends on the pressure and vg of flow. Finally, the energy level of col- 


| tren liquids, and for different flow geometries then the extent of likely 
je damage in any prototype can be confidently predicted. Just asa 
_ designer of a hydrodynamic ayehem can get an idea of the drag forces if the 
coefficient of drag is } Known, or of the lift forces if the coefficient of lift is — 


the cavitation — number. Th¢ ‘The following a will illus illustrate the meth- 


‘Determine the cumulative damage caused inl hr by a enna be of 2. 5 


cm diameter rotating with a velocity of 45m per ‘sec and with a free s stream _ 
pressure of 0.75 kg per sq cm. _ The material is aluminum with a yield stress 
of 786 kg per sq cm. | The liquid is water at 32° oe 


Cavitation parameter: 


45) (3600) (100) 
(200) . 2.34 x 10 
2.50) (2.7 2.34 x 
ueo 


ae a the experiments, and the correlation seems to be satisfactory in the fully de- > a _ 
— 
— 
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Po 


Cumulative volume: 


> 


10-2 x 16. 


hus the cumulative weight loss equ: equals: (5 yy) @a) - 
Eatin = (7.0 06 x 2.7 x 10-2 g) (2.7 gm per cu cm) 


* 
in which Pe is the density of atuminum 


is method: would be very useful if standard charts could be prepared for 

: - the guidance of designers. The generalization of this method is possible when 
a better under rstanding c can be had of the aspects presented in the following sec- 

Concept of Inception of Damage and Incubation Period. _In Fig. g. 8, the E - _N 
diagram was divided into three zones. The work analyzed here belongs to the ; 

& zone, , namely the fully developed damage zone, For the ‘first zone, the — 
i ete ses zone, it is noticed that there is a lower limit for the velocity of pe 

: ~ named as the threshold velocity below which the material will not be damaged 

at all. In any flow system the cavitation starts when the value of cavitation 


par rameter K reaches a critical value K;. ‘Thus, ; : 


- Py 
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given Therefore, for a given the maximum possible 
& pressure is p;, above which there is no cavitation at all. Therefore, the maxi- 
mum energy of collapse for that velocity is proportional to p; ‘R3. "This value © 
m4 of pj R? should exceed the fatigue limit of the material in repeated dynamic 
indentation. But it is known that the impact fatigue and bending fatigue are 
directly related linearly. 18 5 Similarly, if it is assumed that the indentation fa- 


= (or surface fatigue) and the bending fatigue are directly related, then the 7 
bending fatigue limit can be used for the correlation. 


‘Then the stress produced by the shock wave wit be to the val-— 


gerne 


(22) 
pK 


avitating 


Thus, ons any cavitating Rew ayem,, ‘there ‘is a threshold velocity given by 
Eq. 23 which is a function of the fatigue limit of the material, the cavitation _ 
inoagtion parameter, the vapor pressure, and density of the fluid. ee 
_ For zone 2, the incubation zone, certain number of cycles are required be- yy 
fore the plastic yielding starts. ‘This value of the number of cycles is given by 
Eq. 11. Investigations are in progress (as of 1960) to analyze these two zones. 
Effect of Corrosive Environment.—It is well- known that a corrosive envi- 
ronment reduces the fatigue limit considerably. Hence, the threshold velocity. 
is also lowered. ° This e explains the observed increase in damage when corro- 
P sive liquids are used. This is the reason cathodic protection reduces cavita- 
tion damage. Experiments are being conducted to investigate the effect of cor- 
ore liquids and cathodic protection on cavitation damage —_— —— 


: we For a given flow | geometry, if the value of the cavitation parameter is 


given, then the length of erosion, breadth of erosion, anal the area of erosion» 
can be by mo model analysis. 
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CAVITATION DAMAGE 
if the mechanism of damage is: assumed as the dynamic indentation pro- 
_ duced by the short range shock waves arising from the collapse of vapor bub- - 
bles, then a useful non-dimensional number known as cavitation damage num- 
oe ber is obtained. This number is the ratio of the energy absorbed by the mate-— 
rial in deformation to the energy of collapse of the bubbles. 
a _ It is assumed that the number of repetition of the dynamic indentations > 
E 8 s equal to the number of cavities passed during a given interval, and the cavi-— . 
_ tation damage number has been correlated with the number of cavities. er an 
cum The trend of the experimental correlation can be explained from the 
cumulative damage hypothesis used in connection with fatigue failures, | a 
material deforms elastically as the number of cycles are increased toa de- 
fined number and, hence, it does not absorb any energy in the elastic recovery _ 
stage. However, when the plastic deformation begins the material absorbs 
- more and more energy and eventually reaches a stage of viscous flow. At the — 
— same time, the attenuation of shock waves takes place because of the damping — 


damping also the energy of collapse transmitted 
material defor mation decreases. This the 


‘The fol following use in this paper, essentially 
7 _ with “American Standard Letter Symbols for Hydraulics” (ASA Z10.2-1942), = 
"i prepared by a ‘committee of the American Standards Association with Society © a 
representation, and approved by the Association in 1942; 


= = cavitation number; 
Strouhal 


= energy ergy of indentation; 7 
= energy of collapse f each bubble; 
= = fatigue limit indentation; 


Ey = yield limit; 


= average depth of e erosion 


= length of — 
_ = number of indentation; 
= free stream pressure at =m cavitation starts fora a given nvelocity \ ey 


ai 
q 
| 
K = CaVitation parameter q 
— ‘ 
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=n minimum pressure can withstand; 


= time; 


deformed per collapse, of bubble 

Aw = weight loss; 


: relative 1 length of the cavity, 1/4; a 


ible at themoment the bub- 
Zz "Vapor pressure Of 
Ry = maximum possible radius of the bubble at the beginning of the collapse; 
— 

$ 


-HYDRAUL ICS DIVISION” 


a) Proceedings of the American Society of Civil Engineers 


VIBRATION OF GATES DURING OVERFLOW AND 


By Eduard Naudascher,! 


and the basic mechanism of ‘excitation are 


a 


ond une. in low dams or weirs iaiibias under r large depths of overflow in 
combination with underflow may be subject toa violent vibration brought about 
: = a periodically changing flow patternonthe downstream side. Because of the 


rare occurrence of such operating conditions, little attention was given to this 
a few observations 


a written request must be filed with the Executive Secretary, ASCE. "aus paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 

ciety of Civil Engineers, Vol. 87, No. HY 5, September, 1961. | 

al 1 Asst. Prof. and Research Engrg., Iowa Inst. of Hydr. Research, State Univ. of Iowa, 

Towa City, Iowa; on the Institut fur Hydromechanik, Technische 

nik 


= 
_ Violent vibrations may occur at gates being subjected simultaneously to 
overflow and underflow. Model studies have been made to evaluate frequency 
aprons. 
— 
ug 
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have been reported on on. 2, 3,4,5,6,7,8,9 In 1953 Petrikat3 described model ite 7 
of a tainter gate in which the entire test apparatus was destroyed due to 
rt _ vibration. The gate, designed for a tropical river, was operated with simulta- 
neous overflow and underflow Goring flood periods when the head oo, ex- 


more section flap—and any type of Pro leaf spillway gates. With such gates the — 
a flood waters will at first be released by lowering the flap or the upper gate 
leaf. Later, in order not to raise the pool level unduly, both gate sections are 
lifted. During this course of operation, water is discharged above and below 
the gate and vibrations are usually inevitable. The same vibration phenomenon ~ 
8: also encountered during the handling of stop logs in running water. 7,8 ——- 
‘This paper is an abstract of an extensive model investigation conducted at 
the Hydromechanics Laboratory, Karlsruhe Institute of Technology, Germany. - 
Its purpose is to make available the most important aspects of the problem as — 
to criteria for preferable types of gates and aprons, ‘magnitude and frequency | 
of the exciting forces, and means of mitigating the vibration. Greater use has 
been made of dimensionless representation than in the original report, and 
_ emphasis is placed on the explanation of the exciting mechanism of the vibra- 
. tion. A complete description of the investigation is to be found in the writer’ s 
SIMILARITY CONDITIONS 


“In ge In general, , the vibration of a gate with overflow and underflow has to be 


considered as a self-excited vibration with four coupling members | (Fig. 1): 
1. Upper jet (referred to also nappe) 
Space between jets 
4, Elastic gate 
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Webrverschiuessen, "by E E. Naudascher, Techn. Mitt. Krupp, 
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OF F DOUBLE- _LE AF GATE IN A 50- CM GLASS- 
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of flexure and torsion. In order to simulate the dynamic behavior of the gate as 
well as the dynamic behavior of the flow, at least two model laws have to be 
: satisfied: that of Cauchy (for structural elasticity) and that of Froude. This is Ee 
impracticable in most cases; moreover, it would restrict the applicability of 
the test results to a few specific cases. To overcome this difficulty, the model 
P gate was made very stiff and retained by fixed, in contrast to elastic, supports. 2s 


Buc Under these conditions, the Froude conversion laws canbe applied, that . 
for 


which is ratio of the in the model to that in ‘the guetotype. The 
of the four coupling 1 members has te been eliminated by this sim-— 


asa first approximation can be computed as a forced vibration, when 
re excitation due to the unsteady flow system is known from model experiments. a 
- the | case of one degree of freedom, for instance, that in the direction of 
the prototype gate , suspension—the forced vibration of the gate is described by 


in which m —; = inertia force; ¢ x= = elastic force; = viscous damping 

= turbulent damping force, + = Coulomb’s friction = 


gam opposing motion; Pe cos (wt) = exciting force in x- te a to the 

a flow; and x = deflection of the gate from equilibrium position; t = time; m = 
mass of the gate and added mass of the water; c = spring modulus in x- "direc- — 
tion; k, and k, = ad damping coefficients; and w= angular of the excita- 

‘The error in this simplification will become sa as the stiffness and 
the damping of the prototype gate | become larger. The permissibility « of the | 
discussed simplification was checked by exploratory tests witha flexible sus-— 
model gate which are tobe described later. Final conclusions may only 

be obtained by prototype menourements, 


MODEL | OF A DOUBLE. GATE 

. 2) was based On the design of a a double leaf gate in the River Ruhr ee 


As far as the flow is governed mainly by the _ 
inertia and gravity forces of the moving liquid. The complexity of the problem | ‘ — 
Be. __is implied by a question as to how far the elasticity and damping characteristics [i 
— 
ce 
al 4 
— 
— 
— 
| — 
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_~ type resulta existed. The lower leaf of this gate is of the girder-plate ep, 
_ the underside of which was provided with 15 holes, 4 mm in diameter (model — 
dimensions). About half of the clear width was represented in the model . Three 
: ea _ types of lower leaf were tested in combination with five types of apron. Al-— 
though performed on adouble leaf gate, the model tests are also applicable toa 
gate with» an attached flap in the lowered position. A pipe of 1. 2-cm inner 
% diameter served for aeration a of the nappe. 
_ Flow Regimes.—An important quantity defining the flow conditions near the 
Gate is of the wa water si surface beneath the crest, This “inner” 
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. MODEL OF GATE 


e IN COMBINATION WITH APRON TYPE 1 (DIMENSIONS IN ten 
on level depends on the a ‘oft upper nine lower j jet, on the t type ot 
apron, and on the ‘size of the lateral air ‘vent aerating the -happe (Fig. 5). 


: upper limit for dj. Because the edge is irregular (Fig. 3), this limit does Bs 
Rot appear clearly defined in Fig. 5. The lower limit for dj is reached when the Ma 

4 tailwater drops below a certain level leaving dy—and, hence, the discharge un- 
affected by tailwater conditions. Analogous to the terminology for sluice gates, | 

: the corresponding discharge will be called “free” despite the fact that the sis 


— a, ‘Type I ‘Type — 
— 
— 
— _ 
— 


st 


cres 


upper leaf 


tion 


| 


on 


ti 


49 no aera 
with aera 


7 i= 


4 


IN cm) © 


a 


é 


—_— 8356-4 
fi 
APRON TYPES TESTED IN COMBINATION WITH 


rizonta 


he 

32 


MENSIONS 


= 


_ 


LE 

| 


(I 


be 


devez 


m 


the surface roller remains on the water 
da: 


three flow regimes may thus be e distinguishe 


7 1. * The ‘flow with: the nappe remaining on the water surface 


2) keeping the gate supported on a set of adjusting screws 0.1 mm aprat. Only 
oe the force exceeded the weight set by the balance was the model shortly _ 


_ displaced from its supports. In this way the extremes of the hydrodynamic — q 


load that were reached or exceeded about once a second were measured, 
This “concept was adopted | in order to exclude peak values due to sporadic 
shock loads. A correction was made for the difference in buoyancy of the _ 
__ lower leaf. A typical plot of these downward or upward forces against the head ay 
differential between | headwater and tailwater is represented in Fig. 6. The q 
ss difference A F between upper and lower carey indicates the force fluctuation _ 
that is responsible for the gate vibration, 
_ Excitation.— —A measure of the sniguitets of the excitation may be set forth 


of form of a coefficient, similar to the lift coefficient of airfoils, by dividing 
the amplitude of the fluctuating by the product of the 
(00 of the lower leaf A A and the dynamic pressure. Sub- 
ll - d2) for the dynamic pressure, the excitation coefficient is 


‘It must be kept in mind that AF does not a the pe: peak- — ene but, Par : 
ae a double amplitude that is exceeded about once a vos 
ee coefficients for the investigated types of lower leaf (Fig. 3) are 
" - shown in Fig. 4. In combination with a apron type 1 (Fig. 4) the leaf type 0 was 
found to be most favorable. In combination with lower leaf type I, the excitation 
— coefficient was a minimum for apron type 1 and a maximum for ahorizontal _— 
floor ora downstream slope of 1:1 (Fig. 8). With increasing gate opening the : 
elation Increases without noticeable change in frequency. 
= - Influence of the Geometry of Gate and Apron. —The most violent pressure _ = 
» ‘ fluctuations are to be found in the vicinity | of the juncture of upper and lower — 
ele jets. The more the exposed areas of the gate are withdrawn from this zone, 
i the less will be the transmission of pressure pulsation to the gate. To withdraw ig 
_ the gate structure from this zone, the floor downstream from the gate as well — 
as the bottom of the gate have to be inclined as much as possible. A sloped 7. 
floor with a depressed stilling basin was found to be most effective, because 


- bottom and the surface of the lower jet becomes very small, as shown in Fig. 
9(c), the sarees of vortices from this surface and, hence, _ the excitation will 5 


— — 
— 
— 
4 Hydrodynamic Forces on the Lower Leaf.—The hvdrodvnamic forces acting _ 
; 
| > 
“¥ 
& 
: b ra _ this permits a larger range of dispersion of the nappe (Fig. 9). This means + i 
increase in damping for equalheaddifferentials4d. 
ks 
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be disturbed. ‘This phenomenon can be utilized to mitigate 11 12The 
decreased excitation coefficients for the bottom plate type b, as compared with 
type a, are . due to the same cause (Fig. 21). The other exception from the rule | 
f is an excessively steep floor that leads to a separation of the lower jet. For Bees 
ae small gate openings, flows with rising and omens jet are equally possible at —? 


turbance by the » roller occurs, ‘the separated limiting streamlines of 
the lower jet tend to cling either to the floor or to the gate bottom, provided - 
the the gate is ‘not 0 of the open- -truss type. Although the generation | of vortices was — 
prevented: in the case the rising jet, both unfavorable negative 


4 mid position independent of the tailwater. The influence of the crest shape was 
_ investigated with the second test arrangement (Fig. 14), 

Influence of the Aeration. —Without aeration of the nappe the tendency toward 

vibration was greater in every case e investigated. As will be explained in the ary 

- last section, this is due to increased pressure change beneath the nappe. It is 

: _ interesting to note that there was a well-defined maximum of the exciting 

< forces for nonaerated flow at a certain head differential that was independent 
of the gate opening. From Fig. 7 the head differential is found to be approxi- 
= mately do - d2 = 16 cmor the ratio dg/do = 0.57. According to Fig. 5 this ratio : 
a corresponds to the condition at which the inner water level reaches the down- . 


PS stream edge of the gate crest (point A in Fig. 1), that is, when the air space 
_ below the nappe just means of a momentum consideration 


= had no effect within the range of small head differentials do - dg when 
Influence of the Elasticity.—As a supplement to the series of tests with the 
gate in fixed positions, measurements were conducted with an elastic gate 
a suspension, By means of a spring suspension system the lower leaf was free 
-s to vibrate vertically. The upper leaf remained fixed. In the’ following only the 
_ The frequency of the gate vibration is unaffected by the natural gate fre- 
E ae within the accuracy of measurement, as verified by Fig. 11. By plotting 
the frequency against the difference dj - -y, which for a constant do isa meas- — 
os an ure of the velocity at the surface of the lower jet as well as of the mass of 
water included between the jets, a fairly good correlation was obtained, Fur- he 
= evidence that the “frequency of the excitation depdns on these quantities 


explained by the random pulsation of the surface roller and by the 
fener water surface. Both sources of disturbance grow with dropping tailwater — ae 
the inner water surface approaches the gate bottom, only 
periodic shock “loads” are induced by waves (Fig. 16(a) and Fig. -17(a)). a 
shown in the last section, the gate is self-controlled, An 


ait control may apparently be effected by the elasticity of the gate in 


‘Die schwingungsanfachenden Kraefte im Wehrbau,” by K. Petrikat, Man- 
12 “Beitrag zur Untersuchung der schwingungserregenden Kraefte an "gleichzeitig 
eber - und unterstroemten Wehrverschluessen, ” by E. Naudascher, Tech. Mitt. Krupp, 


17, No. 5, 1959, p. 270. 
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_b) Without overflow: h= 0 24.3.cm; 


gate I with heavy 
(in still water) — 
netural Frequency of | 
gate I with light spring _ 
‘(in still water) 


gate I with light 
(in still water) 
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gate type I 


"APRON TYPE 1, do = 36.8 cm, h = 


— 
3 
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that the gate lip will be per sodicalty deflected. om lower leaf flexibly sus- 
E pended, it was possible to magnify the gate deflections appreciably by knocking 
e - lightly on the flume wall with the frequency of the alternately generated vor- 
7 tices (Fig. 12). The reason seems to bean ¢ augmentation of control rather than 
i: an extraneous excitation, because the external impulses had no component in 
‘. “4 the direction of freedom of motion ofthe gate. Of course, an amplification of a 
Be. multiply coupled vibration like the one under consideration depends strongly i.e 
r Ro on the tuning of the individual coupling members that might have been rather 
Ry favorable in these tests. It is planned to check this phenomenon of extraneous 
ae control by field tests because of its possible importance. For instance, a simi- { 


effect may be induced by ship engines in navigable rivers. 
Means of Mitigating Vibration. —Criteria for favorable types of gate and 
apron have been given. In case the remaining excitation would still affect the 


safety of the structure, additional means acai tas vibration have to be con- 


seme been mentioned in connection with Fig. 9. Extensive studies with the 
Second test t arrangement —- 14) indicated that one or two > horizontal plates, 


: ; . gravel and debris. Therefore, it is more advisable to use two lateral ledges 
ae attached to the piers and to slope the floor downstream from the gate at 1 on 2 
- - (but not steeper than 1 on 1.8). By deflecting the upper jet these ledges give 
He rise to eddies along a vertical axis that disturb and break up the exciting vor- | 
tices generated along a horizontal axis. Most effective was the arrangement 
_ depicted in Fig. 13, because it provided additional aeration. None of the flow _ 
spoilers attached to the crest, which were successfully used to cure fluctuating > 


‘a wake either through a controlled cpentag | in the skin plate or by deflecting 
parts of the upper jet back toward oe. had device 
‘ 


In the first test arrangement there were no means to separate wiedenee - - 


4 in the fluctuating hydrodynamic force and shock loads from the main a 


_ precise relationship with time (of the pressure and force 
fae was studied with an arrangement of an idealized gate the base of Fo 
3 which was simply represented by a bottom plate (Fig. 14). yan Pee aia 
Depending on the type of gate, — 
‘distinguished (Fig. 15). 
three zones" which the hydrodynamic forces are transmitted to the 
= structure. The loading of an upper leaf or a flap, for instance, is signified by | + - 
« ‘Fig. 15, III. The resultant exciting force canbe found by superposition for any : 
(of “these cases when the pressure amplitudes within the regions a, b, c mt 
their temporal correlation is known. Hence, it was decided to record simulta- | 
== the pressure fluctuation at afew representative points in these regions, 
+; _ In addition, the torque exerted on the bottom plate was measured in order 
to obtain the resultant excitation on this plate directly. For this reason the 
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13. —DEVICE TO MITIGATE VIBRATION: ‘CURVED SIDE SHIELD ON 
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FIG. 12.-EFFECT OF EXTRANEOUS EXCITATION ON VERTICAL DEFLECTION 

if 
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meter. process was recorded on movie with 

the torque measurement (Fig. 16). Some of the essential findings are reported 

+ Depending on the flow regime, rather different conditions for the vibration — 
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2 FIG. 14.-SECOND TEST ARRANGEMENT. TYPES OF AN IDEALIZED na 
TESTED IN COMBINATION WITH HORIZONTAL 


, B,C WHERE HYDRO- 
DYNAMIC FORCES ARE TRANSMITTED TO THE GATE 


(CASES OF LOADING) 


oF. The similarity deviates because a waved water surface exists and the 

_ approaching velocity is } not an independent variable. For the free discharge, 

_ gravity forces, _ which are negligible in the first-mentioned regime, have — a 

Comainant significance. Furthermore, an air space beneath the = begins to 
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“interact with the flow asa water exists behind 
“the gate. In the case of a submerged nappe the vibration is affected in addition 
x by the turbulent motion of the surface roller (Fig. 16(b)). A further complica- 
& arises from a multitude of possible geometric boundary situations. e 


Fhe steady state of gate vibration and , hence, a functional — ode the 


ion is impeded mainly by three effects: 
The air space the nappe usually disturbs the vibration process = 


irregular mot 
along in slugs by the jouer jet. For high depth of overflow, ‘however, the 
tation is large enough to control the surface Suctuations well as the 
- 2. The motion of the surface roller, being in general aietetie tends: 
. ei put the vibration out of paase, particularly at high submergence of the nappe, 
at small depths of overflow, or for nonaerated flow. At extremely low sub- i: 
-mergence, on the other hand, the roller action can be controlled, andevena 
magnification of the excitation a kind of regenerative process 
oes 3. _ Structural parts of the gate cchenting into the the wake were found to be 
specific sources” of vibration disturbance ‘whenever they» were close to the 
= or the lower jet. As the nappe was brought in contact with the gate, 
(for example, with crest type C), it fluctuated with increased frequency but | 
diminished “amplitude only below the line of contact. The ideal shape of a 
a _ with both overflow and underflow | would follow the shape of the wake 


oscillograms from | torque are presented in 
_ Fig. 17. Before evaluation all superimposed secondary effects, mostly ee 
: _ free \ vibrations of the measuring equipment, were eliminated graphically. 7 
The remaining irregularity in frequency and amplitude was taken into account 


When the air space in the wake was held at Slatin pressure by le 


coincided with the lowest frequency, presumably because the virtual mass 


air vent higher frequency and lower excitation always obtained 


"mentioned by the air space beneath the nappe which grows as 
. _h/w becomes smaller. The critical values of h/w below which vibration could 
= be maintained are shown2 in Fig 18. The « effect of structural parts close to 
_ the surface of the lower jet is demonstrated in Fig. 20. The lower the bottom | *) 
plate was placed, the higher were the measured frequencies . From Fig. 21 it 5 
evident tt that the corresponding was decreased. With the 


‘Modified excitation were determined 
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_ VIBRATION 
which a, the ‘amplitude of the mean pressure difference between the 
upper and lower sides the bottom plate, as from torque measure-— 
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in which ap is t amplitde of pressure at a point of 


“Types A,a ; w=25em; iy* Bem 
air vent plugged 
zs © air vent $= 09cm 
air vent 3.6cm 
h/w = 6/10 
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eo 4 air vent plugged 
© air vent 09cm 


FIG. 19.- _MEAN F FREQUENCY f | OBTAINED FROM susie 
ie pan are sonae distributed o over the bottom plate. Actually the pressure 
a amplitudes were larger toward the ” of the plate. The exact value of Ce could 
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_ All m measurements were performed for constant headwater and tailwater rele- 
vations, Comparative experiments with rising or falling tailwater resulted in 
similar vibration phenomena, mostly with magnified excitation, 


_ ATTEMPT TO ) ANALYZE THE MECHANISM OF EXCITATION 

In Fig. 23 ; an attempt has ‘been made to sketch s six consecutive flow p. patterns eB 


taking iid = bs i, of the vibration process for the case without any 


- 


complete submergence (upper jet remains on »aurtace) 3 
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“Wa Ayah FIG. 20. —MEAN FREQUENCY PRESSURE» 
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of the above-listed disturbing effects. 6. (The 8 - bubble shown in Fig. 
sidered to be too small to disturb the excitation.) sor aahasiaeatey 
‘Tt ean be shown that for the free-surface potential flow around ; a een 
- flat plate (Fig. 22) there exists an unstable state of equilibrium analogous to 
_ the one for a cylinder in parallel flow demonstrated by Foeppl. 13 There are 
— two surfaces of discontinuity or mixing zones. The upper one may be approxi- 
= ‘mated by a cycloid which essentially depends on the velocity Vy and the magni- 7 
q of the negative pressure behind the overflow edge. (M. Cc. apes were 
gives an analysis of three characteristic nappe shapes. ‘These shapes were | 
- also recognized in the case of combined overflow and underflow. For certain _ 
flow conditions they all occurred in periodic sequence.) The configuration of i 
oy vortices contained within these surfaces is unstable. In a real fluid the fie 
_ tices tend to grow until finally they are carried along with the flow. An infini- - = 


-_tesimal disturbance is sufficient to initiate a periodic shedding of 


caveat was demonstrated by means of dye injected on ‘the upstream side 
of the gate.) _ Therefore, the clockwise vortex grows and increases the insta- 
‘ - bility of the upper surface of discontinuity until, in — 22 (f), this upper sur- _ 
-- The simultaneous fluctuation of the nappe can be explained by the super- . 
y posed circulation around the plate > which, in Fig. 23(a), reduces the | velocity — 
a Vy and augments” -v. Consequently, the trajectories of the nappe tend to move > 
foward which | is accordance with equations 


tion to the velocity “change, helps to bring the nappe closer to the plate on 
account of pressure reduction beneath the nappe. ._ At the phases r represented _ 
by Figs. 23(d), (e), and (f) the lower jet is decelerated so that more and more _ 

ad retarded and reversed fluid enters the wake andincreases the pressure there. _ 
_ Together with the acceleration in the upper streamlines this causes the —- 

reason why the aeration of the nappe plays ar an important in 
' vibration phenomenon is that the periodic pressure change described above 
does not develop to | the same ) degree when t the nappe is ventilated. Neverth = 
a less , vibrations have been experienced | with the ventilated nappe as well. 

4 Their violence depends to a large extent on the size of the air vent itnetaniia 
; air demand seeG, H. Hickox. ore aif the nappe submerged by the surface roller 


a2 ala force pulsation exerted on the bottom plate Fig. 14) finds its explan 
- | = in the periodic supply and evacuation of separated fluid, taking place par- 
oe aa Mathematische Stroemungslehre, by W. Mueller, Springer, Berlin, 1928, p. 108. 
fee 14 "Sur Ia forme de la nappe des deversoirs sans aerage,” by M.C. _Woronetz, C. 
hebd. Seances Academie Sci. 238, 1954, p. 1688. 
“Aeration of of Spillways,” by G. H. Hickox, Transactions, ASCE, Vol. 1944, 
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23(b) up to Fig. 23(c) where it occupies the entire wake. As a result, the _ 

“ial _ velocity gradient and the vorticity in the lower boundary of the wake willin- __ a a 
“eis _ ¢rease to such an extent that this surface will gradually break up and form a “a 
s clockwise vortex, To keep the total circulation constant, a counterclockwise 
— 
| 
— 

= 
q 

2 
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= flow is superposed on the flow that is represented in Fig. 23, pero a 
_ downward force on the plate in phases a and b and an upward force in phases 
d and e. The time intervals between the phases of Fig. 23 are not necessarily 
ae ‘The smaller the depth of overflow orthe lower the page oh the —, 


‘FIG. 22.—WAKE BEHIND A VERTICAL mar PLATE F FOR 


E CYC LE OF THE VIBRATION 


DISCHARGE WITH THROTTLED AIR VENT 
is 


"water ities can enter the wake near the toe of the nappe (Fig. 16 (b)). Th 
: ra is the ‘reason that the rate of change in force is greater, sometimes even 
in the upward direction For a large > depth of 
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oscillator, he ty dotted 23, superposed on the main 
flow then produces a transverse f force known as the Magnus effect. Due to the — + 

convex | boundaries of the wake, an acceleration of the adjacent flow causes ‘ 

Mi pressure drop and vice versa. Consequently, inthe lowest position the wake is i 

- acted on by a lifting force (Fig. 25(c)) and in the highest position byadown- 

force (Fig. 23(f). These forces may, hence, be defined as restoring or 

i spring forces, Furthermore, there are self-controlled driving forces due to the 

breaking up of the surfaces of discontinuity that take place with a phase dis- 
placement Of 1/2 against the reversal of wake motion, as explained above. 7 

F These two groups of forces combined with the inertia forces of the wake — 


that 
1. The “spring force” varies linearly with the wake ‘deflection’ x. 
The “driving force” is proportional to dx/dt. 


3. wake mass m remains constant during vibratior 


‘The term parenthesis represents a lift coefficient for 
- displacements a<7 /30, 7 is a form coefficient, and as 0 denotes the posi- 
tion of the profile “axis at which cx= = 0. The slight variation in direction of the a 


16 “Formation of Vortex Streets,” _ by G. Birkhoff, Journal of Applied Physics, , Vol. 
_ 17 Mechanical Vibrations, by J. P. Den Hartog, 4th ‘elites. McGraw-Hill Book Co., 


Technische Stroemungslehre, by B x, Berlin, 1958, p. 
8 


With this concept of the periodic flow conditions, and with reference toa 
_ hypothesis about the formation of the von vortex trail!6, an expres- 
sion for the frequency of the excitation may be dervied. The basic idea is to 
consider the wake behind the obstacle as acylindrical oscillator free to rotate 
+ 
| 
— tie Uriving 1urce, cory as tie hegative damping, 
ae i _ is small with respect to the other forces and smaller than 20% of the critical +. % 
damping,17 then the frequency of the wake vibration may be approximated by 
the natural frequency of the undamped system y 
Because the cross s airfoil 
(Fig. 22(b)), an equation m width in 
— 
“9 
approaching velocity v due to the angular displacements is neglected. 
— 
cal 
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‘The mass per unit width can thas 


in which w* is the cross-sectional ‘area of the mass in of 
height w and p is the at of water. By bia Eqs. 9, 10, and 11 in | 


is” obtained. (At the limit dg > do and he y > ~ the square root of the the last — 4 


“equation corresponds to the Strouhal number. 19 adie 


functional relationship the coefficients 7 and ¢ and the remaining 

ables of the problem. The mass coefficient ¢£, for instance, should account for | 
the: added mass in the adjacent jets and the surface roller, and also indicate 
_ the influence of structural parts within the wake that prevent some of the water 


from Nevertheless, ‘it is remarkable that if one ‘replace v 


are obtained for + various tes test t series in the cas case of nonaerated flow flow. lo sep 


1. In contrast to conditions of gates subject to either overflow ‘or underflow, — 
in which the occurrence of vibration is very limited, vibration under conditions ~ 
of simultaneous overflow and underflow persists through almost any variation ¥ 
of flow conditions as long as the depth of overflow is greater than a critical a : 
decreasing tailwater depth, the exciting forces increase up to 
. stage at which a free water surface begins to develop beneath the nappe. The 
_ disturbance of the excitation thus caused will set in earlier and become wai 
_ intensive with increased nappe aeration efficiency. 
3. Means of mitigating the vibration that are not objectionable from a sii 
“9 tural or functional viewpoint are found to have only limited effect. The a 


tion will be efficiently minimized only when favorable types of -_ and apron oh; 
areselected 


a. 4, In general the structural elements of the gate, exposed to the pulsating 
forces, should be withdrawn as far as possible from the confluence of upper :% 
and lower jet. ' This is accomplished most effectively by a sloping downstream | - 

_ floor. The slope should not exceed 1 on 1.8, however, because a steeper slope 

will lead to separation and instability of the lower jet. On the other hand, there 

~ is a critical distance between gate bottom and lower jet. Reducing this —— 
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= 5. From comparative experiments ona fixed anda flexibly suspended model © 
hi gate it has been found that the complex phenomenon of a self-excited and multi- — 
ia ‘s ply coupled gate vibration can be approximated by a forced vibration, the exci- | 
tation of which is to be determined from tests with the gate in fixed positions. 
_ a . Exterme magnification of the gate deflections can be obtained by means 
: ofextraneous control. This is important with gates in navigable rivers, avd, 
_ pressure waves induced by ship engines may adversely affect the gate vibration. — 
— frequency of the excitation is defined satisfactorily by a coefficient — 


= to the number. — 
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HYDRAULICS DIVISION 


-_ _—" ‘es on cylindrical flip buckets may be ‘computed by assuming an ir 
rotational vortex flow. Correlation between maximum theoretical pressures 
_ thus determined and maximum pressures observed in model studies is ca 
sented. This: correlation shows that, despite the simplifications involved, the 
theory may y provide a good estimate of the maximum pressure, both for « chute 


Filip buckets are often used when large amounts of energy must be disposed 
of at the end of spillways. Flip buckets are not dissipating devices, for dis- J 
_ usually requires the depth provided by a hydraulic jump, and it is safe 
to state that no appreciable energy is taken from the flow by | its deflection due 
to the flip bucket. Erosion in the river bed will produce e enough depth for dis _ 
gipation due to turbulence and air A natural stilling pool thus 
_ formed replaces the man-made dissipator, | which would usually be very ex- 
pensive if the flow and the head are considerable, 
i the design of a flip bucket, prime consideration is given to” the force to ~~ 
which the structure is subjected asa result of effect of the carvi- 
Meer 
. _ Note.—Discussion open until February 1, 1962. To extend the closingdate one month, be 
a written request must be filed with the Executive Secretary, ASCE. This paper is part | aa 
, 3 of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- J 
ciety of Civil Engineers, Vol. 87, No. HY 5, September, 1961. 
Senior Design Engineer, Tippetts-Abbett-McCarthy-Stratton, Engineers and Archi- 


« tects, saat York, N. Y.; Lecturer in Civil Engineering, Columbia University, New York, | 
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_ A first approach to the problem isto ) compute the overall dynamic effect, or 


‘. total force caused by the jet deflection, which can be determined by applying a o 


the momentum equation. For this, the assumption is made that the true jet 
_ deflection is parallel to the slope of the bucket lip. ay a eee age 
In addition, ‘it is also “hecessary to to determine the maximum pressure that 


signer, in many ‘cases, should. seek a laboratory endorsement for his omnia 
- Although the laboratory should have a very early role in the design stage, the - 
_ designer first should resort to previous experimental results on similar struc- 
- tures and to theoretical studies. His” conclusions will then produce a layout | 
suitable for experimental work, When the ‘magnitude « of the structure does not. 
_ warrant a model investigation, the designer must analyze the forces involved 7 
without specific experimental confirmation, = 
: = This paper presents the results of studies made to compare the pressure 
. ps on cylindrical flip buckets as computed by theoretical considerations with that 
measured in model investigations, 
ak Notation.—The letter symbols» adopted for | use in this paper are defined 
a where they first appear, in 1 the illustrations or in the text, and are arranged | 
of reference, in the Appendix. 


4 


Previous studies on flip buckets were presented by D. B. Gumensky, 
- assumed that the angular ite il of the water on the flip bucket would remain — 
H. Douma,’ 3 assumed that the velocity distri- 


who 


pattern, all the streamlines are circular and concentric with the bucket. ‘This | 
assumption has also been adopted for the studies described herein, 
to Fig. 1, the upper streamline of such flow has a radius equal 
R - d; the bucket radius is R, and the thickness of the vortex is d. The differ-_ 
ence in piezometric heads between any two points (1 and 2), is located respec- 7 


ms on the upper streamline and on the bucket, is } given by Bernoulli’s equa- 


in and and V9 a are the velocities at points 1 and 2 respectively. Thus: 


— Side Walls in Chutes and D. B. 


asce, Vol. 119, No. 2675, 1954. 
J. H. Douma on 2. 
4 “Hydrodynamics, * by H. ‘Lamb, Cambridge 1932, PP. 19. 
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bia _must coincide with the dynamic force computed by momentum considerations. q 8 
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BUCKETS. 
Since the pressure at point lis s atmospheric: 
dite 
in which aicaie the difference in elevation between points 1 and 2, 
at point 2, and is the “specific re of the water 


“Fic. id |.—CHUTE FLIP BUCKET 
ems may be noted that Vy ‘should vary for different locations of point 1 on the 
_ upper streamline, since H; is variable. _ However, the irrotational vortex pat- 
ae asusmes that Vy varies only with the radius, and must be constant over © 
any streamline. For a high dam, this inconsistency may be neglected, because 
& variation of H; is relatively small. Therefore, point 1 may be selected at 


the lowest water-surface elevation o over tl as bucket, with Hy =Hy. emule i 
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From ‘Eq. 6, , d may be computed by developing the ‘the logarithm in ‘series: ee 
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+ 
R 


— 


| 


tion of the jet thickness. 

ae The theoretical value of the pressure compared by Eq. 4 for any point of 
the bucket surface should not vary much for a high dam spillway. ™ However, 7 
the maximum value corresponds to the lowest point of the bucket, because then 
is maximum and equal to d. 


Thus, the maximum pressure is given by: 


Once ne: maximum pressure has been the theoretical pressure 
diagram may be determined as shown in Fig. 2. A circle is drawn with center 
at the bucket center and radius equal to R+ h. Then, the pressure at any point — 7 
—_— A may be represented bya: radial segment AB equal to h minus the elevation of q i 


_ point A with respect to the lowest point of the bucket surface. = 
7 ” Since the curvature of the flow is assumed to start at thestation of the point 
: f tangency between the chute and the bucket, the theory introduces a discon- ™ _ | 
tinuity in the pressure diagram. Similarly, a discontinuity results with the | 
- computation of the 18 from Eq. 6. In fact, the flow curvature must start up-_ 
_ stream of the bucket.2,5 For a thin jet, however, it may be assumed that the — 
curvature starts to the tangency point. 
Fig. 2 shows a comparison between the shapes. of the theoretical pressure 
and of diagrams determined by the experimental investigations re- 


ae a From Eq. 4, Po can be de nea, p tiat the thickness of the jet is f is ; 
_known. This is accomplished by integr 
— 
— 
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EXPERIMENTAL 

\~ PRESSURE DIAGRAM 


DIAGRAM 


THEORETICAL 


2. —THEORETICAL AND ‘EXPERIMENTAL DIAGRAMS 


FIG. 3. —TUNNEL FLIP. 
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_ The foregoing reasonings have been einen — in mind the case of achute 7 
zz bucket. In te ont of a tunnel bucket (Fig. 3), like those described by a 
. 3 would be less amen in theory. Its u use 


iy computed with the theory and the values determined experimentally. The the- ey 
ory is in good agreement with the values of the pressures | obtained for the 
similar case of pipe curves. 7,8 
For a tunnel flip bucket in Eq. 10 should be taken equal to the velocity 
P _ just upstream of the bucket. Since the change in flow direction is much less b 
for a tunnel flip bucket than for a chute flip b bucket, the flow curvature should A, 

- begin closer to the tangency point in the case of the tunnel. This also — sts 
the use of the depth just of the bucket as a 

; In the two- dissensiounl case ofa straight chute spillway of constant width, — 
the following simple relationship may be derived by dimensional analysis for 


: in which (Fig. 1) H is the reservoir elevation with respect to the lowest point 

of the bucket and H, is the head over the crest of the weir. Eq. 11 does not 
include the terminal angle of the bucket, this fact being consistent with Eq. 10. 
> The United States Army Corps of Engineers has performed model studies _ 
at the Experiment | , Vicksburg, Miss., on the chute spillways 

These constant - width spillways have 

eytindrical buckets of the type in this paper. For Hartwell Dam 
spillway bucket, radii of 30 ft and 40 ft, and lip angles of 20° , 30°, and 40° 
were: studied. ‘The | chute slope was 1.00 vertical to 0.67 horizontal. | The max- rf 
imum prototype head over the bucket was H = 160.4 ft and the maximum dis- _ 

_ charge approximately 1,000 cfs per ft. For Pine Flat Dam, the bucket i 
was 50 ft and the lip angle 20°. The chute | slope was 1.00 vertical to 0.78 


horizontal. _ The maximum prototype head over the bucket was He _ 3 ft and 


Se maximum ap proximatel 1, ,350 cis | or ft. ft. 


Strausser, Proceedings, ASCE, Vol. 81, No. 747, i955. 
8 Fluid Mechanics Tor Hydraulic Engineers, by H. Rouse, McGraw-Hill Book | Co., " 

"Army Engineer Waterways Experiment Station: Technical Memorandum No. 
2-375, Vicksburg, Miss.,1953. 
10. S. Army Engineer Waterways E: Experiment Station: Technical Memorandum No. 


2 393, Vicksburg, Miss., 1954. 
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‘FLIP BUCKETS 


oa with the model studies, ‘but, s so far, for discharges much smaller than 


_ The Hartwell Dam model studies showed that, for the same bucket radius 
: 


and discharge, the maximum pressures werev verys similar for bucket lip angles — 
: of 30° and 40°, For 20° angles, the maximum pressure was no less than 93% of 
‘ that corresponding to 40°. This departure tended to decrease for the press 
discharges used in the studies. Thus, the ‘assumption that the maximum pres- 
sure is independent of the lip angle appears correct, at least for the angles — 
most commonly used in American practice. In any case, for small lip se 
the maximum pressure determined t by comparison with th model studies 
Ve would on the safe side. plot of 11 with the data © 


ing to this particular layout with those « of two projects that have been tested 4 
‘Fig. Sisa of observed maximum pressures for Hartwell and Pine Flat 


With a standard error of 3.4 ft and a correlation coefficient of 0. 999, This 
correlation | is good, despite the wide range | of variation of h and of the lip — 
- gle a. The coefficient of Eq. 12 takes care of the simplifications made in the 


‘a cal pressures ‘according to Eq. 10, but with d ae using Eq. ud . The best 
fit straight line be expressedas: 


«thes 


with a standard error of 6.4 ft anda correlations coefficient of 0.997 i 
. Eq 13 shows that the irrotational theory provides a good estimate of the 
experimental values of the maximum pressure, a The simplifications made to 
compute theoretical pressures for Eq. 12 result in a systematic larger differ - 
- ence between these values and the experimental ones. When the depth is com- 
_ puted according to the irrotational theory without simplifications, the system- 
atic departure between theoretical and experimental values is only 2%. How- — 
_ ever, the use of the corrective coefficient results in good correlation for r both 
Eq. 12 and ‘Eq. 13, and the former is easier to use. 


7 to depend on the lip angle a, being close to the lowest point for lip angles of 


for the computation of the theoretical maximum pressures, and the rest of the | 
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not enough to draw many contours, it shows a very definite trend. _ This plot 
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1961 
pressures on the bucket are computed according to the theory applying 
roy the same coefficients, the resulting diagram will be on the safe 
. This may be done by increasing the pressures on _ 
the theoretical diagram by the between and , denoted 


MODEL STUDIES ON TUNNEL SPILLWAY F FLIP BUCKETS 
- The Bureau of Reclamation, United States Department of Interior (USBR) 
od has conducted model studies on the tunnel flip buckets for Flaming Gorge, _ 
Glen Canyon, and Whiskeytown Dams, In a six different buckets with var- ae 
_ iable discharges and conditions were studied. 
ia Fig. 7 shows a plot of observed maximum bucket pressures for these pro- - 
jects versus maximum pressures computed to The best fit 


line be expressed as: 


— 


with: a standard error is 8. 3 ft and a correlation coefficient of 0. 995. 
It may be observed that the correlation is good, although the petencten ts is 
2 | a somewhat larger than for the chute spillways. _ This should be expected since 


a _ the geometry of the bucket may vary from case to case aaa 4 than for  - 

ad In these studies, the maximum pressure was observed at a distance of 0. i. 

q R to 0.18 R downstream of the point of tangency for an angle a of 15°. This 
a distance increased to 0.35 R, for a = 35° -= 0.42 R to 0.45 —— a= 38°. * 


_ MAXIMUM PRESSURES ON FLIP BUCKETS y 
CONVERGING CHUTE SPILLWAYS 


charges over the bucket may not be uniform, and maximum dischagre may oc- 
= cur in the vicinity of the centerline. As a consequence, the pressures over the 
center of the bucket are greater than the value computed by Eq. 10 which is 
a ‘based on the assumption of a uniform discharge over the bucket. Conversely, 
in the vicinity of the side walls, the pressure maximum is smaller than the 
a computed with Eq. 10. _ The variation of discharge per | unit width changes © 
with the. convergence, the length, and the slope of the chute, (aa 
; _ Recent experiments were made on two models of the convergent Shihmen 
Spillway (Taiwan) by the Shihmen Development Commission 1 by arrange-— 
ment with the Hydraulic Laboratory of Taiwan University, Taipei. In these 
=. total convergences between the side walls of 40° and 20° were used, as is 
‘measured by the ange of the circular weir. . The final chute ‘slopes 
1.72 and 1.47, naively. ‘The walls had transition curves in plan, so that 
they became parallel to the centerline upstream of the bucket. 
From the results of these tests, the pressure computed by Eq. 10 for a aes 
_— form discharge over the bucket approximates the average maximum pressure © 
observed in the model for the same total discharge but with the actual dis- __ 
— distribution. The pressure distribution along the bucket was more uni- 
form for the model with 20° convergence and 1: 32 ‘final chute _— than for the. 


| 

— 

| 4 


me with 40° convergence and 1:1 final It was also observed that 
eae rise started at about two depths upstream of the bucket tangency | 
“4 point. This fact is of importance to the location of a joint between the bucket “4 | 
_ and the chute, as the joint should be located upstream of the region of et 
: rise to prevent the transmission of high pressures through the joint and crea- 
_tion of excessive uplift pressures. = 


_ More investigations are necessary to draw definite ® conclusions about the the 


Maximum pressures on spillway flip buckets may be easily computed by as- 

7 _ Suming an irrotational vortex flow. The pressures measured on models of both — 
straight chute and tunnel flip buckets show good correlation with and are close 7. 
_to the theoretical values of those pressures. The theoretical formulas may be 

corrected, by the use of suitable coefficients, to conform with the experimental a + 
results. ‘The distribution of pressures over the bucket may be approximated _ 
oy a simple graphical procedure . Although the results of model tests on con- 
_ verging spillways show qualitative agreement with the theoretical pressure . 
computations, more studies are required to establish the relationship between — 


the —*, of the spillway and the pressure e distribution on the bucket. od 
i 


Tippetts- Abbett- “McCarthy- Stratton, Engineers and Architects, New York, 
aa N. Y. Some of the computations and figures were made by Herbert Meltzer, 
7 _ Chung-Li Ling and William B. Tippetts, engineers of the Civil-Hydraulic Sec- 

of the said Department, Arthur R, Luecker and Gordon R. Williams Tame 

_ Associates, made very useful suggestions andcomments, 2 


The is indebted to the United States Army Engineer Waterways 
periment ‘Station and in particular to its Director, Edmund H. Lang, for per- 
mission to publish data pertaining to Pine Flat and Hartwell models. 

_ The writer also ‘gratefully acknowledges the information on tunnel spillway way 
APPENDIX.—NOTATION 


The following for use in the conform 
wt “American Standard Letter Symbols for Hydraulics” (ASA Z10.2- 1942) 


a constant equal to of the voi vortex the cf 


4 thickness of the vortex; 
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= acceleration | one to gravity; 


= head over the crest st of ants ~ illway; : 
° 


’ oon = total head over lowest point of water surface, ignoring head loss an 


— 


= total head over any point 1 on the surface of the —— ignoring =| 


woe = maximum head, g given in euhenn of water, on on the bu cket; es 
exp 1 = experimental value of the maximum pressure in column of water; 7 


htheor = theoretical value of the maximum pressure in column of water; a ‘ 
= pressure at any point of the bucket; 


monty = dischar e per unit width of 


= radius o the bucket; a 


velocity at any strea 


vy velocity at the vortex surface; 


= Peservolr elevation Over the lowest Doin OF th bu Ket: 
| 
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ESTIMATING THE PROBABLE MAXIMUM PRECIPITATION 


ia 


of estimating the probable maximum precipitation is proposed. ‘Empirical 
distributions of the spatial and temporal sampling errors enable the analyst 


used to obtain the extreme precipitation. A comparison of the results from 
the traditional method with those based on the method presented herein. shows — 
differences in both directions which occasionally exceed a factor of two. The 
emphasis has been on an examination of and | rather than onasearch for _ 


to estimate the magnitude of the errors” attached to the statistics which = 


INTRODUCTION 
? ‘The procedure originally used for estimating the probable maximum pre- 
cipitation (PMP) during the 1940’s and still in use in the 1960’s, with slight 
. modifications, uses a combination of a physical model and several estimated 


meteorological parameters, Unlike most investigations in which the a 


— can be compared with actual observations, the quality of the estimates obtained 
from the traditional approach cannot be evaluated. The fact that a PMP estimate _ 
has not been exceeded is no standard for evaluating the quality of the estimate. 
The lack of a dependent va and fact that rainfall is the most intensive- 


_ Note.—Discussion open until February 1, a To extend the closingdate one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
_ of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 5, September, 1961. 
_ 1 Meteorologist, Hydr. Services Div., U. S. Weather Bur., Washington, D. 
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oly of all the 1 elements ts suggest an approach 
ae _ Attempts have been made to by-pass the traditional approach by graduating 
i a sample of rainfall data according to some extreme-value distribution, but 


eoverns serious weaknesses were € encountered. First, there was no reason n for 


ts curve presented a theoretical rather than a real picture of what had ones 
_ or would occur. _ Second, there was no basis for stopping at a particular prob- — 
ae ability level. Third, little information was available for a judgment about the | 
Dat magnitude of the sampling errors attached to the estimated distribution n para - 
1961, there. are which should be considered when estimating» 
PMP; more data and additional statistical techniques. Data for thousands of © 
in stations are now available -ina a form suitable to probability analysis, that is, 
A the annual maximum 24-hr amounts have been extracted from climatic tabu- 
at lations. A much smaller but significant quantity of 1-hr data is also available. 
_ Statistical innovations include several new extreme-value distributions and 
of their relationship to each other. 4 
_ _The primary purpose of this paper is to propose certain sysematic methods _ 
_ for the analysis of these rainfall series with the goal of reducing the range of : 
uncertainty surrounding the estimates. The analyses presented herein 
* are based on 24-hr rainfalls from about 2,600 stations with a total of 95,000 _ = 
station-years of data. A combination of several statistics from each station’ s 4 
data was computed and viewed in the aggregate in order to arrive at an en- 
_ voloping statistic. Because reliability is a major design parameter, empirical _ 
7 z measures of several components of sampling error have been determined for 
- each statistic. Empirical adjustment factors have been derived for such data 
Z deficiencies as small sample size and calendar- -day rainfall rather than the _ 
1,440-min maximum rainfall, Comparisons between the traditional approach 


7 and that developed herein are _—— for two large regions of the all 


"ANALYSES 


PMP.- paper, the PMP is defined as the largest rain- 
fall (precipitation) that a station is ever r likely to experience f for a particular 7 
_ Basic Data.—The data from a key gr group of 198 long-record 24-hr stations 74 
(164 first-order Weather Bureau United States Dept. of Commerce (USWB) and 
34 cooperative observer) were used to develop most of the relationships pre- — Be 
ne Daily rainfall records from an additional 2,400 stations were analyzed in 
Rainfall Sample. —The extreme rainfall experience at a station can be ll 
> resented by a sample composed of the annual maximum rainfalls for a particu- — 
lar duration. Each extreme, or annual maximum rainfall, then represents 
a relatively large | number of actual observations, because each value in the . 
series is itself the extreme of a large number of observations. This series is 
amenable to probability analysis because the items are causally 
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j 
year is Bie is as assumed that the remains the 
same during the sampling period. 
Standardized Variate. : 
= Definition.—Yen T. Chow,2 F. ASCE demonstrated 1 that the only difference 
fa? ‘between | the various distributions which lend themselves to the analysis of rs 


extreme- -value hydrologic data is common statistical variable, K 
Le 


in which xp is the rainfall he's re urn T, in years, when a 


. _— -value distribution is used. The terms ¥y and d SN are the mean and 
standard deviation for a series of N annual ‘maxima. 
if the maximum observed 1 rainfall, xj, is substituted for xT, and Km, a 

‘kK, then Ky is the number of standard deviations that must be added to the - 

; mean to obtain xy. The variable Ky isa common denominator which admits ot 

_ comparison on a probability basis in two ways: (1) it will tell us how rare a 

= observed rainfall at one station is with respect to the maxima a 
other stations; and (2) how rare the rainfall is at a station if it is “et : 


‘that the data follow a particular distribution. The former will be presented _ 
Definition of Ind Independent Ky. —The variable, was determined for e 
of the 2,600 24-hr stations by computing the mean and standard deviation by the i 
conventional procedures but omitting the maximum observed rainfall from the 
computation. This _ computation is equivalent to observing the maximum — 4 
observed rainfall after Xn and snhave been computed, The distributions of the 
Ky’s are given in Table 1. Neither the maximum observed rainfall, the mean ¥ 
“nor the standard deviation was adjusted for such data deficiencies as short | ; 
and observational-day rainfall, The term, Ky, as used henceforth 
this paper refers to the independent computations. ne 
7 a _ Independence and Randomness of Kjy.—In order that any inferences drawn 
from the distributions of the independent Ky’s be valid, the / assumptions of :. 
independence and randomness must be satisfied. Plots of some of the larger 
_ Ky’s show no systematic relationship between rainfall magnitude and geo- ; 
graphy. ‘Examination of the Kyy’s for 68 long-record stations in Iowa, a region 
where the rainfall regime is not influenced locally by orography or bodies of 
. water, shows a range from 2.2 to 11.2. The unusually large storm of a . 
1943 on the west coast of the United States produced the maximum observed | 
cr hr rainfalls for 33 of 50 stations in Los Angeles County, Calif. (5,000 sq 
miles), but the 's associated with» these maxima vary from 2,2 to 7.1. 
“These are indications that the Ky’ s occur in a random fashion for raingage 
networks as dense as the two investigated and previously noted, For practical — 
purposes, the | entire sample of 2,600 Kjy’s may be considered a random aie sf 
7 6h a common population of rainfall data regardless of the proximity of the oe 
stations or the storm which produced the maxima, — 
Relationship Between Ky and Return Period.- —The second characteristic 
of Ky is that it is uniquely related to the probability of occurrence of an po 
‘ at a station if a theoretical distribution is assumed. For example, the 100-yr 


is associated with a of 3.5 for the ‘Fisher-Tippett type I distribution 


i | 2 “A General Formula for Hydrologic Frequency Analysis,” * by V. T. Chow, Trans- a 


actions, Amer. Geophysical Vol. 32, 1951, p. 
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‘whereas’ smaller K K, perhaps 3.1 , might be 


associated with the two-parameter log-normal. However, it should be em- 


i 7 bias in favor of the 50-yr means, This bias is so because the distribution of g 


a phasized that the quality of the data ord‘narily encountered in practice from 
a re _ Samples of the size considered here does not justify using a theoretical dis- 
— to graduate the data to some rare probability of the order of magnitude oa 
a 107 yr. Use of the meanand standard deviation as estimators for PMP pro- 
vide the amount of accuracy required to exploit the data. eB pene 
General.—Because relatively small samples will be used inthese analyses, 
some examination of the adequacy of these samples is pertinent. A series of _ 
annual maxima is assumed to be a probability aamngee. However, it is not neces- 
sarily a representative sample or, in other words, a miniature replica of the i. a 
_ population from which it was drawn. ‘This is especially true when dealing with | « 
a vast and inaccessible though finite population, Itis also quite possible that the | 
results from the nonrepresentative sample provide a good design estimate . 
because during the entire lifetime o: of the hydraulic structure a representative — 
‘sample is not available. The size of sample is not an indication of its repre- 7 
_ sentativeness or its randomness, but size is related to the amount of chance <a 
variation to be expected from one sample to another. 
Consideration of the principal components of sampling variation is an essen-_ 
thal preliminary to the analysis of computing the statistics which are used in 
- combination to arrive at the PMP. Measures of the uncertainity involved in 
estimating both the mean and the deviation are developed in the 


The Mean of the Annual Series .— 
Sources of Bias in the Mean. —A comparison the arithmetic means 
198 50- records with those for 10-yr 1 records shows a slight positive 
extremes is _ Skewed to the right which means that there is a greater chance — 
of getting a large than a small extreme as the length of record increases, 7 ; 
_ Another factor which can contribute to this bias is the introduction of an out- cd : 
iy lier, or unusually large event, intothe series of extremes. Recognition of = 
factors is an important preparatory step to amend of iis the mean and 
the standard deviation. 
Adjustment of Mean for Outlier.—In a study4 oo in 1960, the 
“Fisher Tippett type I distribution was found to give an adequate description of 
the behavior of rainfall extremes. Making use of this knowledge, pene enol 
-_extreme- value distributions for several lengths of record were developed by 
reading suitable values on extreme-value probability paper (Fisher-Tippett 
type at the intersection of a straight line and the probability 
levels, m/(N + 1), at which m is the rank of each item and N is the length of 
record, The slope of the line, or direction of these hypothetical distributions, — 7 
was made equal to the average 198 station 100-yr to 2-yr ratio of 2. 5. a 
- ratios of two means were computed for several record lengths: the numerator 


mean does not include the maximum observed event whereas the denominator im 
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“os a mum event from three times the second highest event to equality with the | _ 
= “Statistics of Extremes,” by E. J. Gumbel, Columbia Univ. Press, New York, 1958, 
and M. A. Kohler — 


second highest. The relationship between the increase inthe mean as a Sunetion 


of the ratio of the two means and record length is illustrated in = 1. For ex- 


1 effect of the maximum observed event on the mean, No consideration was s given 


-_- records from the 198 key stations were first adjusted for an outlier according — 


to the —— of Fig. 1. The average ratios of 50-yr mean to 10-yr, 15-yr, : is 


. 


ADJUSTMENT FACTOR 


2 


20-yr, and 30-yr | means were then ‘determined. A comparison of the small 
number of available 60-yr and 70-yr means with the 50-yr means showed _— 
ver a negligible difference. Thus, the statistics from the 50-yr records 
_ were - used as a standard to adjust those from shorter records, Fig. 2 shows 


a plot of the 50-yr means versus those for the 10-yr means after an adjust- cm 


a had been made for the maximum observed event. The average bias is 
“seen to be about 5% an average sampling error of approximately 5%. 


oe Adjustment of Mean for Sample Size.—The means from both 10-yr and ae | 7 
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upward by 2% to make it comparable to the 50-yr standard. 
Measure of Spatial Sampling Error. —Daily data were examined from 68 
+ i long- record stations in Iowa in order to. gain some perspective of the spatial 
a variation of the mean. Eventhough a slight latitudinal trend exists, the average _ 
around the overall mean of 2.5 in, is 0.2 in. This that 
| mean of the annual series is is a relatively stable statistic. 
“ Definition. — Another common and easily understood statistic is the 
= deviation. n, Essentially, this statistic is a measure of « dispersion of the data 
around their mean. When dealing with extreme rainfalls, it is an unbalanced _ 
measure around the mean of a distribution which is skewed to the right,—the 7, 
Ts average coefficient of skew for a selected sample of 118 24-hr stations for J 
a common 43-yr record being approximately 1.2, 
i The three principal sources of unreliability in estimating the standard ~ 
a: deviation are the ‘sampling errors due to spatial variation, sample size, , and os 
the inclusion of an outlier into the distribution. The amount of confidence o> 
a be placed on a single station’ s standard deviation can be estimated from the 
_-—,- Measure of Spatial Sampling Variation.—The chance occurrence of several 
_—idarge rainfalls at one station and not.at another nearby station can result in a> 
_ large discrepancy between the two stations’ standard deviations. The problem =a 5 
_ herein is to obtain some measure of this station-to-station variation. The 68 © Hl, 
: -24-hr stations in Iowa were again selected for analysis. As indicated earlier. 
: ‘the rainfall regime in Iowa is not influenced locally by bodies of water or 
“orography. There is, however, a slight systematic latitudinal variation ac: across ’ 
the state which was not considered. A measure of the intrinsic spatial sampling © a 
 varignies is provided by the dispersion around the mean standard deviation of 
0.97 in.. The standard deviations range from 0.70 in to 1.30 in. with two-thirds _ 
of them within the interval 0.81to1.13in. 
_ The Effect of An Outlier.—The effect of an outlier is much greater on -on the 
"standard deviation than the effect previously shown for the mean, For in- 
_ stance ¥ in 1955, (12.12 in, . occurred at Hartford, Conn., during a continuous _ 


and increased the standard deviation by 67%, but the mean was only increased : 7 
by 6%. For shorter records, an anomalous event might increase the standard - 


& -hr period. This amount exceeded the previously observed maximum by 80% 


= 


Adjustment of Standard Deviation for Outlier. test statistic was created 


_ Which based on the same hypothetical series and procedures developed 
P previously for the mean, The statistic is the ratio of two standard deviations — -_ 
and is sensitive to the kind of anomalies expected. The numerator does not 1. 
ss include the maximum observed value in the computation of the standard» = . 
_ deviation whereas the denominator includes all the observations. The graphical _ 


relationship is shown in Fig. 4, An average ratio of 0.85 for a 15-yr record | st ' 


_ indicates et the extreme observation occurred according to a probability P 
distribution defines a straight line on extreme- value probability paper 
with slope equivalent to a 100-yr to 2-yr ratio of 2.5. means no adjust- 
_ ment for the standard deviation. A ratio of 0.70 indicates that an outlier is in 
the series of extremes, that is , an extremely rare event, perhaps, a 1,000-yr 

event occurring during a short period of record, The | standard would 


"then be adjusted downward or reduced by approximately = 
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Distribution of Standard Deviation Ratios. —The distributions ns of 
vs Bod ratios for the key group of 198 24-hr stations are illustrated in Fig. 5 
a in the form of histograms. The ‘dispersion of the ratios around their mean 
_ ‘varies inversely with record length, that is, the random variation about the — 
- = value is considerably smaller when the frequencies are computed from 
long-record groups than when computed from short-record groups. The 15-yr 
4 record histogram shows that approximately two-thirds of the ratios are to the 
deft of 0. ratio or that the for two-thirds stations 
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‘FIG. 4. DIAGRAM FOR ADJUSTING THE STANDARD» 
DEVIATION FOR THE MAXIMUM OBSERVED 


before their average | ‘return period, This is consistent with the probability 
theory which holds that there is approximately a 2:1 chance that an event will 
occur before its average returnperiod. event wl 
_ Adjustment of Standard Deviation for Sample Size.—The effect of record 7 S| 
=a is much more pronounced on the magnitude of the standard deviation 


than the mean, The average ratio of the 50-yr standard deviation to 
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standard both have been adjusted for an is 1.29, 
. compares to the 1.05 ratio for the mean, The comparison of the two standard 
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“FIG. 6.—COMPARISON OF FROM 50-YR AND 10- 


. average between the 50-yr standard deviation, which was 
used as a standard, and the shorter record standard deviations is shown in 
_— diagram B of Fig. 3. Tests comparing the standard deviations from the small 
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ence standard for the estimation | problem of this paper. 
Definition. —The numerical value of the coefficient of variation (Cy) 
Obtained by dividing the standard deviation by the mean, This statistic ex- 
presses the variability of a distribution as a percentage of the mean. In other e - 

words, the larger the Cy, the larger the dispersion of the ‘data around their 
mean, A graphical illustration on log-normal paper - would show the 3 
associated with a large Cy having a steeper slope t than the data 

Geographical Variation. Because the | Cy does show some systematic geo 

_ graphical variation, it is an especially valuable guide for estimating the stand- _ 

ard The mean, or denominator, showed both little temporal vari- __ 

ability for the 198 stations or spatial variability for the Iowa stations. The 

a average Cy for the 2,600 stations investigated for this study is approximately | 

‘i % 36% with 90% of the stations falling in the 30% to 40% interval. The C,’s . 4 
northwest Washington state and on the windward side of the Sierra mountains 

ae are generally less than 30% whereas Hawaii, some parts of midwestern United a 

a States, and the desert areas of California and Arizona, are generally about 50%. _ “is 
pigs This geographical effect is pronounced enough inthese regions to transcend a 4 
“* some of the sampl ling fluctuations described in connection with the standard er j 
A deviation. Also, this effect canbe explained partially by considering the meteo- "s 4 

eS rology of the particular region, For example, the most intense storms in the i q 

_ desert region of Arizona are associated with thunderstorms, Because thunder- _ 

_ storm rainfall occurs in a spotty pattern, the chance that a heavy rainfall will a 

ak occur at aparticular station is quite small. However, the occurrence of several 

extremely large rainfalls and a large number of relatively small rainfalls 

A * over a period of years results inalarge Cy. On the other hand, the stations in # 

— northwestern Washington have a relatively small Cy. Here, the rainfalls which © 
- make up the annual series are associated with storm systems which cover 7 
_ thousands of square miles and occur many times a year. Combined with this bs 
@ - is the orographic effect which helps to extract a large amount of moisture from 

_ each passing storm. The result is a series of annual maxima which vary little _ 

7 any _ Adjustment for Sample Size.—It has already been stated that the mean and > cS 

v standard deviation increase with sample size. Because the standard deviation — a 

_ increases at a much faster rate than the mean, the Cy also increases with sam i 


ple size. This is illustrated in diagram C of Fig. 3 which was prepared from es 
- the: data used to construct the mean and the standard deviation diagram of the — 


Observational-Day Series Versus 1,440 Minute —Most rainfall obser- 
vations are made in nonrecording gages with the subsequent deficiency that the e- 
amounts are for an | arbitrary beginning and ending rather r than for the period ae 
re. containing the maximum rainfall. Adjusting a2-yr rainfall, or any other return- 
period amount, derived from a series of calendar-day or observational-day 
rainfalls by the factor 1.13 makes the new result equivalent to the 1,440-min 2 
for the corresponding return period. The small uncertainty. involved 
in using this relationship is exhibitedin Fig. 7 for a group of twenty — 
SG dispersed stations with a wide range in rainfall magnitude. _ 
storm,—Quite of ‘Day Amount Versus 1,440-Minute Amount for an Individua 
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of two observational- day a amounts. It is necessary to how much 
of the two-day total actually occurred in a continuous 1,440-min period. The | 
= relationship « established for this purpose is composed of the maxi- — 
mum observational-day and one-half of the maximum pment day rainfall as 


“124 storms from widely scattered stations. 
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OBSERVATIONAL - DAY RAINFALL 
FIG. 7. _—RE LATIONSHIP ‘BETWEEN 2. i. 440-MIN 


General. ~All presented thus far have — on n point: 
rainfalls but depths over areas are requiredfor the design of hydraulic struc-— 
tures, A large number of areal depths based on a handful of stations have been 


mates is ‘not sanamensttle A ‘small number of dense raingage networks provide ig 
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“Storm Rainfall 1 th Corps. of Engrs., Army, February, 


rainfall = max ob-day + 1/2 max adj. ob-day @ 
— 
= 
| 


might reliable areal estimates but, unfortunately, 
networks have not experienced an unusually large storm. 
_ Variation with Gage Density.—The data from a 40-gage network, covering / 2 
approximately ~400-sq miles in West Virginia, were analyzed to gain — 
_ measure of just one component of sanipling error, average depths from selected | 
geographically balanced subnetworks were compared with depths” from the 
40-gage network. Even with some built-in correlation asa result of comparing 


with the whole 9 show a deal of of scatter. In Neu 
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MAXIMUM OBSERVATIONAL- ‘DAY * ADJACENT MAXIMUM RAINFALL (INCHES) 


3. 8.~ _RE LATIONSHIP BETWEEN i, 440- MIN RAINFALL AND ND MAXIMUM OB- oe 


SERVATIONAL-DAY + 1/2 ADJACENT MAXIMUM OBSERVATIONA L- 


of anything better, the numerical results the 24-hr area 
curves ¢ of Technical Paper No. _ 38 and presented in Table 2 can be used for 

Unofficial Observations... immediately a after a flood 
uncover some unusually large rainfalls in an unofficial gage or a container | 
Pte 6 “Generalized Estimates of Probable Maximum Precipitation for the United States 


“¢ _ West of the 105th Meridian for Areas to 400 Square Miles and Durations to 24 ecotupens = 
‘Tech Paper No. 38, U. 8. Weather Bur., Washington, I D. Cc. ++ 1960, p. 
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"rainfall maxima, Eighteen of the largest unofficial observations are = oer 
in Fig. 10 and will be presented in a subsequent section where a comparison # 
| made between the results of the traditional approach and that of this paper. 

Empirical Estimateof PMP, 
_Rationale.—Enveloping Kyy from the large number of computed Kjy’s would 

seem t to provide ¢ an acceptable index for estimating PMP. From Table 1, it is 


sidered the appropriate order of magnitude for Ky, because it envelops all 
maximum observed 24-hr re rainfalls at official stations which have a series of 
Method for Computing 24-hr_ PMP.- —There are ‘six fundamental steps in the 
‘g computation of the PMP by the method developed herein. First, the mean ll 


; the standard deviation a are computed by the conventional procedures. ——, i 


‘rainfall according to the relationships of Figs. 1 and 4. Thirdly, the mean and 
_ the standard deviation are onan upward, when necessary, to meet the 50-yr ; 


‘aie step the mean to ‘the product of fifteen times the 


0.93 
— 


‘Finally, and extremely the Cy’s at nearby stations should be 
pared with each other. When large discrepancies appear, the data for several - 
as ‘stations should be plotted on some form of extreme-value probability paper _ ra, 
to determine if the Cy is the result of a trend which has been operative over eS) 
i long period of time or the result of one or more anomalous events. If deemed a 

- necessary, the Cy should be adjusted by changing the standard deviation only. 

: 7 The mean, as shown previously, is a relatively stable statistic. ive fk 
‘a Comparison of (% + with Traditional Approach. — 

Comparison With Hydrometeorological Report? No. 33.—Hydrometecrolog- 
¥ ‘ical Report No. 33, prepared by the traditional approach, presents generalized 

5 7% charts and diagrams for obtaining PMP estimates for several combinations of 

areas and durations east of 105° W. For or comparison purposes, the 10-sq mile 
24- hr depths of the latter report are considered equivalent to the 24- hr 
station values of the procedure outlined herein. Actually, the 10-sq-mile values 

ie: obtained by drawing a circle around the maximum observed 24-hr a 


“Seasonal Variation of Probable Maximum Precipitation East of the 105th Meridian © 


for Areas from 10 to 1000 Square Miles and Durations of 6, 12, 24, and 48 Hours,” Hy- 
drometeorological No. 33, U. S. Weather Bur. D. 1956, Pp. 


le as a bucket, bottle, and so on, The quality of these observations is, of Pe, om) Ae 
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values or by us ‘using some e other simple smoothing pr process. For onputente, 
100 first-order and 100 long-record cooperative stations were used to define ; 


along with some of the large unofficial observations 
Examination of the two sets of isolines shows close agreement in the Gulf of 
Mexico region but then they begin: to ‘diverge: rapidly, the statistical map 

showing a steeper gradient and smaller values from south to north, They 

os in magnitude by a factor of about two in the north and west. Except for 


the three unofficial observations at Cherry Creek, Colo., Boyden, Towa, and dq 
- Smethport, ~Pa., the mean plus fifteen times the standard deviation envelops 
¥ the data. Approximately twenty-five standard deviations plus the mean 
would be necessary envelop» the three large unofficial observations. 
re Comparison in California with Weather Bureau Technical Paper No. 38. —The 
foregoing comparison covers only the relatively flat portion of the United _ 
_ States. It appeared worthwhile to make another comparison in a region where f 
' orography has great influence on the rainfall regime. California was on 
mo the orographic factor and the large amount of data available 
Fig. 11 shows two maps, one basedonthe traditional approach from Technical — 
B Paper No. 38 and the other basedonan analysis of more than 450 stations. The 
= were positioned by interpolating between the stations (solid dots) with 
c= attention paid to topography. Examination of the (x+ 15sy map reveals a 
range in rainfall magnitude which is approximately 50% larger than the Tech- , 
_ nical Paper No. 38 results. The range forthe former being from 8 in. , to 60 in, 
the latter | ranges from 10 in. to45 in. with the discrepancies between 


— study has been based on that valuable is 
; buried in each station’s sample | of extremes, and that this sample can an provide 
_ estimators which can be used in conjunction with enveloping statistics to esti- 


‘mate PMP . The results based on this method have been compared with those > 


obtained from the traditional approach and reveal that the estimates from the —s a 


_ latter approach, which rests heavily on meteorological parametersof unknown © 


quality, tend to be much larger than the statistical estimates for the — on 


o States east of 105°W - In Califronia, the differences are in both directions with iG 
the Technical Paper No. 38 estimates, in n general, being slightly smaller, 


‘The seriousness of | the sampling variations disclosed for the statistics used 


to estimate PMP can be appraised only with reference to other procedures ca th: 


* Stenting 0 . It is hoped that the material presented herein will provide a 


&g 
— 

— 
ye 
| 
| 
4 
— 
shal 


September, 1961 


Journal of th 


Proceedings of of the American Civil | Engineers 


"TERMINAL SHAPE APE OF A A SHALLOW LIQUID FRONT 
E, Roy Tinney, A. M. ASCE and D. L. Bassett? 
4 


7 
sides liquid front attains s a terminal shape and velocity a short time © 
ee! the liquid is introduced into anopenchannel. The equations for the shape oy 
of such a front are derived for laminar flow from Navier-Stokes equations. 7 
"simplified differential equation for a turbulent front is then presented with 
: — solutions obtained by a high speed computer for various values of 
‘the relative roughness parameter. It then shown these profiles 
the same as for a laminar front when a modification is made in the length | 
parameter to account for channel roughness, 
- Concomitant experimental investigations with oil verify the equation for the 


shape of a laminar front l-in, deep, and similar experiments 


fold variation in ‘slope, and roughness, Both the and 


the tests are for 2 an impervious bed with zero infiltration . ae 7 
enctniines of a Front.—When a liquid is introduced slowly onto the upper 


flume), the wetting front appears to. stretch out into a gradually tapered 


- Note. —Discussion « open until February 1, 1962. To extend the closing date one month iss 
" . written request must be filed with the Executive Secretary, . ASCE. This paper is part — 


= 


_ ciety of Civil Engineers, Vol. . 87, No. HY 5, September, 1961. 
_- 1 Head, The R. L. Albrook Hydraulic Lab., Div. of Industrial Research, and Pro 
= of Civ. Engrg., Washington State Univ., Pullman, 
8 Asst. Prof., Dept. of Agric. Engrg., and Asst. Agric. — Washington Agric. ix- 
Stas., Washington State Univ., Pullman, Wash. 
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large, ‘but it soon decays as the slope of the water surface decreases. Finally, 
the velocity of the front becomes constant. From this moment on, the process 
is simply one of extending the length of the uniform section behind a front | 
whose shape remains the same. (Figs. 1 and 2). 
Bi Notation.—The letter symbols adopted for use ‘in this p paper are defined and > 
arranged alphabetically, of reference, in the Appendix. 


“@atance the: upstream depth becomes equal to the normal 
i depth corresponding to the discharge, slope, relative roughness, and Reynolds 
humber and, as stated previously, the front adopts a terminal shape and con-_ 
stant. velocity of advance. If, in an element of time, ¢ dt, the front moves a dis- 
tance, ds, it follows from continuity that the normal be exe 
_ tended by an equal length. Thus for the 


for a unit width of upstream eection 


‘Thus, the velocity of edvence average velocity in the up- 
stream uniform section. This holds for both laminar and turbulent 
Shape of the Laminar , Front.—The pre of a laminar front may be derived 
directly from Navier-Stokes equations if the following or 
Curvilinear are that is, the front be 
tapering. (This condition obviously does not hold right at the tip, but the region 
of marked curvature 2 D back from the tip. 


«Fig. 2) hese? 


4 3. ‘The velocity distribution is similar throughout the front or, alternatively, 
the effects of non-similar velocity distributions small. (This obviously 


ignores the process whereby the boundary layer is developed. awl > 
‘For tw two-dimensional flow of an incompressible flow, the Navier- —- 


In the the af advance of tho frantic rolotively 

| 
q 
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a 

irs _ Since the flow is very nearly parallel to the bed except in the immediate vicinity ye . 
cle ___ of the tip, the vertical component, v, of the velocity vector is small compared to — 
horizontal componen ermsofEq 4 are therefore small excent 
— 
hy | 
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— PIG. 2.—TIP OF LAMINAR FRONT 


reduced to a statement of the hydrostatic pres-_ 


sure law. ‘(This is the same assumption stated above that curvilinear effects | 


_ Eq. 3 can now be rewritten as follows: 


In selecting the velocity dhetrictich: tant is assumed to hold throughout the a 


ie front, it is logical to choose one that is applicable to the uniform flow up- > 
Stream of the front. This well-known distribution,3 which can be derived from i 
tz Navier-Stokes equation for steady flow down a plane, is given by saab 


and at t = 


The function, = f(x), simply describes the pe of the surface of the con- 
_tainer holding the liquid prior to its release into th the flume, — 
Substituting Eq. 6 into Eq. 5 and letting = 1. 


My 


This the transient phenomenon that precedes the establish- 


- ment of the terminal shape of the front. A stepwise solution of Eq. 8 is now 


| being attempted on a high-speed digital computer to obtain the shape of the fl 
quia. and velocity of advance in the first few moments after release of the . 


a oh After the front has acquired its terminal shape, the derivatives of U — 
"art respect to x and t become very small as a consequence of the long, gradually —_ 
== profile. In order to arrive at the shape of the front, these derivatives 7 
are Eq. 8 along with the term involving (ad/: x 


3 Te Chow, McGraw-Hill Co., Inc. , New York, 

1959. 
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normal depth, D, for laminar flow down aplaneis™ 


This can be by separation of giving 


(5) ) 


+ . (12 


Ines the constant of integration it is more convenient to 


variable & (Fig. 1) so that the shape of the front is defined in terms a the dis- 
tance back from the tip. Eq. 12 can be then er in the form 
tan 


At = 0 so tha 


Finally, th the terminal shape of a gradually tapering two- lami- 


2 _ Experimental Investigation of a Laminar Front.— -In order to check validity — 
of Eq. 2 and 14 a series of tests was made with oil in a glass-sided flume as 7 
shown in Fig. 3. The oil used was a Mobil Oil Company product called Ambrex — 
710 with a kinematic viscosity of 8.5 x 10-4 ft2 per sec at 70°F. and 5. 8 x 
10-4 ft2 per sec at 80°F. and a specific gravity of 0.85. ie et 
Meg Prior to each test on the advance of the front, a test was made to ndiaun a i 
the normal depth and to set the adjustable weir in the headbox at the precise i: 
level of the liquid surface. The advance of the front tests were started with ao 
= _the quick-acting gate closed and the quick-acting bypass valve opened suffi- @ a 
@ | ciently to maintain the level of the liquid in the headbox at the previously set a 
=i ‘jgeveiion of the weir crest. The bypass valve was then suddenly closed. A link- _ 
age connecting the bypass valve handle to the gate caused the latter to open 7 
_ suddenly at the same instant, permitting the flow to proceed down the flume. 
The velocity of advance down the flume was measured by means of a peal 


were triggered simultaneously with an electric relay Scales were 
_ fastened | to the outside of the glass wall of the flume in the field of view of the 
camera. By reading the photographic negatives with a 30-power microscope, 
& accurate determinations of the depths (within + 0.0003 ft) could be ee 
_ ‘By keeping the depths shallow, both the terminal shape and velocity could 7 
be obtained for several slopes ar and discharges. In some instances, , the cameras 
_ were moved closer to the upstream end of the flume to cain the shape of i 
front a during itstransient phase. 
_ An auxiliary test was made to determine the surface velocity. This | informa- 
tion was obtained by timing the motion of a small surface float that was first a 


behind the men. caught to front and was as entrained ained 
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q recorder as the front passed each 2-ft station. The shape of the front was 
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different slopes are reported herein as shown in Table 
= _ The data on velocity of advance are plotted in Fig. 4 for the four runs. Notice 

; Ef that the velocity of advance in all cases becomes equal (or at least within the : 
range of experimental error) to the average velocity in the upstream uniform | 


_ The data on surface velocity are also plotted on Fig. 4. Notice that this 

not decrease significantly, thereby justifying the assumptions 


in the derivation of Eq. 14 that 9U/8x and 9U/ét are small. 


E _ The data on the terminal shape of the front areg given in Fig. 5. The fit of 
the data to the curve defined by Eq. 14 appears to verify the theory for the a 
_ terminal shape of a shallow laminar front. However, only a limited 1 of ' é 
TABLE _—LAMINAR FRONT EXPERIMENTS WITH OIL 


Tempera- 


0,838 
0.794 


TURBULENT FRONT 
of Equation . Shape of Front. the basis of 
7 the laminar front, the following assumptions are made for the e derivation of an 
equation a ‘Shallow: turbulent front: 


4 Surface tension and air can be ignored. 


Curvilinear effects are 
. The surface: velocity remains nearly constant throughout the front so 


4, The velocity distribution is similar throughout t the front, or, alternative-_ 
ly, the effects of changes in velocity distribution are small. (Here oe the 


‘On the basis of these assumptions, the s forces on an element in the front can 


ai assumed to be in Thus, referring to Fig. 6 


— | 
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inft |Vinfs} R | fF Gf 
— r depths, the flume was too short to obtain 
| the depth was tested. At impaired, 4 
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If the von Karman- Prandtl logarithmic velocity | distribution for rough open — 
channels is assumed to hoid throughout the front, then 

Gabris Keulegan4 has found from an analysis of Bazin’ wide e chan- 


Now in the upstream section where the depth is constant 


be 


It follows, then, 


. 19 into Eq. 16 
ale ao 

“s Ea. 20 has been integrated numerically with an IBM 650 computer for values” 
ot D/k = 10, 102, 103, 104, and 105 for k<d<0.99 D assuming that § = 0 when 
d=k The results are presented graphically in Fig. 7. A recomputation with - 

= 7, 25 did not produce a significant change in the shapes of these curves, 
In attempting to develop one universal profile for all fronts, it was noted 
that each of these theoretical turbulent front profiles could be transformed to = 
fit the laminar profile quite closely by asimple change of scale of the ae 
for each value of D/k. Thus, it follows that 


wat d 
ia) o(2)= tanh mAb 


“The function, ), however, canbe rewritten interms of the Darcy- ‘Weisbach 
coefficient, f, or Manning’s n, For one particular value of ‘Manning’ S| n, say 


on = ng, the turbulent profile and the laminar profile will be identical. = 


in which ch no and b both require experimental determination. 


- Experimental Investigations of a Turbulent Front.—The tes test p: program dis- — ® 
_ cussed in this section is part ofa larger research program intended to develop ; 


_ 4 Laws of Turbulent Flow in Open Channels,” by Gabris H. Keulegan, Journal =f 


Research, U. S. Natl. ‘Bur. of Standards, Research Paper-R.P. 1151, Vol. 21, , December, 
1938, 738 
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a criteria for irrigation The phase is 
an investigation of shape of the advancing front in shallow overland flow, which 
characteristically unsteady and nonuniform in its initial stage. 
___ In order to achieve control of variables not possible under field conditions, 
this test program is being conducted in laboratory facilities in which flow 
characteristics can be observed and recorded. The principal facility is a 
< = flume (Fig. 8) 69 ft long and 3 ft wide with appropriate equipment 
for control and measurement of rate of flow and channel slope. The flume is 
- 7 formed from steel plate with observation windows inone side, and is supported 
by large trusses to minimize deflection. The floor of the channel is covered — 
with one inch of concrete which has been sanded to pro-_ 


vide a smooth, uniform surface, and has with a coating 


Because local disturbances such as and virtually preclude 


Ny conventional laminar characteristics in atypical irrigation stream, these tests 
_ are being conducted inthe turbulent range of flow. Values of the variables were © 
- selected as being appropriate to both border irrigation and laboratory investi- 
gations. ‘Tests have been made for nearly all of the 


&g 
— 
8.—LABORATORY FLUME USED IN THE STUDIES OF 
a 
il 

> 
— 
— 


i di t foll 
n ical ed ed as 


1 Channel slope _ 10%, 0.25%, 0.50%, 1. 
oa Unit stream - 0.033 cfs, 0.10 cfs, 0.167 
rate Zero for these tests _ 
4 Channel (1) n = 0.0081 D-9.031 - 
(2) n= 0.0119 D-9.247 (metal 
(3) n= = 0.0042 D- “1. 37 (crushed rock ~ 
convenience, surfaces used in these tests are described terms 
of Manning’ sn for the range of depth involved in this study and for steady flow 
conditions. The equations for Manning’s n as a function of the depth were | 
determined with tests on uniform flow prior to the tests on advancing fronts. 
¥ These equations are valid, of course, only over a limited range of the depths. 


— a Surface (1) is the sanded concrete floor of the flume with a plastic seal 


ing ath (1/2-in. openings, 5 Ib per sq yd) over the smooth concrete. . Surface 
(3) achieved by placing a single and complete layer of crus:._d rock 
; (3/8 in, to 3/4- in. sieve size) over the expanded metal lath. Water tempera- 
: i: ture for the majority of these tests was between 60°F and 65°F. Test conditions : 
“a thus encompassed a ten-fold increase in channel slope, a five-fold increase in 
. stream size, and approximately a six-fold increase in surface roughness as q 
characterized by Manning’s n for average depths ™ 
3 : Prior to each test of the advance ofa front over a given surface, the channel _ 
slope ‘and rate of ‘flow are established. During the time of these adjustments | 
_the entire stream is by-passed around the channel and back to the sump. The 
test is started by closing a quick-acting valve, causing the water to filla nar-— 
_ row secondary forebay and then to flow through dampening screens into and 
‘ bak down the channel. Full flow into the channel is thus achieved immediately fol- 7 
} lowing water entry, and depth in the forebay ‘may increase with time if down- Z 
a _ When the advancing front ‘reaches the end of the 50-ft test section in the 
flume, four cameras photograph simultaneously the water levels in each of - &§ a 
12 open-tube manometers connected to the bed of the flume, The depths are = : 
- gubsequently read from these photographs with microscope aid to the nearest 
1/1000 ft and adjusted as necessary for channel slope, etc. These depths per- 
mit the drawing of a best fit water surface profile for each test. Points from 
curve are used in the subsequent analysis. 
_ Experimental Results on Turbulent Fronts. —Typical water surface profiles 
are contained in Fig. 9 from which it is obvious that a, d, is a function of | ; 
_ p at least the distance behind the front, x; channel slope, sin 9; unit stream size, 
_ and the ; roughness of the surface. Notice that all the profiles approach uni- 
_ form flow conditions except the 0.0010 slope in Fig. 9(a). The seme was not 
a long enough for this front to develop its terminal —. gf 
In order to evaluate ng andbinEq. 21, fronts were investigated whose pro- 
_ files fitted the theoretical profile for laminar flow. Certain tests 1 with water — a a 
flowing over the metal lath (surface 2) satisfiedthis requirement. These tests 4 - 
established the value of np tobe 0.020. When the adjustment for relative — 
ness suggested earlier, (n/ng)>, is applied tothe abscissae where ng is takento 
_ be 0.020 and b = 2/3, asingle curve as shown in Fig. 10 results, indicating that ‘ 
_ the terminal shape of turbulent fronts, at least for the flow regimes tested, can > 
by the use of Eq. 21. Fig. 10 is is useful in 


HYS 
— 
cin 
ae 
i 
| 
| 
Ta 
al 


| DISCHARGE 


SURFACE |(3) CRUSHED ROCK 


% 


is 


033 


: 


— 
as 
; 


gue 


LATH 


3- CRUSHED JO 
ROCK 


URBULENT FRONT 


‘FIG. 10 


FIG. 10.—TERMINAL SHAPE O 


~ 


of 
| Ff 
= | 
— 
— 
— 
ae 
— 
ey 
— 
— 


a ‘September, ber, 1961 


the length of the front. It will be observed that the relative depth has attained _ 
es value of 0.99 when the length parameter (abscissa) has obtained a value 
It now appears to ‘predict the water surface profile for turbulent 
- conditions and to determine the distance behind the front beyond which essen- _ 
tially uniform flow conditions exist. These results, however, are for relatively 
shallow fronts, Short fronts may not develop their terminal and could 
_ be expected to deviate from the shape presented here 
should be noted that in addition tothe work on the turbulent water 
shape reported here, investigations are under way on the following phases of : 


a. surface. shape where extends throughout the entire 


surface shape where infiltration is and of the ‘the form. 


k t@ 


When a liquid is introduced slowly into it decolerates and 
quires” a long, gradually tapering profile | which becomes stabilized and there- 
after advances at a constant velocity equal to the average velocity in the up- 
a The terminal shape of the front for ial flow is a function only of the — 
— slope ‘and the normal depth corresponding to the discharge | and slope. For =: 
F bulent ee the shape is also dependent upon the se roughness of the 
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A PPENDIX. 


= constant of — 


= Darcy- -Weisbach resistance coefficient; ~ 


ya = = Manning’s n for a particular turbulent Profile (No =O. 02); 
pressure; 


7 
per t width of flume; 


= Reynolds Number Vv D/y; ¥ 


= component of velocity parallel to the bed; 


= terminal velocity of advance of the front; 
= = average veenenche of flow in uniform section; i, 


= body force per unit mass in x- direction; 


* coordinate axis parallelto the bed ( waatiies inthe downstream direction); 


body force per unit mass in y- -direction; 
— y aa coordinate axis per pendicular to the bed (positive uy upward) ; 


7 “di a axis wi i at ti f the f 
ey in at ante-with origi nz e ront (positive 


shear on the bed in the. 


70 _ = shear on the bed in the upstream uniform sec 
= kinematic viscosity of fluid; and 


im 
 ‘=normal depth of flow in upstream uniform sec 
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Ottis: 


_A liquid jet deflected by a device that is composed of straight seg- 
ss but, if these segments are short compared to the thickness of the jet, a 
the deflection angle is unknown. It can be computed, however, for the two- q 4 


dimensional. case for symmetrical deflectors ‘by applying free streamline 


_ The segmental deflector is first mapped by conformal transformation onto 


integrals of are > then evaluated with a high 

_ speed digital computer for a wide range of the dimensionless parameters -_ 7 
fining the geometry of the he deflector. The of this analysis are pre- 
in graphical form. 


ways of dams are than compared to this theoretical analysis and it is shown © 
that the exit angle of the jet can be accurately predicted by the theory, despite ~~ 
the fact that boundary layer effects and the influence of gravity on the flow in 


“the deflector are ignored. 


For the application of the theory to flip buckets on spiliways,a set of curves’ 
is presented for designers. | This set of curves gives the dimensions of the 


_ Note.—Discussion open until February 1, 1962. To extend the closing date one month, *) 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
_ of the copyrighted Journal of the Hydraulics Division, = of the American So- 
of Civil Vol. 87, No. HY 5 1961. 
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4 
that ; are required to produce ‘the desired angle as a function of the 
_ face angle of the dam and the thickness of the nappe. Thus, flip buckets may Pe : a 
be designed to give the desired exit angle of the jet fora particular ein - 


— shown on Fig. i, , could be used if the exit angle of the flow could be pre- oo 
for each ‘discharge. problem arises of determining the exit 
in- 


angie the spillway face, the dimensions “a” and “b® of the buck 
et, and the inclination angle a» of the lip of the bucket. If the exit angle is 
known the trajectory of the free jet can | be be readily « determined by the laws for 
_ In the absence of gravity, the flow would leave the bucket and attain an nangle — 

Bas shown by the dotted lines in Fig. 1. To find this angle it is assumed (a) a 
_that the influence of gravity is small within the confines of the bucket (it is, - ; 
course, significant when the flow leaves the bucket), (b) that the influence of 

° _ viscosity is negligible within the bucket, and (c) that the flow approaching the 

_ bucket has essentially a uniform velocity distribution or, alternatively, that az 

- the non-uniformity of the velocity distribution is not significant. The first as- 
sumption, of course, restricts the application of the theory to high velocities _ 

-atwhich inertial forces are to forces (i. e large 


flip bucket on the spillway of a dam generally has atypical ‘ski jump? 


i A values of the ow number, , V/Vgd). The second assumption requires ene 

ia inertial forces as compared to viscous forces (that is, large values of Reynolds» 

number Vd/ v). Both of these conditions are usually satisfied on the spillways 
s all except | very low dams. As far as the third assumption is concerned, it . 
is known that a non-uniform velocity distribution exists in the approach flow 
but the area of large velocity gradients is confined to a small region near the | 


- of the spillway despite the fact that the boundary layer extends to the free 


surface. It appears, therefore, that none of these three assumptions seriously 91 


disagrees with real conditions. The problem | is thus reduced to a two- 


dimensional potential flow problem to which classical hydro-dynamics canbe 


_ A general theory is developed in the following section for the design of jet 
iacers composed of straight segments without particular reference to spe- 


cific applications. The analysis follows directly from the application of com- 
on variables, in particular conformal transformation theory, for which the : 


7 
reader may refer to many 


CONFORMAL TRANSFORMATION OF THE FLOW REGION © 


Notation .—The letter adopted for use in this paper are defined 


"where they first appear and are arranged alphabetically, for convenience of ref- a 


erence inAppendix™ 
_ The boundaries of the flow region under consideration are composed of two 
types, ‘straight segments and free streamlines. ‘Si Since along the straight seg- 


"ments a and free streamlines the inclination and magnitude, respectively, of the 


a closed polygonal figure whose horizontal portions correspond to straight ; 


and whose vertical portions correspond to free streamlines. ‘This 


polygon can then be mapped by the Schwartz-Christoffel transformation onto ae 


= line where the complex potential can readily be determined, rT 


st 


results in 


the and inclination of the vector Now define 


by Lamb, Sixth Edition, Dover Publications, New York. 


: a 7 Fluid Dynamics, by V. L. Streeter, McGraw-Hill Book Company, Inc., New York, - 
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anoles and are taken and denoted by a. Letting w be the poten- 
— 


— 
and flow re; region of the z-plane (Fig. 2) onto indicated region on of 
- Q-plane. This being a closed polygonal region the Schwartz-Christoffel trans- 


formation can be applied to map the boundary onto the real axis of the t-plane. 


sis 


e 


to the syn symmetry introduced by making ay , and as equal, 


i, 0, +1, may be assigned to the points E, P, and C, respectively, and -k and 
+k to the and of D is taken to bef. 
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es THEORY 
-¢, 


woking the principal value of att= N=O 


-isin- 


From the source and equal sink at t aed and t = f, respectively, the following 
complex potential is obtained: ‘ 


oie (in (t 1) 
the. definition of Fa ‘Eq. 8 is ‘now 


Vv 


— 


 eangeating from A Bee and noting that in this range of t 
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— HYS 
J GI) = 2 a 4 
— 
4 
: = am 
— 
@ 
a 
— Similarly, integrating from B to C 
lz 


so that the expression in 1 parenthesis is now real, , Eq. 13 is is cotained af 
+ (f 


+ 


By assigning v values to f,and it is now now possible to ‘compute ute correspond - 


that the ordinate scale is for each 


FIG. 3.-THEORETICAL EXIT ANGLES AS A FUNCTION 
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integrals, ‘as descirbed in Appendix I, 


Eqs. 11 and 13 have been evaluated with an IBM 650 lesa 7 eae 
a = 30°, 45°, and 60°. For each of these angles specific values of 8 were se- 
lected, in “each case less than | a, and the corresponding values of the dimen- 
BS sion parameters a/d and b/d were computed, ‘The results of this computation — 


| 


AB= 


al 
w 
a 
w 
= 


4 


_ the exit angles of ‘the flow for each bucket for ‘various aa of the | approach — 


flow. The data from these six test series are plotted in dimensionless - _ 


5 
partial fulfillment of of the for the of of Master of Science, 
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4 STREAMLINE THEORY 
The exit angles were at first measured directly. 
. greater success k by measuring the trajectory of the lower nappe of the "jet and ae Fi 
the exit velocity a and | computing the exit angie to produce the tra-— 


a cee comparison of Copp” s data with the > theory is given in Fig. 5. Because wiih c 
the limited number of tests, only 9 points could be directly compared. The 
maximum deviation is 2.5° | and the mean of the deviation is less than a ear. 


“ter of a degree. Copp states that these are within the range of ex- 


In the obvious application to flip buckets ber an example of a seg- 


mental bucket see Noxon Rapids dam9), the angle § is first determined on the — 
basis of a safe trajectory length, L (Fig. 1), using the velocity V as calculated 
for the available head H: The dimensionless parameters a/dandb/dare se- __ 
oe lected from a graph similar to Fig. 3. The dimensions of the bucket are then _, 
fixed by choosing the most critical discharge with regard to scour and com- >) 
a puting the corresponding depth d. _ Once the bucket has been thus designed for 
a certain discharge, the angle Bal and the length L must be recomputed for all 


wou 


A numerical solution of the integrals of transformation shows effect 
of the ‘geometry of the deflector on the exit angle of the jet. 
3. Data from previous experiments on flip buckets for spillways ‘show that 


. ae The curves in Fig. 3 relate ‘the quantities a/d and b/d for fixed angles a we 

B. They were obtained by choosing values of k =0. 0.06, G.40, 0.20, 
0. 90, the of Ea. 18, evalualing the 


P 
4q method has been developed in the foregoing sections for designing seg- 
eg 4 mental deflectors for prescribed values of the angle a and the exit angles of _ “ 1 
— 
_ 
— | 
The flow ir yments can = 
— 
«=e 
— 
The intercepts of the curve g to a/d 
puted from the separate formula 
9 “Spillway Model Study Noxon Rapids Hydroelectric Development,” = | 


aid 


s cos 8 in -- 


te Be 
B) 


16 was obtained ‘i Siao and P. G. “Hubbard, 10 M. ASCE, ‘using ar an 
= = argument very similar to the conformal transformation presented previously. A 
_ The numerical evaluation of the integrals Eqs. 11 and 13° 13 is complicated b nd 
4 the fact that the integrand of Eq. 11 is singular at both end- points of the inter- 
of integration and the integrand of Eq. 13 is singular at the lower end- 


— ee for values of k approaching 0 or 1 the integrands develop 7 x 


evaluate integral Eq follow — and write 

of 


(t+) @- all 


k)1-(a/m) y (x) + 


10 “Deflection of Jets,” by T. T. Siao and P. G. Hubbard, Free Streamline 
" of Transition Flow and Jet Deflection, ed. by J. S. McNown an iu-Sheen Yih, Studies a : 
in Engrg. Bulletin 35, State Univ. of lowa, 1953.00 
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2 For small values of k eantatass are encountered when using an integration | 2 
scheme involving starting values of at k-h, Consequently, the writers have 
chosen to integrate Eq. . 20 using the Runge-Kutta Fourth Order Method which 


involves no special starting formulas. én (For a discussion of this method, see 


ak Eq. 2la gives th the initial condition and Eq. 21b can be used to obtain y (k + h) 


s without the indetermination that would arise from a straight forward applica~ 


tion of the Runge-Kutta formulas. Thus denoting k + nh h by ‘Xp an and Ly ¥ (Xn) by tay 


+ 


a tiating Eq. 20, upon set- 
‘ 
* | 
| 
| 
* 
F (x + (1 - a) (v, foo) |. — 
12 Numerical Solution of Differential Equations, by Miine, vonn W ey & 
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_ Notice one of objections to the numerical integration of 


5. the preceding formulation, since the complicated evaluation of the function 
F (x) need be done only twice for each forward step of the integration. An equiv- a 


alent of can be y ( x) = (x -k k) w (x). 


(x) = x =) 
the differential Eq. 20 becomes" 
a 


with the initial condition y (k) =y 


Ge: 
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wt 
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=— |F (x) + 
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to 1 since error in the evaluation of F (x) is minimized by Reet by h B« 
and the necessity of F' (k) is eliminated. a/d is eval- 


1. 


_ These two integrals can now be evaluated by the same technique that — used 


with 1, 


2 ‘length of segment of 


length of ‘downstream m sloping segment of eer; 


- depth of approach fl flow; 


base of nat 1 ithms; 


dinate in t- amend corresponding to infinity; 


— 

a 

— 
— — 
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g = acceleration of gravity; 


= coordinate in t-plane to ends horizontal ‘segment o of 
flector; 


| 
= horizontal distance of of trajectory; 
3 N = = - constant of integration; : 
variable e of integration; in plane; 
= velocity c in x-direction; 
Ve = velocity vector; “ 


v = = : velocity component in y- -direction; 


x = real axis in Z- lane; 
y = imaginary axis in z- -plane; 


5 
= complex coordinate in Z-plane; 


a = angle ¢ of inclination of deflector faces; 


ar... 
pe exit angle of of jet at at infinity; 


derivative of of complex ve velocity potential; 


= inclination of velocity vector; 


= kinematic viscosity of fluid; 


- we complex velocity potential; a 


= complex coordinate in Q- -plane. = 
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‘SYNO PSIS 


A open channel was investigated to determine 


the characteristics of uniform flow in an idealized highway gutter. be nt 

"velocity contours, and transverse velocity distributions are presented for var- 

ious slope configurations and flow rates, 


INTRODUCTION 
Inc with. an ananalytical investigation of flow gutters: into curb- 
opening inlets, 2 extensive experimental work was done with uniform flows in 
a hydraulically smooth triangular channel (one side of which was nearly ver- 
head placed on a continuous grade with no depression in the channel floor (Fig. 7 
= The data taken were almost entirely limited to flows at supercritical velo- — 

' cities because the theoretical development was predicated on such a = 
___ Notation. .-The letter symbols adopted for use in this paper are defined 
ea they first appear, in the illustrations or in the text, and are —— 

"alphabetically, for convenience of reference, in the Appendix, 


_ Note.—Discussion open until February 2 1962. To > extend the c closing: date one one 2 month, : 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
4 the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
= of Civil Engineers, Vol. 87, No. HY 5, September, 1961. a ves ; 

1 Research Asst., Stanford Univ., Stanford, Calif. 

2 “Hydrodynamics of Flow into ‘Curb- -Opening Inlets,” by R. J. Wasley, thesis pre- 
"nel to Stanford Univ. at Stanford, Calif., in November, 1960, in partial fulfillment of — 
the requirements for the degree of Doctor of Philosophy; also available as Technical © 
Report No. 6, Dept. of Civ. Engrg., Stanford Univ., Stanford, Colm, 
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‘TRANSVERSE SECTIO 


“SECTION B-8 


ors 


INATE SYSTEM 


«SECTION 


COORDINAT 
TRIANGULAR CHANNEL 


of 
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G. 1.—PERSPECTIVE DRAWING 


SECTION A-A 
STILLING RESERVOIR SECTION 
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- : oe Equipment. —The channel as finally constructed was an idealized full-scale 
prototype, 50 ft in overall length by 6 ft in width, _ ‘The actual triangular chan-— 

_ nelwas 32 ft-3in. in length, the remainder being the curb-opening inlet (Fig. D 
a: The rigid frame construction was achieved with 1-1/8 in. plywood panels ~ 
bolted to transverse 315.7 beams on 4 ft centers. The plywood was surfaced _ 

_ with 1/4 in. tempered masonite fastened with waterproof cement. The joints 

= were filled with a waterproof sealer, sanded, and four coats oh a waterproof a 
7 epoxy resin were applied, | each coat being w wet- -sanded, The resulting surface 
was almost perfectly smooth visually and remained so throughout the experi- _ 
‘The I- beams were fastened t to threaded 3/4 in, standard galvanized pipes” 
by means of nuts and washers on top p and bottom of flanges welded to the bottom _ 
of the beams at each end. The pipes were attached to 3/4 in. steel rods set in . 
the concrete floor of the laboratory with iron cement. A bracing system of 
tension rods and turnbuckles was added to the pipe ‘substructure for rigidity - 
The stilling basin at the upstream end was a reservoir 10 ft by 8 ft by 4 ft 
in height. A baffle system was utilized to quiet the flow. The entrance of the 
water to the channel itself was effected through a stainless steel “transition 
section. A rubber diaphragm connected this section to the reservoir and allowed [ 
any changes in slopes to be easily accomplished without causing leaks. == 
A surveyor’s level was used to accurately position the apparatus. Trans- 
4 7 __ verse slopes, 6,, were possible within a range of zero to 10%; the range of © 
slopes, , was nearly as flexible, with a range of zero to ) 8% 


_ posed of dots at one- quarter foot intervals painted on the channel surface. The 7 
origin and coordinate system are shown in Fig. 1 
i. horizontal datum was used for the support of the measuring apparatus by 

"fastening I-beams to the pipe s standards at an arbitrary el elevation from 
floor. Aluminum tees were bolted to the beams and an aluminum cross- 
-member was placed longitudinally so it could slide transversely on the beams, | 
To this was attached either the point gage or the Pitot tube, , capable of moving 
longitudinally on the cross-member. eer wads 
_ Certain appurtenances were used waves 
ances that developed. An adjustable sluice gate was used in a transverse po- ; 
sition at the entrance, consisting of a series of sheet aluminum panels, each of © 
_ which could be placed at any desired height above the channel. Slightly down-— 
stream, two groups of aluminum vanes, one behind the other were positioned : 
to straighten the flow and lessen or eliminate any surface irregularity that oc- 
In addition (for certain flow rates and configurations), to aid in =. 
omen action, a screen was placed in front of the vane system, and a piece 
>. of cloth was attached to the downstream end and allowed to ride on the | surface 
of the water; it was determined that these latter modifications were not pa: par- 
= —, effective (Fig. 4). it was found that the lower flow rates did not re- a : 


Th 
300 hr of Ac gi. noticeable had in the surface. 
: _ During changes in configuration, some slight cracking at the joints occurred, 
but this was quickly and easily repaired with no ; appreciable effect on any m mea- 
surement, Toward the end of the experimental program, some very slight sag 
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45 
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between became noticeable. However there was ‘no 
Flow Measurement.—The desired flow rate Q, was obtained from readings | 
_of a manometer that had been previously calibrated by standard weighing meth- — 
ods. Two Venturi meters were avail25ie for use in this experiment, one with 
a4in. throat diameter and another with a 2 in. diameter. The total range - 
_ flow rates in the experimental program was from a minimum of 0.0032 ee . 
ae a maximum of 2.98 cfs. _ The measurements below 0.040 cfs, however, were x 
"obtained | from the actual weighing ofa quantity of the fluid ina specified in- 
J _ Depth Measurement .-A point gage was used in all the measurements of 
pe _ The point itself was constructed with a 45° angle bend near the end to 
allow easier measurements at or near the curb face and directly under the — 
i _ Supporting I- beams. _ The apparatus was suspended from the datum plane : 
_ described earlier. The difference between the channel elevation and the water 7 
an surface at a given point constituted the depth at that point. —™ 
ae It was found that any depth measurement could be satisfactorily ie 
to 0.001 ft with the apparatus previously described. Although there were slight _ 
“rapid fluctuations on the water surface itself, probably caused by turbulence, 
_ these could be averaged quite adequately by eye and identical measurements = 
- _ Velocity Measurement. -A pitot probe with a vertical tube was used for ral 
velocity measurements. The velocities were generally great enough to create | 
an adequate velocity head to insure sufficient accuracy in reading. Three dif- e 
_ ferent diameter probe openings were used (0.010 in., 0.025 in., and 0.065 in.). J 
_ Itwas found that all three probe diameters gave essentially the ‘same readings; 
» _ hence, it was decided to use the intermediate opening size (0.025 in.) to reduce — 
the necessary response time compared with the fine opening size. 
cs 
In addition, an internal check was used to allow a comparison with » the mea- 
sured flow rate. This was achieved by the integration of the transverse files) 
city distribution (initially obtained from the integration of the velocity profiles) 
to determine the quantity of flow. A reasonable agreement was obtained for 
oa computed flow from the Pitot measurements that generally gave —— 
higher results over the measured flow (see the section entitled “Results”). 
7 The capillary | correction was found to be negligible for all lym ~aepesorell 
eo; The most critical possible correction was connected with the fine probe — 
opening in which the internal diameter of the vertical standing tube was 0. 015 
oR (both the other standing tubes were larger in diameter). The capillary rise3 < 
= for this diameter using tap water is about 0.008 ft. The minimum total ey 
oe - reading taken was 0.280 ft (accurate to about + 0.004 ft), causing the agueey 


‘This section presents | the actual experimental data in reduced form, where- 
as the v: various comparisons between this investigation and experimental work — 


from different sources are given in the section entitled “Results.” ~~ 


=. “Handbook of Hydraulics, ” by H. W. King, McGraw-Hill Book Co., Inc. ., New York, a 
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the maximum 1 transverse slope distance, and ‘the depth of flow a at 
- face, Zo are given as functions of the flow rate, Q, in Figs. 2 and 3. A a 
verification “of steady uniform flow can be seen by ‘noting the slope of the re- a 
| 


— line. Manning’ s formula a triangular c channel of this type in this” : 


J 
‘FIG. 2.- VS yo FOR ALL CONFIGURATIONS 
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ye dove 


FIG. 3.—Q VS z, rom - 


‘nee (neglecting curb shear) can be derived readily! and is found to on 


“Elements of R.K. Linsley and J. Franzini, McGraw- 
Hill Book Co., Inc., New York, 1955. 
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& line indicates the uniformity of flow conditions, 


ene, f, as expressed by® 


inwhich g is acceleration — to iis is the channel hydraulic radius, and s 


a is the average channel velocity; ; this equation assumes that the 
The Reynolds rR, is given 


in which p is the kinematic viscosity of water. The data and the Stanton diagram 
will be given in the appropriate section under “Results” and the latter Aten ad - 
>. The experimental velocities are plotted in the form of contours for each ach : 
_ configuration and flow condition investigated (see Figs. 5 through 9 in which 
_ the vertical ‘scales are greatly exaggerated). _ The contours are drawn inter- 
polating among average depths and velocities at a given position, obtained from 
_ superposing the various profile sections. The error of averaging the veloci- 
ties and depths at a given point or section was not significantly different from 
any - extreme value for that point or section. The superposed velocity. profiles 
fort two typical situations = 1.00%, 99 = 1: 12 ,andQ = 2. 98 cfs; = 5. 


- the averaging of the panyeryaryee and depths at a 1 given | position, can be : ‘accom- 


_ The various x-coordinates of the transverse sections s (applicable to all ) all con- 
_ figurations and flow rates) are given in the figure. 

Surface tension 1 phenomena were observable throughout the investigation. 
To insure consistency of these effects, the surface was completely wetted be- on 
fore each run and allowed to dry before taking any values. The increase in the 

_ measured transverse plan dimension, y, after this was done and after flow was 

et established was attributed to surface tension causes. For the flat cross slope, a 

2 this difference was noticed to be as much as 20% for yo. For the steep cross” 

- As slope, this difference was almost negligible, with gravity being the stronger 

= influence. a The longitudinal slope had relatively little effect on these changes. a 

From Fig. 2, it can be seen that there is a difference between the calculated 
Se = 9 0)» shown by dashed lines, and the experimental y,, shown by sol- 
id lines. For all configurations and flow rates (except for very slight opposite | 
4 results for 9 = 5%, 9, =1: 12), the e calculated values were greater than the ex- “i 


perimental. it is to expect that surface tension between the liquid 


“Elementary ‘Fluid id Mechanics,” by K. Vennard, John n Wiley and Sons Sons, New 
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air interface would exert a “pull” te ‘toward the curb causing this to occur 
7 simila ilar to to the way a droplet of water forms on a solid surface), 


ona triangular channel 
= similar shape (except with the curb slope shallower than 1 on 2), but with | 
a concrete surface, that the distribution could be expressed 


7 Cou 


sinc 


in which y y is the particular transverse dimension, oi and C$ are weragevelo- 


city in in any given t transverse se section. — 
= this does not fully represent the shape . of the actual distribution. as 
_ Another method that theoretically describes the distribution without any co- ; 
- efficients, except the utilization of the Darcy friction factor, f, can be given. 
"experimental information (or handbooks) an average friction factor 
= the channel can be given by Eq. 2. - All these quantities can be readily | 
calculated from the known upstream parameters 
a certain on this transverse section, and 


Channel, ” by C. L. Larson, thesis presented to the Univ. of Minnesota, at Minneapolis, 
Mi nn. in 1940, in partial fulfillment of the requirements for the degree of Master of 


i _ 6 “Transverse Distribution of Velocity and Quantity of Flow ina Shallow cali, 
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Comparison of Investigations.—On the basis 
— of each profile section (originally obtained fro po 
Fig. 4), transverse velocity distributions c 
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with the we substituted for | the R. 
_ Because of such a simple representation of the distribution as_ Eq 2. 5, a4 
; ad be expected that complete agreement | with | experiment would be more _* 7 


that consequence of the theory. It is when compared with 


th 


=(2.75K0. 


TRANSVERSE ve LOCITY DISTRIBUTIONS FOR 


the results of the > Pitot measurements, th: that the shape of the distribution is: 
quite well defined, in general, but the ‘magnitude is somewhat in error nearer 
the curb (assuming the Pitot results arecorrect), 0s 


- It is noted that Eq. 5 does not account for the curb shear that would cause * | 


velocity there to approach zero, This could be accomplished by represent- 


aad by a one-half power expression similar to 


— 
— 
— 
1a = 4 
— 
3 
op 
a 
a 
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the section and would curtail the direct this method. ‘This 
, _—— of the distribution is considered below in this section, in connec- 
tion with another topic. this constitutes only a 
It could be argued, that because the friction factor in Eq. 2 was oie ahd a 
by using the bulk quantities of the channel section (that is, an average friction 
factor) and then applied directly to Eq. 5, that fair agreement should be ob- a 
_ tained. ‘However, the variation of the friction factor across the channelis quite 
= ” and an “average” found in the preceding manner well be 
Another approach to describing the friction factor, f, can be 
bs ag vanced. _ ‘There should be a systematic variation of the friction factor To 


_ the flow, approaching a large value as the depth decreases toward zero. To Bah 
_ determine this, an incremental slice of fluid could be considered in the trans- a a 
ba verse direction and the experimental | mean velocity and depth (as the onal 
radius) in that element used to compute a friction coefficient for that small _ 
aoe width. In addition, the Reynolds number associated with this factor could be 
e *. calculated 1 from Ea. 3. This would | allow a curve to be plotted on the a 
From this the possibility of there being a definite ‘relationship 
.< could be investigated in the manner in which f varies with y or with R that 
_ would allow an analytical expression for | Cy in Eq. 5 as a function of either of - 
these quantities, Besides requiring more detailed information on the incom- 
ing experimental parameter, the end result from an approach of this kind would 
not be greatly different from the previously noted “average” 
_ The “equivalent diameter” (d=4 R) is also suspect. _ For smooth channels" 
a that are similar in shape toa circular pipe conduit, or even ae tri- 
; angular flumes not badly deformed, this works quite satisfactorily.7? However, 
for badly distorted channels, as in the present case, far from a circular cross- y 
- - section, and for smooth surfaces, the “equivalent diameter” has been shown to © 


not equal t to four times the hydraulic radius. 


a 4 thausationl representation of the distribution and the actual variation, it is only — 
necessary to integrate the theoretical distribution (accounting for the triangu- = 
LC “lare cross- section) to obtain the flow rate and compare it with the measured 
in the channel. _ This has been done and is given in Table 1. 
aha It can be noted that a surprisingly constant ratio of these two flow —— es 
(measured to theoretical) exists. The average value of this ratio, N, not con-— 
sidering the 1 on 96 | cross-slope, is 0. 856. Disregarding the longitudinal slope | + y 
_ of 5% and also the 1 on 96 transverse slope, the ratio is 0.855, The total mean - ; 
value is 0.829. The average value of the ratio for the shallower transverse 2. 
To account for the constant value of N, an of the va 


7 Channel Flow at Small by L. E. Silberman, 
and H. C. Nelson, Transactions, ASCE, Vol. 123, 1958. 
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If this is done, and the distribution again integrated, it is found that a certain 
portion of the total discrepancy between the theoretical and experimental flow 
_ rates can be accounted for. For example, in the intermediate configuration — 
(90 = 3%, 8 =1on 24, and Q=0.725 cfs), N is only increased from 0 0.848 to to 
873 (N equal to unity being perfect agreement). 
— - fact that the friction factor in a rather complicated manner, in n addition i + 
: > to the uncertainty of the value of the “equivalent diameter,” probably accounts i 
for: the major portion of the discrepancy. . However, as previously s stated, the 
representation has the tremendous advantage of being simple and directly use- 
f= able while being accurate to within about 85% for computing the true flow rate ae”; 
when integrated. The friction factor, for any given practical problem, , would - 
_ probably be ascertained from the estimation of the physical condition of the _ 
channel itself and would be susceptible to any error of judgment. — Hence, any <7 
refinement should be regarded in this light as it might not be - consistent with 
R. W. Powell, F. ASCE and C. J. Posey, F. ASCE® have done some experi- 
‘mental work that has some bearing on the preceding problem. They have in- — 


= vestigated th the variation of the friction factor in a 


Q, , (Pitot) in 


- 
cubic feet per feet per second 


0.29 


:96 


for various slopes, configurations, flow rates, and roughnesses. They” > 
x. , found, for flow at supercritical velocities in a smooth channel, the curve a 
- best fit lay above the omeeth pipe curve (using d=4 R). The Powell and Posey 4 
curve is seen in Fig. 15. Their curve is, , generally speaking, consistent with 
+ the experimental data taken in the } present KS nt ofvelo- 


To guarantee the avoidance of any gross errors in the ofvelo- 


“ical of Figs. 10 through 14 (using a large scale), the rates 

mes Pitot measurements were determined. The results are given in Table 2. _. 

is acceptable for the "higher flow rates but n¢ not as for 


“Resistance Experiments ina ‘Triangular’ Channel, by R. 
Posey, Proceedings, ASCE, Vol. . 85, No. HY 3, ‘March, 19589. 
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OPEN CHANNEL FLOW 
the lower flows. This could result from the fact that the vetostiies were — 
_ in the latter case and somewhat more difficult to measure | . with the ramen 
1. The main body of information is summarized in Table 1, as far as the im- 
portant parameters and values are concerned. Stanton ‘diagram, discussed 
earlier, is given in Fig. 15 using the data shown in Table 1. . In addition, other 
applicable information from the Corps of Engineers,10 Larson,11 and Powell 
4 and Posey, !2 is The resulting curves are seen to be of the expected 
oT _ Range of Measurements. .—A matrix of apparatus configurations was inves- 
_ tigated. ‘This information is most easily and clearly presented in a plot = 


as Fig. 16. The solid circles indicate that the measurements were taken for 
om 
| 
a s- -Corps of Engineers - 7 
*-Corps of Engineers - Broomed 
— Smooth pipe curve 


FIG. |. 15.—STANTON DIAGRAM 


| Ser % Ordinary hydraulic critical depth criteria were applied to obtain the 
= critical slope, (So) crit.- at using the expression that allows the calculation of 


“Airfield Drainage Structure Investigation,” U. S. Corps of St. Paul Dist. 
Sub-Office, Hydr. Lab., Report No. . 54, April, 


- Minn., in 1949, in ‘partial 1 fulfillment of the voquivemente for the degree of Master of 

eB “Resistance Experiments in a ‘Triangular Channel, i by R. Ww. Powell and C. J. 

Posey, Proceedings, ASCE, Vol. No. HY 3, March, 1959, P- 31. 
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_ By: wine a representative computed value for Manning’s n n of 0.009 forthe sur 


face of the channel (see Table 1), the various critical slopes were 


compared with the actual slopes, as given in Table 3. 


_ Although in all cases the critical slope, @o)er it. Was less than the actual &§ 
slope, %o, the former was quite near the latter | for the V/ 2% configuration, es- _ 
_ pecially for the lower flow rates. 
_ that raised the suspicion that flow at subcritical conditions was occurring due > 

vs to the lack of absolute e certainty inherent in the previously noted critical aia 7 


Experimental velocity | contours and traneverse velocity 
smooth, shallow, triangular open channel have been presented for uniform flow : 
‘conditions for a wide variety of slopes and flows at supercritical velocities. 
Velocity distributions were determined f1 from graphical integration of velo- 
“city profiles. These were represented mathematically using other experi- 
- mental results (Eqs. 4) and by a method based on Darcy’s equation (Eq. = j 
The latter seemed to coscre. e the shape of the distribution for | smooth par © . 


am 


> 


af ‘the Reynolds number plotted reasonably « on ‘the Stanton diagram rte. 15) indi- 
mt cated that the the application o of ee results to rougher surfaces would ras 


ry under whose direction this research was conducted. Several 
discussions were held with Carl F. Izzard , Fe ASCE, Chief, Division of Hy-— 
_ Research, Bureau of Public Roads, who also provided much useful _ 
information. assistance of F, G. J. Duncan, R. 
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equation in the f 
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_ It was shown that there was a fairly — 
rate to the semi-theoretical (Darcy) flow rate 
q 
— United States Department of Commerce, | 


September, 


ge The following symbols, adopted for use in this paper, conform ied 
with: “American Standard Letter Symbols for Hydraulics” ‘(ASA Z10. 2- ase), 


Co » coefficients; 


=a ratio defined byt measured flow rate to rate; 


= an average of maui ratios defined above; 


the coefficient in Manning’ s formula; 


Reynolds 


_ = the x-component of velocity (loca 


A, 


= the x- of velocity (maximum ex experimental— average 
= the : x- of of section), in feet sec: 


ie - ' the coordinate ax axes and distances measured parallel thereto; big 


¢ 


maximum distances, in feet; 


= the transverse 
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“ 4 aan _? AQUIFER TESTS ON PARTIALLY PENETRATING WELLS 
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to make use of complete penetration formulas with or without adjustments, de- 

used, are based on the steady- state solutions of partial p Such 

cedures may give fair results if the aquifer penetrated is of small and —— 
"thickness, provided that the assumptions made in each caseare obtained. nee 

ever, in the case of thick aquifers or aquifers of unknown thicknesses, ol 

for determining tl the hydraulic properties, a as well as the 


= _ Note.—Discussion open until February 1, 1962. To extend the closing date one month, 
7 a written request must be filed with the Executive Secretary, ASCE. This paper is part — 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 5, September, 1961. 
ia 1 Senior Hydrologist and Prof. of Hydrol., New Mexico Inst. of Mining and ie 
ani N. Mex.; on leave from the College of Engrg., Univ. of Baghdad, Baghdad, Iraq. 
» 1957. ad 
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The theory of nonsteady ilow toward artesian wells partially penetrating | 
finite nonleaky aquifers is used to outline methods for the determination of the 
formation coefficients as well as the thickness of the water-bearinoe forma-— 
— 
“2 
> 
— 
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September, 19 1961 
partially penetrate the are outlined herein. The are 
based on the theory developed by the writer3 for the nonsteady flow toward e 
‘Steadily discharging well that penetrates an artesian aquifer of infinite areal 
Notation. .—The _adopted for use in this" moe are defined 


artesian wells and that are ‘useful sing data from 
_ tests, will be presented. Details of other equations are available. 3 ae 
Drawdown in Piezometers. —Piezometers are small-diameter pipes driven 
- into an aquifer, so that entrance of water into the ~ is solely from the bot- | 


| ortesion —r>l. 


FIG, 1, _DIAGRAMMATIC REPRESENTATION OF PARTIALLY 
bi | 


PENETRATING WELLS IN ARTESIAN — 


“ing a distance * from the center of a pumped well that is screen between a = 
depths d and es and whose discharge Q is constant, is, for the case of a homo- 

geneous, isotropic, and nonleaky elastic aquifer, infinite in areal extent and of : 

uniform thickness b, as given below, and is valid only d during t the ‘period of — 
specified. (See Fig. 1 for the coordinate te system. Bi nce ok 


For a short period of ‘pumping, the equation of for 
-l- z)2 — 


E 


3 “Drawdown Around a Partially Well, by M M. Ss, 
ASCE, Vol, 87, » No. HY 4, 1961. 
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u, + 2) 


t is the time since pumping ¢ started, Ss ‘ond K are the specific : storage ie 

of water released from stora e ina “unit volume of the aquifer under a a 
head decline, of dimension L- 1) and the hydraulic conductivity of the aquifer, 

7 ‘respectively, and M (u, 8) denotes an infinite integral (see list of symbols), 
sufficient values of which have been tabulated’ for a practical range of the © 
parameters involved. A more detailed table of this function is available4 

Depending on the of the flow ‘system, Eq. may take the form 

a which case the equation is valid for t < c and being con- 

_ stants that depend on the parameters of ‘the flow system under consideration 


For example, if the pumped well is screened throughout its depth of penetra- Ba: 
tion (d =0), and if the is of a equal of ot the 


(, 


M 


If, on the other hand, the enometer is of zero penetration (z =0), the draw- 


case the equation of drawdown is 


For ‘relatively tong periods of p pumping, the equation of drawdown 

4 Professional Paper 102, Div., , New inst. of Mining: and 


ry 5 “Nonsteady Flow toa Well Partially Penetrating an Infinite Leaky Aquifer, ” by I M. be 
* Hantush, Proceedings, Iraqi Scientific Soc., 1957, p. 10; also reprinted by New Mexi- ; 


Inst. of Mining and Tech., Socorro, N. Mex. | 


| 
hd 
| = 


September, 1961 


wt 


K (2 


is known a as s the negative exponential integral of c u). The function is available ‘a 

- Drawdown in Observation Wells.—Unlike a piezometer that registers the = 
= at a point in the aquifer, an observation well screened in the aquifer © 


will reflect the average drawdown in the aquifer profile that is in contact with 
. the well screen (or perforated casing). For an observation well that is eeten: 
a between t. the he depths d’ d* and 1" the he required equations are a are as Seews: 


An approximate equation ofa average fort < 


that gives results sufficiently accurate for practical ; that 


om in which 5 is the average e drawdown in the cheervation well and = is the value 
- of the functionE of Eq. 1, ~ ame the value of z is replaced by (1' + d’) ee 


If > 1 and (1*/1) < 1, average drawdown in the o 
f can, for all practical purposes, be taken as that given by Eq. 1, with the value < 
of z arbitrarily chosen between l' and zero. The choice is generally made so 
= as to ) simplify the ‘equation which, in | certain cases, may take the form of — 2 


= U. S. Geol. Survey, "Werer_Suvpt Paper No. 887, 1942, p. 88; also ma lg 
by O. and E. F. Brater, John Sons, New York, N. Y., 1951. 


— 
is the zero-order modified Bessel function of the second kind, and W (u) is 
— 
1 
— 
— 
mm 6s» For example, if 1=3d, a choice of z =d will reduce Eq. 1 to Eq. 5c. Also, if © — 
im =0, a choice of z =0 will result in Eq. 5b, whereas a choice of z = l gives  - ii 
Eq. 5a. All these equations are of the type of Eq. 
For a pelatively long period of pumping, the equation of average drawdown 
i. 


TESTS 


a Eq. 11, as well as Eq. 6, shows that in this range of time, the rate of ae: 
4 drawdown is the same as though the pumped well completely penetrated the — 


won 
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2. —DIAGRAMMATIC REPRESENTATION OF PARTIALLY PENETRATING W ELLS 


the effect of penetration on on 1 on the drawdown has 
Drawdown in Piezometers or Wells for (r/o) 1.5.—For relatively large 
_ distances, that is (r/b) > 1.5, the equation of drawdown has been shown? to be : 


F In fact, Eq. y. 12 § gives s results. sufficiently a accurate for practical purposes a, é 
(r/b) as small as one, provided u < 0.1 (r/b) )2. The equation is the same 
as it would be if the pumped well completely penetrated the aquifer (Theis | _ 


s 


"by M. 
 : Hantush, Proceedings, Iraqi Scientific Soc., , 1957, p. 10; also reprinted by | New Me B20 
Inst. of Mining and Tech., Socorro, N. Mex. 
_ 8 “Groundwater Hydrology,” by David K. Todd, John Wiley and Sons, Inc., , New York, , « 
1959, p. 90; also “Arid Zone Recent Developments” | H. Schoeller, UNESCO, 
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Drawdown i in vedi Pumping Well.—The equations of drawdown that have been 
_ presented in the previous sections are derived on the assumption that the flux 
asi entering the pumped well is uniformly distributed along the water entry face of © 
the well, Theoretically speaking the hydraulic head, rather than the flux, along ay a 
a the face just outside of the pumped well screen should be uniform. "This ties 
head distribution can be achieved by a lengthy process involving the dis- 
tribution of varying flux elements along the well axis. The drawdown at points 
ae in the immediate vicinity of the pumped well, obtained by the simpler equa- 
tions that are based on a uniform flux across the well screen, does not deviate | _—. 
a appreciably from that obtained by the assumedly more exact, but lengthy pro- | _ 
cess, _ In the actual problem neither a uniform flux nora uniform head is real- a 
= conceived along the face of the well, “due to several involved field and — 
ational conditions. The drawdown, in this case, will have a value between the 
_ two theoretical extremes. _ Thus, “each of the two theoretical expressions can 
used to very closely represent the actual case, simpler equations 


de 


The theoretical drawdown equations, derived on the a assumption of uniform 
flux along the well screen, will obviously give a variable head distribution along 
the face of the pumped well. It has been shown, 9 however, that the maximum 
‘drawdown (least hydraulic head) at the face of the well that these equations 7” 
give, is very closely equal to that which obtains if the hydraulic head along the 
face of the well is maintained uniform . The ‘point along the face of the well at’ 
- which the least hydraulic head occurs depends on the space position of the well 
_ screen, as the well is screened throughout its depth of penetration, the least 
: bi head takes place at the top of the aquifer (z =0). If, on the other hand, the well 
; ae completely penetrates the aquifer and only its lower part is screened, the least _ 
_ head occurs at the bottom of the aquifer (z =b). But if the well is screened be-— 
_ tween the depths d and 1, the least head develops at a depth somewhere be- 
4 tween | and d, being closer to d if most of the well screen is in the upper half 
ct ¢ of the aquifer, and closer to 1 if the reverse is true. Computation in the draw- 


iffer appreciably from that obtained for z = (1/2) (1+ Consequently, if 
called well losses are excluded, the water level ina pumping well can be cal-_ 

_ culated from the equations of drawdown in piezometers (presented in the pre- 

7 ii vious sections) by substituting therein Tw y for r, and the value of z at which the © 


down equation shows, however, that the value of maximum drawdown does not 2. : 


least hydraulic head takes rw being the effective radius of the well. 
PROPERTIES OF DRAWDOWN EQUATIONS © 
Some characteristic features of the drawdown equations that are of impor- 
tance in graphical computations for the formation coefficients will be given. 
Equations for Special Cases. .—Equations of drawdown in piezometers 
deeryation wells that, during relatively short periods of pumping (t < [( (2b 
79 


on - rp)2  S5/ 20 K] ), are of the type of Eq. 4 or Eq. 9, are here called’ versus the 
special cases. Fig. 3 gives theoretical graphs of versus the 


oh de 9 “The Flow of Homogeneous Fluids Through Porous Media,” by M. . Muskat, } eee. 
‘Hill Book Co. New York, N. Y., 1937; or J. Edwards, Inc., Ann Arbor, Mich., 
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ort 


In the t < 20K is given 


‘is, small v values of time [the approximation of M has been 
_ Eq i indicates that in the initial period of pumping, the flow behavesas _ 
though the aquifer ended at the bottom of the pumping well. _ This is shown by 
Pak the dashed curves of Fig. 3 (I' and II'), that give the time-drawdown variation 
for wells I and Il, respectively, if the aquifer is assumed to end at the bottom 5 
of the ‘pumped well. Except for the pair of curves I and I', in the example of 
*Fig. 3 (prey ssumed values of the the other pairs (I and 


ges 


atiection 


—| 


mit 


drawdown. ‘The curves of each pair r deviate from each other, although 

the general trend of variation prior to the formation of the inflection point ap-— 

4 pears to be the ‘same, sie deviation | being greater for larger distances from the | 
an inflection point that develops within 
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TABLE 1. VALUES OF THE FUNCTION erf(x) & & £9 = ae erf(x) 


282(-7) 0. 34207 0.6448 

1134-6) 0. 35007 

254(-6) 0.35825 

451(-6) 0. 36656 

705(-6) 0.37497 0. 8528 
0.00102 0. 38357 0.8548 
0.00138 0. 0.39227 
0.07181 0.40114 
0.00204 . 0.41011 
0.00283 o. 0.41929 
0.00342 fo. x 0.42855 
0.00407 | 0.43801 
0.00478 0. 0.44763 
0.00558 0. 0.45738 
0.00637 0. 0.46727 
0.00819 0.48760 
0.00919 0.49797 
0.01026 0.50858 
0.01136 fo. 0.51931 
0.01253 Jo. 0.53023 
0.01384 0.54132 
0.01505 0. 0.55260 
0.01640 fo. 0.56402 
0.01782 0.57568 
0.01928 Jo. 0.58754 
0.02082 Fo. 0.59950 
0.02241 Fo. 0.61173 
0.02406 0. 0.62412 
0.02577 Jo. 0.63675 
0.02755 Jo. 0.64952 
0.02938 Fo. 0. 66256 
0.03128 Jo. 0.67574 
0.03325 Jo. 0.68921 
0.03527 Jo. 0. 70288 
0.03736 0. 0.71669 
0.03951 0.73079 
0.04i71 0.74514 
0.04401 Jo. 0.75965 
0.04636 0. 0.77446 
0.04877 Jo. 0.78949 
0.05125 Jo. 0.80477 
0.05380 fo. 0.82025 
0.05642 Fo. 0. 83606 
0.05909 Fo. 0.85203 
0.06183 Fo. 0. 86827 
0.06466 . 0.68480 
0.06755 J 0. 4 0. 90166 
0.07050 0.91869 
0.07354 Jo. 0.93612 
0.07664 Fo. 0.95372 
0.07981 Jo. 0.97161 
0.08305 Jo. 0. 98985 
0.08637 Jo. 1.0083 

0.08977 Fo. 1.0271 

0.09324 0.780 | 0. 1.0464 

0.09678 Fo. 1.0657 

0.10039 Jo. 1.0855 

0.10409 Jo. 1.1055 

0.10788 fo. 1.1259 

0.11173 Fo. 1.1466 

0.11566 1.1677 

0.11967 1.1091 

0.12376 1.2108 

0.12794 1.2330 

0.13221 Jo. 1, 2554 

0.13654 1.2783 

0.14098 1.3015 

0.14549 Fo. 1.3251 

0.15009 Jo. 1.3492 7.8827 
0.15477 1.3722 8.0271 
0.15956 Jo. 1.3970 9456 | 6.1751 
0.16441 Fo. 1.4236 0.9464 | 8.3252 
0.16938 1.4492 9473 | 8.4790 
0.17443 0.875/ 0. 1.4753 8.6349 
0.17956 0.880 | 0. 1.5018 

0.16479 Jo. 1.5286 

0.19014 Fo. 1.5561 

0.19556 Jo. 1.5839 

0.20109 4.6122 

0.20671 J 1.6410 

0.21243 Jo. 1.6702 

0.21826 0. 1.7000 

0.22419] 0. 1.7303 

0.23021 Jo. 1.7631 

0.23636 0. 1.7924 

0.24259 Jo. 1.6242 

0.248% 1, 8565 

0.25542 Jo. | 1.8894 

0.26201 Jo. | 1.9229 

0. 26868 0. 1.9569 

0.27546 0. 1.9915 

0.26237 Jo. . 2.0267 

0.28943 Jo. 2.0625 

oes | 0.29658 Jo. | 2.0990 

| 0.30385 Jo. 2.1360 

0.31124 2.17468 

0.31874 | 2.2121 

0.32638 0.995 | 0.8406 | 2.2509 

0.33416 | 1.000 | 0. 8427 | 2.2907 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
2. 
0. 
0. 


SL 


PS SSS SSIS 


— 0.0 1.525 15] 0.995) $89.78 0.990, 
0.100 1.575 65| 0.9965| 146. 0 1747.2 
1,580 70 | 0.9966 38 65 1796.0 
‘ 0.115 1.590 80 | 0.9967 139.2? 0.9997) 1697.9 a 
135 1.610 00 | [9970 2 2120.7 > 
170 $| 20 0. 9 
if 50 213.09 0.9998 | 28 4 
185 4.660 9976| 218.2 45| 0.999 
1.670 60 | 997 5 0.9998 | 2971 
i. 90 258.22 80 | 0.9999 6 
225 -700 | 2 85 3526 
a) 235 7981 | 277. 95 | 0. 99 
1.715 oS 98. 71 00 99 | 3844 4 
0.280 1983 10 4072 
i 20 | 298.78 0.9999 | « 
» (0.885 25 1983) 306.19 15| 0.9999 
4 1.740 0 | 444 
= LL? 4) 321.59 25| 0.9999 | 45 | 
0.270 1984) 330.12 30 | 0.9999 
1.790 | 337. 35| 0 on 
1.760 9985 45 ©. 9999 | 4991 
— 
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and t; is of t at the point. 


ble l for awide rangeofx,. | 

drawdown at the inflection point is given t by 


oa. ie 14 is obtained by —_ to zero the second derivative of Eq. 4 with res- 7s 
pect to logjot, and Eq. 16 is obtained from Eq. 4 by substituting uj for u. After 4 
obtaining the first derivative of of Eq. with ‘Tespect | to logit, uj is substituted 
The drawdown curves depart from those that obtain for an ‘infinitely thick 


aquifer (as described completely by Eq. 4) an 
given, respectively, 


> 


> a 


_§ r 2 
_inwhich x and are the of 14. Eq. 19 is obtained by noting 


"Regardless of the location ef the wells, the curvesare straight lines 
ing the same slope at relatively large values of >(b2 2K 
slope m of these lines is givenby 


a 
4 
4 
V4 
t 
a A Th inflection of the curves is clearly defined for small values of (r/b). a _ 
‘on the geomeetry of the flow system, the inflection disappears (see 
— 
| 
| 
+2 
— 
1 
&g 
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September, 1961 
which ‘is differentiating Eq. 6 or Eq. 10 with to and 
“noting tl that u becomes very small for large values of time. 


9), the properties previously examined still apply qualitatively. 


however, only Eq. 21 and, provided (r 8) is replaced by (1/2) (21+ 1" + d*). 


Eq. 18 are still applicable. In addition, as the length of the screen of anob- _ 

‘servation well is increased, other variables being held constant, the curve in- roe 
_ flection becomes less distinct regardless of the well location. _ The curve ap- ~ 


ae; proached is that of the Theis formula, which is obtained when (1"/) and (d'/b) 
become one and zero respectively (see curves Il and II’, Fig. 3. 


RECOVERY EQUATIONS 
the residual drawdown s' ina a piezometer during ey 


it: 
and end of pumping, 

‘similarly, the average residual drawdown in: well 


equation ‘corresponding to. any of the ‘drawdown equations dis- 
_ cussed above can be formulated readily, subject to the same timecriteria. For 


example, if the equation of drawdown is of the type that reduces to that of Eq. a 


short period of the recovery equation can writ- 
4) 
ae 
M 8) + B) - fy (ut, 8B (24) 
> 


we (u, - 
‘in which u' is the value of u after t by and 


u, ( + 8) (2 


17 can be 
If both to and re greater than s./2 K), 24, as welt as the general 


. 24 can be approximated 


— 

2 
a 
til 

— 

a 

| 


preceding give the recovery in “The recovery 
“4 in observation wells can be similarly obtained by using the appropriate 
age drawdown equation. — Eq. 28 gives the recovery in both cases. 


PR 


use is made of ‘some features of the equation. 
Recovery Rquations for Special Cases. .—Of importance is the equation 
water level recovery after shutoff in a ‘pumped well that is screened throughout 
its depth of penetration, or in nearby wells that may be as piezo- 


Residual drowdonw, 4° 


«FIG, 4, —LOG (t/t!)—RESIDUAL DRAWDOWN VARIATION IN WELLS PARTIALLY | 


on Equations of Drawdown), the during s short of 
_ pumping of the type of =u. 4or 9. Consequently, the equations given in the -_ 7 
preceding section apply. are here called “recovery equations for 


; “Fig. 4 represents theoretical plots of the residual drawdown st (for the set “| 
of parameters shown in the figure) versus the logarithm of (t/t), as as given by 
Eq. 24. These curves have three main characteristics: 
curves” have a straight portion at their lower left end; that is, 
q straight line passing through the point of zero residual drawdown and (t/t )=1, 
: _ This line may be clearly defined, as in the curves (1/ ug) = (4 K to/r2 Sg 106 
and 10*; barely recognizable, ,as in the curve 103, or entirely indistinguishable, 


as in the curves 10? and 10. The range of time i in a which this line forms de-— 


4 
— 

— 

— 
— 
1 A 4’ re in 
— 

tog 
— 


1 


pends on of the of pumping to, the formation ‘coefficients, and 
the thickness of the aquifer. That is, ee straight line forms in the r region for — a 
which (t/t') is less than |1 +(2Kt K to/b? 3) The slope of the straight line, 


which: is by differentiating Eq. 28 with respect to (t/t') 


serving that u and u' become very small for larger values of time. 
a 2. The curves have, in general, an inflection point, in the neighborhood of 


all practical ‘purposes is equal to the slope of the pe dng the inflection | 


@ 


The infl inflection point forms (if it does) in the period 


the r in the neighborhood of of inflection point is even 

by Eq. 24 after dropping out ut fo differentiating this equation (orany any 
twice with | respect to (t/ t) and id equating the results 


_ The slope of the curve mj at the inflection potnt can be obtained from Eq. 24° 
(or other equivalent form) by differentiating with respect to (t/t*) and 
ET uj and u'; for u and u', respectively, in the resulting expression. | . 
- ‘This slope can be approximated closely by either of the following two ee 
on the of (t/t' ) in point occurs: 


a if the inflection | point occurs it in the range (t/t!) < 100, the the shape 
mj = As'/cycle = 2.3c exp (- u "ye erf (x) 
If the inflection point occurs in the range (t/t") > 100, the second term 
of Eq. 31 becomes slope is. then elven 31 — 


a "holds, the curve under consideration will not exhibit an inflection point. ‘The “ 


an 
bea 

shown to be given approximately by the relation 

— 
4 — 
— 
j 4 
— 
4 
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ra _ itatively for other equations (see the section on Other Drawdown + 
also. Quantitatively, however, only Eq. 29 is rigorously valid. The other re- 
lations can be used only to give an estimate of the values therein. 
APPLICATIONS OF THEORY OF AQUIFER TESTS — 


the ‘effects of partial penetration resemble. the effects. of leakage from storage 
oe ina weap semipervious confining layer. 10 Also, if the curve inflection is ap- 
Lal parent, but the period of observation is not long enough to establish the ulti-— 
aa mate straight-line variation ona semilogarithmic time-drawdown plot, the ef- 
fects of partial penetrations resemble the effects of some kind of recharge 
boundary, such as induced infiltration from beds of streams or lakes, 11 or -« 
recharge from water - -bearing strata ‘supplying leakage12 through ‘semipervious 
a beds. The same general effects are observed if the wells completely 
penetrate a ‘sloping water-table aquifer or a wedge-shaped aquifer (nonuniform 4 
_ thickness). = The existence of such conditions should be reasonably eliminated 
prior to performing an analysis, based on the present theory. ae | 
Outline of Methods.—Semilogarithmic time-drawdown or logyg {t/t')—re- 
“ferve drawdown plots or both are an essential part of the analysis of the ob- = j 
oe data. AR the first Ssttiaile in n the analysis is the construction of ; 


under consideration, the methods are as follows: 


Inflection Point Method. —It the equation | of drawdown for the 
well under consideration during relatively short periods of pumping is of on 4 


2; of Eq. 4 or Eq. 9 (or can be 0 en ae by such equations), the follow - 

(1) On the semilogarithmic plot, construct the tangent in a the yn of me 


curve inflection and measure its slope mj = As/cycle. 


(3) Obtain the value of x and erf(x) from Table 1(or from a graph construct- 
(4) Compute u; from (x/8)2 and obtain the value of exp | -uj), — a slide 
“ Compute K, using Eq. 17 with appropriate expression for c. = Se 
(6) Obtain the value of M (uj, B), using tabular values of this function and the ; 
_ computed values of uj and 8. Then compute the value of sj from Eq. 16. cael _ 


ae “Modification of the Theory of Leaky Aquifers,” by M.S. Hantush, Journal of Geo- _ 


phy sical Research, Vol. 65, 1960, p. 3713. 
4@ 11 “Analysis of Data From Pumping Wells Near a River,” by M. 8. _Hantush, ‘Journal 


of Geophysical Research, Vol. 64, 1959, p. 1921. 

*Nonsteady Radial Flow in an Infinite Leaky Aquifer,” by M. S. Hantush and c. 
Jacob, Transactions, Amer. Geophysical Union, Vol. 36, 1955, p.95. = 
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September, 1961 
() From the semilog plot, the tj to to the computed 
- value of sj. Then calculate the specific storage from S, = (4 K tj uj ;/r? 
(8) Construct the ultimate straight line of the ‘curve (if it is 
_ in the observed data) and measure its slope m = As/ cycle. ~ Then compute 


o ) Compute b from the calculated values of K and Kb. 


Type- -C 


urve Method. the time-drawdown curve of the 


te observed data exhibits an inflection of the type shown in Fig. 3, and if the num- ay 
@ ber and distribution of the observed points are such that the details of the curve 
ee prior to the attainment of the ultimate straight-line variation are discernible, 

a type-curve method can be used to determine the formation coefficients and 
= often the thickness of the aquifer. The method is essentially the “Theis graph- 


bss offs tor 3) versus t on a logarithmic paper of ‘the same scale as ‘that of the 


type curve. The is as follows: 


ee. (1) Construct a ty type curve for each well, using 1g tabular values of the function 


 M Mi (u , 8) (ora graph prepared therefrom) and the E(or E) function of that well. 
— @ Plot the observed data on a logarithmic paper of the same scale as that _ 


_ Of the type curve. . 
6 Superimpose the observed data on the type curve, keeping the coordinate 


axes of the two curves parallel, and adjust until the best fitting position is ob- = - A 


7 tained. ‘ In matching the observed curve to the type curve, one should keep in 
‘mind that the observed points for relatively large values of time may deviate 
upward (having larger values) from | the type curve. _ The deviation is to be ex- 


pected, because the type curve is s for drawdown \ values slants g relatively short - 


ow Compute K from Eq. 1 or Eq. 8, ‘whichever r applies, and S, from th 


"imposed sheets and record coordinates; 
(or rE), 1/u) and [s (or 3), t 


(6) If the observed curve departs fr 


for piezometer, or 


e 


at the point of “departure. “mote the relations 


(1/2) [i + (1/2) (I +a") +r 


ae “Groundwater Hydrology, ” by David K. Todd, John | Wiley and Sons, Inc., New York, 


1959, p. 90; also “Arid Zone aes Recent Development, 
Le Paris, France, 1959, p. 37. om 1 


H. Schoeller, UNESCO, 


4 Ow eva under consideration ic or ction o 
p 
— 
= 
tm 
— 


= observation we wells or, in case E (or E) reduces to theM function of Eq. 4, from 
which ug is of uat the “departure” point 
(7) From the computed values of K and b, calculate the traseuntastyity of the 
a (8) If the observed curve does not ‘depart from the t type | curve ve within the ob- 
‘Served data, the value of (1/ u) of | a in the vicinity of the last ob- 


(9) It is to be observed that whether the data under analysis are pnd 
. 3 from piezometers or observation wells, and regardless of the location of the 
>) ot point of observation, step (8) of the ‘inflection-point 1 method for ‘computing the a 
_ transmissivity of the aquifer (Kb) always holds. Carry out this step and obtain 
Ps —* values of T ‘obtained from step (7) and ve (9) may not agree closely. — 
a This may be because the best fitting position of the two superimposed curves __ 
had not been obtained on the first trial, or because the ultimate straight por- P 
_ tion of the curve had been drawn either flatter or steeper than necessary, or 
= because the apparent transmissivity of the aquifer (in case of water- table — 
= ‘ aquifers) had not attained its uniform value within the period of observation, 


If such agreement is perme the procedure is repeated again with the neces- 


m is more likely that the apparent transmissivity has not attained a uni uniform val- of : 
For certain flow systems, both of the preceding methods may apply, in which we 
: case they serve to check each other. _ G Generally, the results of the computa-_ 
¥ tions do not agree on the first trial. The disagreement may bea resultof 7 
under - or over-estimating the slope at the curve inflection, as well as the slope 
of the ultimate straight line; it may also be due to a poor choice of the best | a 


matching position, or to all these factors together. A second trial with the nec _ ; 
essary adjustments generally | results in a satisfactory answer 
+S _ Theis' or Jacob's Method or Both.—If the semilogarithmic plot of the ob- 


or served data does not show any inflection despite a relatively long period of 


= pumping, three possibilities | may arise: (a) the observation well , although — 
~ close to the pumped well, is completely penetrating and screened throughout 7 
4 the aquifer. Such information is generally known a priori; (b) the observation _ 
well is relatively distant from the pumped well; that is, , (r/b) > 1.5. To aid % 
in deciding whether this is the case, information about the approximate thick- Gg 
ness of the aquifer should be available. ‘Such information can be obtained from ; _ 
previous geological or geophysical studies, or from analysis of data from a 
- nearby well through procedures such as those previously ehguytua orat least 
< (4/2) [ie + + (2. 4 V3 x) a (35 


q 

q 

— 
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x being the s solution of Eq. 14. 


_ If the first or the second possibility obtains, the method of | Theis and noo 
vided the ultimate straight line of the” ‘semilogarithmic plot forms, that of 
: me 14 can be used for calculating the formation coefficients. if the mt 
_ possibility pertains, the type-curve method can be used, provided the number 
and distribution of the observed points on the curved part of the semilogarith- 
mic plot oeielonty define a curve on the logarithmic plot needed for the type- 
curve method. Alternatively, the following procedure can be used, provided 
an estimated of the aquifer thickness is available: 
— Jacob's Method Adjusted for Partial Penetration. —During the time period 
‘tn which the ultimate semilogarithmic straight line forms, the drawdown is 


on the well observed, by either Eq. 6 or Eq. 10. Because the 


a fae (1) On the observed semilogarithmic plot, construct the ultimate straight 


line and extend it to the zero-drawdown axis. 
“7 _ (2) Obtain the slope, (m= As/ cycle) of this line and its time intercept, | 
cs) _ (4) Compute the value of f, (or fs) from their appropriate ex expressions i=) 

few terms of the series involved are generally sufficient), 
he (5) Compute exp (fs) or exp (f, s)|> using a slide rule or appropriate tables. 
Then calculate the storage coefficient 1 from 
Analysis of] Recovery Dats Data.- —If the equation of drawdown ! in ‘the observation 
. well is one that, ‘during short times, reduces to the type of Eq. 4 or Eq. 5, the ' 4 


) Construct the tangent to this” curve in region of the curve it mn 
and measure its slope m;=As'/cycle. 

(3) Locate the inflection point on this by inspection and its 
value of (t/t), which hereafter is designated by (t/t');. Over or underesti-— 
"mating the location of the inflection point will not materially alter mam calcula- — 

= tion to follow, as will be pointed out subsequently, 
(4) Compute the value of (to/ /t' i), which is = (t/t'); - 4 and ‘that of (t/to), 
(5) Compute f(x) from =f (x) and obtain the corresponding values of x 
: anderf(x)from Tablel. = 


“Hydrology,” by C. O. Wisler and E. F. Wiley and Sons, New York, 
1951, p. 234. 


| in which x is th econd possibilities are § im 
ruled out, the flow parameters are such thatthe relation 

| 

2 
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using 

(2) If the value of (t/t); > K neglecting the sec- 


(8) If the value of (t/t'); < calculate the values of ( (t'/t); andx y{(t'7t 
7 and obtain the values of exp /t), Jand ert rf (x t); i): "Then calculate 
(9) From the known values of t, and (t' i/to),, calculate th the 
Calculate the specific ‘storage from Sg = (4 K t'; u'j 
(11) If a straight-line variation between s' and logio ( (t/t") is discernible 


through the latter part of the data , including the origin point, t/t") = A, gts 0], 


Eq. 29; then, knowing the values of K and (K compute the value of b. 
It is to be observed | that if (t/t'); > 10, the second term of ‘Eq. 31 is com- 
paratively \ very small. . 4 Thus, a large e error in | estimating the value of | (t/ ™ =< 
will not affect significantly the value of this second term. The largest error, 
ce any will be in the first term of Eq. 31, owing to the factor (t; /t 9): This er- 
ror, however, rarely exceeds + 3%. Ev en if one misses the true location of the 
“inflection point by a factor of 10 (a remote possibility) and if, for example, the 
true value of (t/t"); is 10 and the estimated value is 100, the error in the val-_ 
ue of (t; /to), | and consequently in the value of K, will ‘be of the order (1.01 
= -1.11 1. 11 = -9%. The error in the value of Sg is approximately directly ov 
— to the error in estimating 
_ An exact graphical solution for the case in which the: short-time a 
, is not available. ‘The procedure outlined 
"previously ¢ can be followed to obtain a rough estimation of the formation coef- 
ficients. However, if the examined in step (11) forms, me solu- 


act, ‘because in the present case Eq. 29 obtains also. “ad 
_ "Examples of Application.— Aquifer test data obtained in the Rio Cenettn Val- 
{ Poi ley near the town of Socorro, N. Mex., will be used to illustrate the application ' 
sof the methods presented. . The we'ls tested penetrate a water- -table aquifer. j 
- However, because the observed drawdowns are small compared to the depth of | 3 
penetration of the pumped wells, the theory presented herein, although pre- es 
_ sented for artesian conditions, is still fairly applicable to water-table condi- 
bea tions. The computations in what follows were made with a slide rule. aes =. 
i ‘The Faulkner Farm Test.—A schematic representation of the wells tested 
is shown in Fig. 5. These wells are located in the SW 1/4 NE 1/4 sec. 35, T. 
,R. 1 W., , Socorro County, N. Mex. Both the pumped well and the observed 
- 4 well are perforated throughout their depths of penetration. These depths are 
116 ft and 60 ft, respectively. The observed well is 17 ft distant from the 7 
pumped well. An aquifer test using these two wells was run by J. F. Waldron!5 
during September 1952. During the test, the pumped well discharged at a uni- 


ee ‘New Mexico,” by J. F. Waldron, thesis goreceeteny to Stanford Univ. , at Palo Alto, Calif.. * 
1956, in partial of for the of ‘Philosophy. 
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< 2, the average drawdown s in the observed well during» 
relatively short = is given Eq. 9, which, with d' =0, = 60ft, 


Analysis of Drawdown Data. _The observed ‘semilogarithmic c curve of Fig. 5 


= clearly shows an inflection, and because Eq. 24 is not of the form of 1 5, the 


au 


m7 7.4 ft/cycle 


— 


avn 
‘sh 


FIG, 5. ~—TIME- DRAWDOWN VARIATION IN TE THE. FAULKNER TEST 


| Type- e-Curve Method.—From Table 1 and E (u) of Eq. 38, . 3 type curve is a 
prepared on logarithmic paper, and the observed data are plotted on logarith- 
_ mic paper of the same scale. The two curves are superimposed and adjusted 
3 to obtain the matching position shown in Fig. 6. The valeus E , (1/u), 5, andt 


at the matching point shown in the figure are 5. 6, 500, 10 ft, and 100 min (6, —_ 
. sec), respectively. Eq. 38, is computed ed as 6.8 86 x 10-* and from 


iy, 


s ‘elon the period of the test, the observed curve does not depart from the - a 
t 


is, the time ta at the eventual ‘point of departure is greater 


4K 1/(2) 22 , Sg is computed as 1.14 x 10-4 


— 
= 
1 


AQUIFER 1 TESTS 
than the period of observation. . Thus, fror from Fig. 6, tg and are greater 
than 500 min and 2, 500, , respectively. Consequently, the thi thickness ¢ of t the aquif- 


(1/2) fl + (1/2) I" S/u = (1/2) [116 + 30 - +17 x 2,500] = 1023 
— Inflection-Point Method.—It has been pointed out that the application of the - 


_ inflection- point method on data obtained from flow systems for which Eq. 4 or “a a 
9 | does not obtain will generally give an approximate solution. ‘This, of 
course, is due to replacing the actual equation by an approximate equation of — ; — 
: the form of of Eq. 4. ‘4. For r instance, | if the drawdown in the observation wel well of this 


Fm. 6. ,—OBSERVED DATA AND THE THEORETICAL TYPE CURVE FOR | 


_ THE OBSERVATION WELL OF THE FAULKNER TEST wie 


% shown in Fig, 5. "The tangent to the curve in the region of the ‘er inflection -- 
is constructed; its slope m; is measured as 7.4 ft per cycle. 5 =% 


= (6.82)2//m = 26.2. From Table 1, x = 1.666, and erf(x) = 0.982. (x/B)? 20 
=0. 06, and exp(-0.06) =0. 94, From 17, K is computed as 7x 10- -4 fp fps. h- 
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a September, 1961 5 

- terpolation in tables of the function M (u, / 8) ¢ gives M (uj, | 6.8 2) = 2. ne = . 
= Eq. 39, s; is computed ; as 8 ft. For sj =8 ft, Fig. 6 shows t; = 3.1 min = 186 sec. 

f ‘Then, from Sg =4 K tj uj/r2, Sg is computed as 1.08x10-4 ft. 

The ultimate straight-line portion of the curve is not discernible within the 
= data (see dashed line of Fig. 5). Hence, the transmissivity (Kb) and, | 
consequently, the ——— of the aquifer cannot be determined from the pres- 

Olson Wells” Test. lage 7 (a) shows» semilogarithmic | time-drawdown 

curves for two observation wells prepared from data obtained from a pumping = 
test on what is known as the Olson Wells. The mcoryes.deg of these wells is shown Bt 


tL 6 i, The other well, ‘designated a as s observation well (1) is s 182 ft east of the 
_ pumped well. A third well, designated as observation well (2) (known as the 
Lopez dug well), is 384 ft west of the pumped well and is assumed to be zero 

penetration (actually about 1.1 ft penetration). 
_ Analysis of drawdown data from well (1)— Eq. cs will, of course, give a 
"average drawdown in this well. However, because (r/)) = (182/75) and 

= = (75/75) <1, 9 can be approximated by Eq. 1, in which z is chosen 
arbitrarily between 0 and 1 (see section on average drawdown for short times). 
_ For z = 25 ft (this choice is made because Eq. 1 will reduce to the form of _ 


4), the drawdown in this well for me s will, from Eq. 


= 
4, with 


and B= =0, 55. Because the in- 


of the semilogarithmic curve the data is clearly manifested 
(see Fig. 7 (a)), the inflection-point method and/or the type- curve a will 
provide a sotution for the formation parameters. 

(u, 0. 55), prepared « on logarithmic paper, ,and the observed 
data are plotted on logarithmic paper of the same scale. The two curves are 
superimposed and adjusted to obtain the matching position shown in Fig. * 
The values E [M (u, , 0.55) )], (1/u), 5, and t at the matching point shown are 0. se, 
10, 1 ft, and 1,000 min (6 x 10-4 wel —- From Eq. 40, K is calcu- 


3.5 x ft- The value of f (1/uq)at t the point t of departure is Then, from 


the relation b= (1/ 2) [re r b is calculated as 695 ft. “The transmis- 
_ sivity of the aquifer is T = Kb =0. 319 sq ft per sec, and the storage ——— 
Inflection -Point Method.—The observed ‘semilogarithmic plot is prepared 
as shown in Fig. 7 (a). The tangent to this curve in the region of the curve in- 
flegtion is constructed. _ Its slope m; is measured as 1.55 ft per cycle. f(x) 
= (0.55) =0. 17, From ‘Table’ 1, x =0.371, and erf(x) =0. 40. uj 
=(x/8)2 = 0.454, and exp (-0.454) = 0.635. From Eq. 17, K is computed as" 

78 x 10- -4 From tables of the function (u, ‘B) [or of 


| lated as 4.84 x 10- 4 fps, and from Ss = 4K t/(1/u) r2, S, is calculated as. 


| 
a 
ge = campus of the New Mexico Institute of Mining and Technology, Socorro, N. Mex. _ 
— 3 _ Each of the two Olson wells has a 75-ft depth of penetration, of which the bot- | 
— 
— 
4 
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_m=090- 
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Residual drawdown 4, ft 


~TIME-DRAWDOWN VARIATION IN THE OBSERVATION WELLS OF THE 


OLSON TEST; b) LOG (t/t') - RESIDUAL DRAWDOWN VARIATION INWELL 


 Elu) = M (u, 0.55) 
forwell@) 


FIG, —OBSERVED DATA AND THE THEORETICAL 


: FOR WELL (1) OF THE OLSON TEST 
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0. 55)], M 0. 58) 0.33, From E Eq. 40, 5; =0.88 ft. From Fig. 7 (a), 
ty = 220 min or 13,200 sec. Hence, S, = 4 K tj uj/r2 = 3.46 x 10-4 ft- The 
“is ‘geese straight line of the curve in Fig. 7 (a) is constructed. — Its measured 
slope is 0.90 ft per cycle. From Eq. Kb is ft per sec. 


Kb\ 
Consequently, b= =(. 326/ 4.78 x 10- 4). 680 ft. The storage coefficient 


The calculations are those of a trail. In the first trial, 
hydraulic conductivity as obtained by using the inflection- point method was 
a underestimated. Apparently, this resulted from drawing “the two tangents” of _ 
wa the observed curve steeper than was necessary. _ ‘The coefficient obtained by 
peat the type-curve method was overestimated _ because the best matching 

— was not obtained, it appears, on the first trial. . The two tangents were 

redrawn with flatter slopes in the first case, and another fitting position was 

obtained in the second case, , because the observed data permitted these adjust- 

| ha The eaheetaene were repeated, with the results being those previous- 


discharge within probably the first 100 n min of the pi pumping period (because of. - 
Bt ess characteristics), the observed points within that period were rejected. — 
_ In addition, one should note that during the very early period of pumping, the 
observ ed data generally do not follow the theoretical drawdown equation, be- ,? 
cause of the probable variation of the formation coefficients, especially if the 


; _ Analysis of L Drawdown Data from Observation Well (2). - — The observed sem 


Tacob’s. method does are ruled out, although the general tendo of the 
time. -drawdown variation resembles that given by the Theis formula. 
_ Jacob's Method Adjusted for Partial Penetration.—Although there are enough 


fair solution by using the type-curve method (drawdown equation is Eq. 1, with 


a points on the curved part of the observed semilogarithmic plot to warrant a 
z= =0), , Tacob’s method adjusted for pi penetration is followed for = 


and tp =270 min (16,200 sec) Fig. (a). 

calculated as 0.326 sq ft per sec. From Eq. 6, f, is a te pep a value 
—of.b = 700 ft (only five terms of the infinite series are needed), as 0.81. a 
exp (0.81) = 2.25. From = = [2.25 25 T ty exp (fs)/r? the storage coefficient is 

Analysis of Recovery. Data from Well (1).—Fig. 7 (b) is is the observed curve of 
residual drawdown s' versus the logarithm of (t/t'). The tangent in the region 
| of the curve inflection ‘is constructed; its “measured slope mj; = As' per cycle 

= 1.65 ft per cycle. By inspection, the inflection point is located at (t/t"); = 30. 


Because (t/t 100, all terms of Eq. 31 are computation is 


= 
* 
— 
4 
curve inflection. The well is known to be of practically zero penetration, and 
Se &§ the aquifer is known to be deep, its thickness, as estimated from analysis of 7 
ee : _ data from well (1), is about 700 (actually an average of 670) ft; consequently, Z — 
t 
7 


(tite) = (t/t)4/(to/ t's) = (30/29) = 
BNF = (0. = (0.55 “0. 


nd 
= 0. 454) = = 0, 635; 


- 


we 


The straight portion of the lower left nities of the observed curve is not 
clearly discernible. | ‘However, for the purpose of illustrating the procedure, — 
the straight line is as shown by the dashed line in Fig. 7 (b). 
measured slope m is 0.90 ft per cycle. Then, from Eq. 29 K b is computed as 
Finally, b=(K b)/K = “(0 326/4. 63 x 10° 4) = 705 ft, and 


The following symbols, adopted for use in this paper, conform essentially 
with “American Standard Letter Symbols for Hydraulics” (ASA. Z10.2-1942), 
- prepared by a committee of the American Standards -_—a with Society 


= thickness of aquifer; dimension | 
: a = constant, depending on the well discharge, length of well screen, and > 
the the hydraulic conductivity of the dimension 


depth from the top of the aquifer ¢ of the unscreened ‘portion (unpe 
forated secti f the casi of the ed well, ‘iL; 
forated section of the casing) of the pumped well, 
_ depth from the top of the aquifer “ ‘the seinen portion of the ob- 


f (x) The error functio 


— 
» 
: (See 4 alues OI p al 
0.40. Then u'; = (x/8)2 = 0.454; exp(-u';) 
exp (-0.454/ 30) = 0.985; and eri(x/ 
 &g From Eq. 31, K = [(2.3) (1.6)/(400m) (1.65)] [(1.035) (0.635) (0.4) - (0.985) (0.077) 
(1.77 x 10-3) (0.264 0.003) = 4.63 x 10-4 fps. From (to/t'))=29 and 
gg 
| 
— 
— 


- 


a/b, 2/b) ="2b/7 (1 - [sin 1/b) 
cos (na z/b) W(u, na r/b); 


‘ 


with 2 Ko (nt r/b) replacing W (u, n7 
= hydraulic conductivity of the aquifer, 


= the zero- order modified Bessel f function of the second kind; 


ervation hole, ag 


= = radial ines ‘measured from center - of well, 


Be 

= drawdown of piezometric surface at any ‘time and at any _— in the 
average drawdown in observation holes, L; __ 

‘residual drawdown in piezometers, 


L; 


ig 
= transmissivity of the aquifer, L? T 
=time since nates, 


— 
meric . 
iE 
— 
— = 
— 
— 


exp (-y /4 fant 


= vertical | coordinate measured from the top of the 


downward and 


: ao as = one of the parameters of the function M (u, B), also it is a a constant 
depends on one or all of the variables, I',d, a » 2, 


of 
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‘studies 
five steep and creeks in 1 this area, the Florida, LaPlata, 
: . and Animas Rivers, and Junction Creek, a tributary o of the "Animas, ea 
bo nine of the ten locations there were sections at which critical velocity 
governed. Numerous other sections required checking to determine if critical 
conditions were present. As dams are buildfarther into the mountains, criti- 
cal velocity conditions will become more common. 
In Table 1, p. 84, “Hydraulic Characteristics, ” in order to cover critical — 
a mentioned on p. (82, the headings, the ‘Durango office has found helpful 
are: Elev., A, T, 7 top width anda starting or check point for wetted peri- i 
‘meter, p) n, Kg, A/T, V cr (which equals 5. 
elsewhere? , and curve versus 5.67 (A/T)1/2, 
Even if aad or ‘Ge alone is used in computations, depending on ‘the ieee 
; to this difficult subject, the use of either or both precanyy should be permitted. — 
If there is a strong proof that n varies with he. ann— stage “ve! should be 
drawn and Kg values computed from it. 
_ In November 1960, the Durango office was s advised that Table e 1 had been 
_ revised to give the following headings: Elev., a, p, r, r,2/3 Ky = 1.486 ar, 2/3 
an, and Kq = Kp/n. These headings were recommended for use in cases inwhich 
me Table 2, p. 91, “Water Surtace ‘Profile Computations— Method A,” the 
Durango ‘Office uses column for which should be logically placed be- 
fore Column 6, which equals (Q/ Ka) - . This new - column may be called 5A, 
is used also in Computing Column 9, ‘steps 2 and 3, Pp. 86), when there are 
ca It would be logical and helpful to put ; as many ‘step directions and fixed rules 2s 
_ of procedure as possible, in fine print if necessary, at the column headings, so. 
a that: the inexperienced | computer would not have to refer to some other source 
When 


2. He might make mistakes here column 18 would provide him a place to 
Veet 4 April 1959, by J. Lara and K. Schroeder (Proc. agi 1997). 7 
1 U.S. Bur. of Reclam., Durango, Colo. 


Colo., November, 1957. 


— 
— 
EDWARD N. WHITNEY,! F. ASCE.—Methods A and B and four other 
2.4 
4 
Investigations, Bur. of Reclam., U. S. Dept. of Interior, 


record the. « , plus or minus, and awide Column 19 would space 


of adverse slopes which occasionally occur; 
3. Possibility of change in n due to sediment deposition or timber being 


logged or burned off; 


_ 4, Accuracy of estimate; and a 
Int within 0.1 ft water surface eleva- 
ot ‘tion, Column 2, and water surface elevation, Column 17, it has proved a j 
ry ; to plot a curve, called a Try-Get Curve, the Try elevation on the vertical scale, : - 


= and the Get elevation on the horizontal. i alt 45° line is drawn between points of > 


= letter from the Denver r office in n November, 1960 gave the information that 
Lara used this ; method, when the computations were begun. 7 His computation 
was also shown for a straight line intersection between two points located on 


 ‘glaanie that will locate a curve, whose Distal with the 45° line will give = 
the desired answer. The curve solves the varying relationships between the y 


dimensional and characteristics of the two sections in When 


erence may well ‘depend on the _ of the discharge, RY 


and the 
i. aa The shape and slope of the curve developed for one Q can be used as an aid a 
aE in estimating water surfaces for lower or higher Q’s. ‘Critical velocity cond- 
| are shown by failure of the Try-Get curve to cross the 45° line. For the 
“a last upstream section only, if pressed for time, and if the estimator is sure of — 

his points and | curve, he can pick the correct elevation without further om 


understood that some computations for all or parts of Method A have 
heen done on electronic computers. Has it been possible to eliminate the cut 
* and try method now in nto in this gay would undoubtedly be ape 


? 4 
the Bureau of Reclamation 


M. LARA,3 A, M. ASCE, AND. KENNETH 
to Closure) The report,2 referred to by was pre- 
pared in more detail than the paper by the writers. 
= Critical conditions that can occur in mountainous streams do require special 
analysis. Whitney presents a helpful suggestion in | tabulating: the various hy- 


draulie characteristics 

__versus 5.67 (A/T)! 
_ believed he meant the plotting of Ver versus A/T. _ The form on which the hy- 
or modified to suit the type of 


analysis being made. 


referred to by Whitney is apparently in error. It 


Hydr. Engr., Bur. of Reclam., Denver, Colo. 
Hydr. , Bur. of Reclam., Denver, Colo. 
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_ DISCUSSION 
Whitney” ‘comments on n suggested a additional columns 1 in the water surface’ 
* prot computations form will be considered when the next revision to the — 
The Try- Get Method of reducing the number of trial assumed water surface _ 
: _ elevations is an informal standard procedure often used in these computations. — 
a "9 It has been found, however, that the curve traced by a plot of the Try-Get ele- 
a vations is generally shaped slightly concave downward. A plot of this curve 7 
aa does serve to indicate the occurrence of critical conditions by | not learned 
ing the 45° line as mentioned by 
oer The Denver Office of the Bureau of Reclamation is currently checking a> 


program of computing water surface profiles by Method A using the IBM 650 


che 
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. Because of the experience of making these measurements, the writer felt that 


there is opportunity for correction and expansion. star 
‘Thus, the discussions presented are a valuable addition to the material herein. 
ex Bush and Ho have emphasized one of the key problems « of settling analysis. 
Boreas can ‘adequate initial mixing be provided before settling starts without leav- 
ing enough turbulence to affect the first few minutes of settling? During | the 
xtc described herein, ‘considerable time and effort were devoted to this” 
a question. In one series of tests a suspension was made up of waterand a small — 
a, of the standard asbestos used for preparing gooch crucibles. This sus- . 
_ pensionwas mixed by the same method used for the clayand alum suspensions. 
At the end of the n mixing, the asbestos 8 suspension \ was poured into th the ee 
by means of a large funnel. The n motion of the asbestos particles v was ob- 
he served. Within about 45 sec only a few particles were observed to have up- A os 
ward motion. By the end of 120 sec, any upward | motion of particles was ‘slow 
similar to motion caused by ‘convection currents. Even this motion dis- 
‘appeared in min or2 min. Thus, the turbulence did not seem to be significant 
_ for 10 min as reported in the study. 
3 et During the test with asbestos particles | an attempt 1 was made to | guess the 
~ order of magnitude of the coefficient of diffusion, ez. A number of cane 


analyses had already been “made ‘in at turbulence jar of the type described by 
Rouse. 15 As part of these analyses, several values of had been measured. 


in the tests with asbestos particles the value « of ez , would be less than 10 — * 


tube and less than 5 cm2 per sec 2 min later. 


7 per: sec immediately after the end of pouring the suspension into the ae . 
‘The magnitude of e z is not the only factor to be considered in meet & - 


in the analysis shows that the slope of the concentration profile must be known. a 
4 From Figs. 7 and 9 of the paper, it will be seen that this slope is small for - 
small values of except near the surface of the suspension. Similar results» 
7 5 were observed during settling analyses of domestic sewage. " The e reason for 
a _ these small slopes has been presented in Section 4 (e) of the paper. =~ 
om fe The slope of the concentration profiles in Fig. 7 can be measured directly. ee 
At t = 250 sec ‘this slope is. negligible for depths greater than 40 cm. 
2. gm per cm4 


a December 1959, by. Ronald T. McLaughlin (Proc. . Paper 2311 
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SETTLING PROPERTIES OF SUSPENSIONS? 
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$= 0.0120 ero 
250 sec 
c 


36 | ©. 012 0.013 cm per sec 


Therefore, the estimated at depth of 20. cm a time q 
sec is equivalent to a mean settling velocity of about 0. cm per sec. 

For smaller values of t, the value of e, would be greater, but the slope of 


. the concentration profile would be smaller. As a rough estimate, it might ~- 
that the first 250 sec of the analysis the mean vas 


greater degths the effect would be less because the slope of the is less. 
For example, at = 80 cm and t = 500 sec, , the slope of the profile is about 
0014 @ (z, 0) gm per cm4 | whereas e, might be only 2 cm2 per sec. Since 
= 0.925 9 (z, 0), the equivalent settling “velocity would be only 0.003 cm per sec. 4 : 


it is seen that the increase because of decaying turbulence { is small pint is 

Bier, writer feels that whereas the analysis is correct in principle, the effect | 

f the turbulence invalidate the conclusions for the particular set of 


Re ai: Hall, in his study, is concerned with | using settling analysis as part ofa 
a practical design procedure, Therefore, he has developed a simple test that - 
_ feels is more practical. It is not, however, correct to assume that a test is 
& more practical because it is more simple. Pa Whether a test is practical or not 
a is determined by the useful information or understanding supplied for a given 
amount of time, money, and effort. Ifa test is so simple that it does not bring 
to light some important information, it may be quite impractical. The writer 
= spent as much as a whole day talking with engineers who are trying to 
solve a particular sedimentation problem. — They are usually certain that any 
a, a but the simplest | test will be impractical. In the same amount of time, one 


aii 


ax it has often been possible to obtain a sample of the suspension . bring it to the 


laboratory, and perform a settling analysis. Once a simple multiple depth tube 


has been constructed and calibrated, the analysis can be performed in any lab- 
_ oratory that is set up for determination of suspended solids. 
The test developed by Hall is based on three assumptions about how the 
. thing of sewage proceeds. _ Steps 1 and 2 and the related assumptions, a and b, 
of the study can be accepted. However, step 3 and the resulting assumption, 3, 
do not seem reasonable. Step 3 states that floc have negligible velocity until » 
: they reach a critical size and then they settle rapidly. -Floc, once formed, —- 
simply particles; they may have complicated structure and shape, but they are 
4 particles nevertheless. If they behaved as assumed by Hall,a graph of settling 
velocity as ordinate against particle size as absiissa would show no velocity — 
until a critical size were reached, At that : size, the velocity would suddenly ; 
jump toa large value. als a graph is not to be found among the curves that 


es = 5 cm2 per sec gives (2) ] 
Mince =0.96 oz, 
— 
q 
— 
me 
— 
— 


; bg fs _ the same velocity until they reach the bottom of the suspension. When these 


tion. Ifa particle i increases in size only lo teuinn' with particles | small enough 
7 to have no velocity, the increase in size must be gradual, Hence, the increase 

<> ’ 

. invelocity must be gradual. There could be a sudden increase in size if a floc 

a ul particle joined another particle of significant size. But then both of the parti- 

i cles would be of significant size before joining and would both have a settling © Z 
velocity, which is contrary to the assumption. The same argument holds for 
particle shape and orientation. th 

Hall also states that samples from the surface of | quiescent sewage reflect 5 . 
‘a settling characteristics of the sewage. This statement can be examined by © 
ie considering the rectangular settling tank shown in Fig. 11 (a ) of | of the paper. - 
Taking samples at intervals of time from the surface of the suspension ina a 
‘Settling tube is equivalent to taking samples at the surface, z = 0, of the i. 


i » successive values of x. The removal of particles due to settling is de- 


termined by the complete. concentration profile at each value of x. Therefore, 
"measurements at the surface will be related to removal only if the concentra- 
a) - tion of particles at the surface is the average of the concentration throughout 4 j 
the depth of the tank. - Unfortunately, the surface is the one place where the 
ae is least likely to equal the average concentration; it will almost _ q 
_always be less than the average. If,as suggested by Hall, the efficiency of set- . 
; tling is based on measurements at the surface, it will usually be impossible to _ 
obtain 100% efficiency, even in an ideal tank. 
- inp Fischerstroem questions the way in which the linesare drawnon the various 
concentration profile diagrams. He suggests, correctly, — the lines of con- Bs 
ee stant t should all pass through the origin, 9 (z, t) =0, Z= 0. In Fig. 5, where 
the diagram is explained, the lines of vd saan t are drawn in this manner. 
_ are not in Figs. Tand 9 because there were not enough data nearthe ‘sur- ‘aa 
_ face to give the exact shape « of the profile. | For purposes of computation, each a ~ 
_ profile was assumed to have a er line section from the origin to the up- i 
It is also suggested in the analysis that the lines of constant z/t should all Y 
pass through the point 9 = = 100%, z =0. ¥ In order to settle this question, it is 
necessary first to consider a suspension of particles that settle without floc- a 
- culation. . The particle concentration in this suspension is assumed to be small ae 


enough not to affect the settling velocities of the particles. | In such a suspen- 
7 _ the particles of settling velocity, wj, at the beginning of settling will have + 


‘ oi ik -particles settle in a tube of the type shown in Fig. 5, their « concentration will ar 


a function of only z and t. Hence, it can be written (2 (z, t);at t=O this func- 
Because the particles” settle without flocculation or hindrance, the spatial 
‘rte of i- particles settles as a unit, independent of the particles of on an 
ae velocities. In time t, the whole unit settles a distance wjt. It follows that a ae 
the concentration of i- particles at the depth z and time t is the same as 4 & 
concentration was at the depth z - wit and time t =0. : it is possible © 


one usually finds in the literature. Therefore, the sudden jump in velocity 4 
| 
4 
| 
: 
at 
Die: 
os = 
— 


ped Both sides of Eq. 4 can be summed over all values ra i to vege : 


greater than z in time t, and it will no longer be in suspension at that point - 
Therefore, the summation need only include values of i up to wj = z/t. Eq 5 


7 ae han 


Eq. 7 is the equation of a line of constant z/t ina concentration profile di- 


_ It is now possible to examine what happens when both z and t approach 


zero ) along one of these lines. . In order to avoid mathematical complications at ae 
the limit, itis ee to substitute z/K for t in Eq. 7 and let z approach 


| 


As z goes to zero, Bq. becomes 
e of Eq. is ‘simply the pavtictes with 
_ settling velocity equal to or less than K at z = 0 and t = 0. It follows that the 
points where the lines of constant z/t ¢ cross the axis Z = =0 on the a a 
profile diagram. are determined by the frequency function of the settling velo- 
: cities at the surface of the suspension at the beginning of settling. In fact, the 
‘Spacing of the lines of constant a along the axis: will | give directly the values 
7 
=0, t=0. the 
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7 
— 
— 
— 
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— 
—— 


ension at t = 0, this cumulative distribution at the surface is also the cumula- an 
‘This result, obtained for settling without hindrance of flocculation, also” 


shows what happens when these two processes occur. At the surface of the 


_ suspension at the beginning of settling, there has been no time for flocculation | a _ 


to take place. The settling velocity d distribution is determined by the nature of 
the particles, the properties of the fluid, , and “the concentration of particles. 
_ For that initial instant, before change in concentration or flocculation have be- 
a _ come significant, the settling is similar to that for which Eq. 9 was derived. 
initial, instantaneous frequency distribution of settling velocities 
_ mines where the lines of constant z/t meet the axis. - Immediately after this — 
"instant, flocculation and change in particle concentration have produced some 
effect. Thereafter, the spacing of the lines of constant are no longer 
lated to any frequency distribution of settling velocities. nde, 3 i 
“ a Because it is not yet possible to measure the ) frequency distribution at z =0 
F and t = 0, the lines of constant z/t i in Figs. 5, | * and 9 have not been extended 
aa to the axis. ;. Otherwise, if the results of a multiple - depth : analysis indicated 
iid that the lines of constant z/t could be extrapolated to the axis with te eg = 
the result would be the initial instantaneous frequency distribution at z=0. If, | 
Re in addition there were good reason to believe that the particles of eachsettling — 


velocity were uniformly distributed at t = 0,the frequency distribution obtained 
would represent the whole suspension at t=0. This is only one of several 
7 ways in which the lines of ci constant z/t can be \ used to analyze the results of = 
multiple-depth settling analysis. Others are explained in the thesis 
ic. that describes the research on which the paper is based. eae 


indicating settling tanks could be designed solely on the results of set- 
tling analysis. The intention was to show how the settling properties of the 
oe suspension ‘should be considered as part of the design study. In this regard, it : 
_ was stated that for certain suspensions, the removal of particles during set- | 
-tling will depend more on detention time than on overflow rate. For other 
Pensions, the removal will depend more on overflow rate. _ The mera ; 
oe profile diagram will give a good indication of which one is more er 
The writer also concluded from Fig. 9 of ‘the paper, that at depths greater 


time. _ It is to be emphasized that this conclusion pertains only to the suspen-— . 
sion: used by Fitch and is not general. Fig. 1 of the discussion actually sup- 
_ ports the conclusion, — The lines of removal at constant time are approaching ~ 
.; vertical as the depth approaches 5 ft. There is good reason to believe that 
at greater depths these lines will be essentially vertical and ‘Straight. For 


_ The analysis attempts to establish the superiority of shallow shai tanks * 
me This point will be considered in order to show the error of making a general 
ze about settling tanks without considering the suspension. Fischer. - 
stroem states that for a given detention time, a shallow basin will always give | 

a greater removal than a deeper one. _ This is true oer the particles | have 

Otherwise, it can also” 

be said that when. flocculation occurs particles a are e continually i increasing their 


~ velocity as their size increases until a } Seniting size is reached. Therefore, an 


settling tanks. It should be that there was no intention in the 4 
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, settling velocities are distributed 
— 
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— average cneaten velocity in shallow basins is less than the average in deep- 
- erbasins. In decreasing the depth, only the slowest part of the particle | motion 
Fischerstroem makes a similar statement in another way. The only possi- 
a im he states, of increasing volume load is decreasing the overflow rate. 


“This cost ine 


struction, as well as ‘operating costs. “The cost of excavation depends on tank *. 
- volume unless there isa change from | soil to rock within the depth of the tank. 
eS The cost of land depends on the horizontal area of the tank, while the cost of 
- construction probably depends on the area of the inside surface of the tank. 
For equipment such as scrapers to remove sludge, the ) cost probably i — 
with area more than with depth. ~ on! gt 
To show how these items are related to the settling properties, the 
a — ‘used by Fischerstroem will be used: Assume that 90% of the suspended solids © 
\BLE 4,—COMP RIS ON OF TANKS OF 3 DEPTHS 
; 
ft 


7 


0 1 


settling tank by Fig. 11 (a). of the 
cussion, the overflow rate, depth, and detention time can be determined 


_ depths. ‘These are shown in 1 the f first four columns of Table 4; two of the pf 


_ flow velocity, U. These are shown in Table D1 . The greatest velouity is 0. 035 
‘ per sec and it occurs in the tank of least depth. - Work by A. C. Ingersoll, F. s 
ASCE, and McLaughlin36 indicates that if the mean velocity in the shallow tank No 
_ were much higher, the settled particles at the bottom would be resuspended by — 

wt: the flowing ‘suspension. Hence, the assumed value of flow is an upper limit for 


a a The volume, L B D, the horizontal area, L B, and the area, S, of the inside 


Surface | can now be computed, and they are ‘shown in Table 4. For the purposes | 
a s comparison, the ratio of volume to the volume of tank number 3 is shown in 


-36 “The Resuspension of Flocculent Solids in Sedimentation Basins,” by A.C. Inger- 
or soll and R. T. McLaughlin, enn < Inst. of Technology, Sedimentation Lab., May 1960 0 


= 


— 
— 
< 
i. 
zpimarrry objective in the design of settling tanks is to use the smallestvol- 
ume for a given flow of suspension and percentage removal. The real 
ee s _ tive, however, would seem to be the smallest cost for a given flow of suspen-_ ; - ; 
— 
145 1,00 1,00 1,00 
OVErTIOW Fates are those chosen by rischerstroem, and the third was Se- 
‘the purposes of co tation assume that the 1 h-width = 
tanks will be 20:1 and that the rate of flow will be 1,250 ft’ per hr. Thereis  =§ 
- 
q 
| 
= 


- decreasing the tank —y from 5.0 ft to 0.8 ft the volume is reduced by | about 
30%. At the same time, the horizontal area is increased by: about 330%, and the 
7 inside surface area is increased by about 140%. Except for cases where ¢ all a 
er vation is expensive, , it seems that the tank of depth 5.0 will be cheaper. pow =| : 
. - Using Fig. 1 of the analysis, it is possible to explain why the example turns 
out as it does. Between the depths of 0.8 it and 5 ft, the curve of removal for — 
a constant detention time of 2 hr is steep. The residual increases from 10% “ae 
_ about 17%, whereas theparticles have to settle almost more than five times a 
7 distance in the same length of time. As the depth is increased, the flocculation 
and the settling velocities are increased, allowing the particles to almost cover 
ca the increased distance. It is the last part of the settling that is the fastest and 
increases the distance covered. Consequently, reducing the tank depth 
- from 5 ft to 0.8 ft reduces the distance that the particles must settle, but it _ 
makes use of the earliest and, hence, lowest settling velocities. . Therefore, the — 
effectiveness of the settling is less in the shallow tank. 
iS _ The analysis by Gunnerson emphasizes that = suspension analyzed should — 4 : 
pee be the actual suspension that will do the settling. The settling analyses were ; 


performed on effluents from treatment plants, whereas the suspensions to ” « 
_ considered were the effluents mixed with sea- water. The addition of sea- -water 
produced a suspension of different characteristics , apparently with ell 
tendency toward flocculation. As mentioned by Gunnerson, the effect - sea- 
water should be studied as part < of the analysis. 
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UNSTEADY | FLOW OF GROUND WATER INTOA SURFACE 


WILLIAM HAUSHILD!® and GORDON KRUSE,29 4. M. ASCE.—The discus- 
sions indicate interest in the material presented in the paper and have con- aa 
“tributed toward expansion and understanding of the ne for th the nonlinear 
case of flow of ground water into a surface reservoir. 
The discrepancy shown in Fig. 2 between the experimental data and Solu- 
_ tions I and Ill, which were based on Eq. 7 is indicative of the error resulting . 
from neglecting the nonlinear terms in Eq. 6 ‘if dis small. Solutions Il and IV 
are based on the assumption H << d, inasmuch as they were developed from | a a 

Solutions I and III; however, they are not limited by the assumption because - | 


 - 


they consider the exact saturated thickness of the aquifer through which the 
flow occurs. The small amount of data available for comparison with Solu- 
tions II and IV indicates that their use might be preferable to the error- -function 
‘solutions, Solutions I and Il. . The experiment to which the solutions were ap- 
‘plied was a ‘complete drawdown case and dwas not large compared to H. Solu- an 
ry tions II and IV agreed well with the experimental data although the assumptions | § | 

were not completely satisfied. . When the drawdown is small, Solutions II and 

IV give results which are nearly. equal to those obtained from Solution i There- . 
a fore, the use | of Solution II or IV should g give good results over the full range - . _ 


‘Harr, in considering the steady s state form of Eq. 7, ‘slilaiiaalia the water a - 

a table condition before the reservoir surface is lowered. Therefore, the solu- — 

. tion would be expected to be linear, Eqs. 24 and 25. Eq. 28, the solution t 
Eq. 26, corresponds with Soltuion IV, presented by the writers , only for the 

case of complete drawdown. The consideration of Eq. 7 in the form of Eq. 26 _ 

4 seems questionable. If 7 is the difference in elevation between the water ta- 


ble and the underlying impermeable layer, the flow through a vertical section | if 


ag « 


: & July 1 1960, by William | Haushild and Gordon Kruse. (Proc. Paper 2551) Sa 

Hydr. Engr., Quality Water Branch, U.S.G.S., Fort Collins, Colo, 

20 Agric. Engr., Soil and Water Conservation Research Div. we Agrte. ne Serv- 
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Now considering 


Comparison of Eq. 48 with Eq. 26 shows E Eq. 4. 26 not represent the physical 
_ conditions of the problem. . Also the a a used in Eq. 30 is te nat the same as the a 
TABLE 3.— -BANK STORAGE DATA— ASE" 


K = 6. 034 


an _ Kazmann emphasized that the paper of Keller and Robinson8 contained only 
steady state data. The only non-steady state laboratory data collected by Kel- 7 ; 
ler and Robinson were obtained directly from Robinson’s files by the writers. x ie. 
- Unfortunately they did not note that the data were not published. However, the 
ner procedure and equipment is described by Keller and Robinson. - 
7 ‘The non-steady state data shown in Table 3, with permission of Keller, were > 
collected at a time of 6.5 min after the start of the experiment. The elevation 
of the water table in the aquifer was being lowered at the time of the -measure- 
. _ ments; however, the water surface in the headbox was still unaffected. ie 
_ Although the writers believed that the error-function solution of Eq. 7 was i - 
well known by the they are indebted to Kazmann for the reference 
g nll,12 of the error-function solution 
in connection with ground water studies. 
_ The writers have not assured the profession that the results may be po eal 
in the field. Rather, the application to each particular field location must be 


; 
. 4 at = q 
— ‘4 


considered with to the e relation field conditions and the experi- 


ait, field hardly be used to prove or. disprove the so- 
- lutions presented in the paper. However, field studies will indicate the useful- te 
ness of the solutions. Any applicable field data with which the formula — > 
-Moody’s solution is a valuable contribution and should be useful and easily 
in the tabulated form of Table His solution agrees well with Solu- 
tion II presented by the writers who agree that sufficient wenureny is obtain- 
by second degree interpolation in Fig.4 and Table2, 
- ‘The experimental data of Keller and Robinson are approximated moreclose- _ 
a by Solution IV of the writers than by Eq. 37. _ This may be due to error in ; 
the determination of a. Ina laboratory flume, with an appreciable 4 


te 
The writers a are grateful to the ‘method of 
tion approach to approximating the exact solutions of Eq. 6. It is interesting 
_ that this method and the approach used by the writers lead to identical results. a 


Cor rections .—The ‘second term on ‘the left side of Eq. 6 as presented 
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\ - Dis cussion by Pa LF. Keim and Vit and 


PAUL F. KEIM,” F. ASCE.—Einstein’s observations bring to mind the fol- 


7 In 1939 the writer, then with the Federal Power Commission had occasion © 


for navigation, “flood control, ‘and power. . No less than : a 1 dozen professional 
: a menand local farmers while standing with the writer on various bridges across 
tributaries in the headwaters of the river sadly related that they once fished 
= in n crystal clear and quite swift water. 
up to the lower level of the trusses with practically no channel space weeteo- : 4 
ever for the muddy flow. Certainly this is repeated in many other places as a 
‘result of the profligate use and unscientific development of lands and water, 
As an aside, it is not only silt which has clogged clear streams; witness ” 7 


disposal of waste in many areas, 


_ The author rightly states, “most of the artificial storage for water lel . 
is provided in the form of surface reservoirs. These reservoirs change the | 
regime of the rivers drastically particularly with respect | to their sediment 4 
characteristics: ” That being the case, it is a certainty that the relative + 
stability of channel shape and slope that has been reached under natural con- 
_ ditions will be disturbed. How much disturbed and what the resulting condition “ 
: -3 occasioned by the disturbance will be as of 1961 is largely guess work. The “A 
= of Einstein’s conclusion is s unquestioned if it isdesired to obtainthe _ 
“i In Egypt it is not one but a series of detention and diversion dams on the © 


Nile above and below each of which sediment has created aeae problems — S 


have been if Einstein’ s had been formulated and asa con- 
part of the punning wad river system development 
"equilibrium is results of which cannot be anticipated unless a 
are made to determine the range of possibilities. 
ee _ -The engineering approach philosophy as “developed by Einstein is a contri- 
- bution to both scientific and engineering activity. . if heeded, this philosophy ~ 
can be used to great advantage in the planning of an orderly and ae de- » 
; 2 Prof., Civ. Engrg. Univ. of California, Berkeley, Calif., on leave to Dept. of State % 
‘Tech. Advisor, USOM/UAR, Cairo, Egypt (UAR). 
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VITO. VANONI, 3 ASCE. — he author two important sets 
of sedimentation problems in connection with Conegnt of water resources | 
df in general, on reservoir development in particular. _The writer agrees with 
aa (1) these are important problems, (2) existing knowledge ‘is 
ae inadequate to on with them, and (3) research is needed in this area. These — 
ce, are problems of of sediment transportation under ‘unsteady and non-uniform con- A : 
_ Transportation under steady uniform | conditions, which has received 
a the major part of the total effort in sedimentation research, is farfrom solved, ne 
and the more d difficult non- - steady and non- uniform transport problem has hard- 
a been touched. "Major « dams are being built at an increasing rate, and new . 
onesare being authorized at costs running into many tens of millionsofdollars. 
As the dams are completed, the sedimentation problems grow more | numerous. 
° bake and the need for information to. cope with them more pressing. Activities in 
igen sedimentation research, although increasing in the past two decades, show only © oe 
modest increase in effort, compared to dam construction and the importance 
of the sedimentation problems accompanying their construction. It is urgent 
= § that adequate research programs be started at once to deal with these and oth- i ; 
author has indicated that in dealing with the sedimentation problems of 
4 reservoirs, he was assuming the position of the engineer. The writer agrees pe 
that if the answers needed to guide engineering work are to be obtained in time, ii, 
the engineering approach must be pursued. . However, the writer would like — 
argue that it would be unsound to followthe engineering approach at the exclu- 
- ston of basic research into the transport phenomena. Such basic research es : 
yields information slowly, but it forms the basis for improving methods for batt 
a te dealing with streams. It will provide a basis , not only for dealing with | prob- y! 
: lems that are known today, but also for the multitude of new ones that are cer-_ 
The ‘question arises, who is responsible for seeing to it that research is | 
_ done so information will be available to solve these sedimentation problems? - er 
ra Who shall take the initiative to get the research done? This wirter believes ~~ 
that to initiate research is the responsibility of the same people and organiza- a 
_ tions who build the projects, operate them, and cope with the oe ] 


problems. These organizations are the ones which know the magnitude and 
nature of the problems, and the success they are having in dealing with them 
_ Therefore, they are the best judges of the adequacy of their methods. — 
¥ ‘further within these organizations, it is clear that the professional civil engi- 
Bo! neer is intimately involved in this matter . He is the one the administrator de- 
pends on to estimate funds needed for construction and operation of projects, a 
and forresearch,aswell 
“ The funds will be provided by the organization, which in most cases is the _ 
public, but the program must be made out by the technical expert, in this case, | 
the professional civil engineer. Funds needed for research are nominal, com = 
pared to expenditures for project construction and operation. One item that 
_ is important is time. It takes time to train research engineers, and to do re- ‘'. 
search. Crash programs with large budgets are not the answer. For success, Raa, 
a support of ‘research is needed, which means that it is urgent that a 
Start on adequate programs of research on the outlined the i 


Assoc. Prof., Hydr. ., Calif. Inst. of Tech., Pasadena, ona, Calif. 


: 


ar 
&g 

| 

&§ 
| 
ay 
— 


P PARMAKIAN,”2 F, ASCE.—It is hoped that the following aditiona 
oe - examples where vibration occurred in certain hydraulic structures will am- 


‘merged for long periods. To ‘illustrate. extreme case, at on one installation 


Se _ there are three 12-ft 7-in. wide by 23-ft 5-in. high fixed wheel gates which are a 
ik used to shut off the flow in each penstock intake. ‘There are twenty-two 15-in. 
a or diameter forged s steel wheels on each gate. The wheel bushings are of the self- 7 
lubricating, bronze type, with graphited inserts and are not considered to re- es ; 
he _ quire any further maintenance ¢ after the initial installation. The 4-1/2-in. di- _ 
ameter wheel pins are chrome-plated SAE 1040 steel. These gates are stored ¥ 
under water above the penstock intake openings. After being submerged for © 
a4 pcs 10 yr it was necessary to make an emergency shut-down of a set of these 3 
intake gates | with nearly full flow passing through — the turbine. The first two 
intake gates closed completely without difficulty. However, there was little 
_ unbalanced pressure on these gates since the flow through the remaining i 
_ing was sufficiently large to supply water through the turbine. As the final in. 
take gate was lowered under unbalanced pressure conditions, the gate stem and 
_ crane started to vibrate when the gate was about one-third closed. The vibra 
soon reached violent proportions as the gate continued to | al 
_ The gate finally stopped when it was still approximately 8 ft open. = 
__ The gate was later hoisted to the surface for inspection. An ciaiiinaaiic . 
‘the twenty-two wheels indicated that only one of the wheels had rotated during 
_ the partial gate closure and that all of the other wheels had skidded on gal 
_ track. By using pinch bars and sledge hammer blows it was possible to free — 
- most of the wheels so that they could be made to turn. However, several of the 
, waaets could not be freed at all ~— had to be disassembled. A similar condi- 


os timating the coefficient under field conditions, especially when the gate is sub- a 


was due primarily to the failure of the wheels to rotate and roll on the track 4 
after having been submerged for a number of years. . The graphited inserts — 
- were subsequently scraped off and grease fittings installed on all wheels. Pe- — 


-riodical greasing 0 the pins is now required to that the 


Tech. ‘Analysis Branch, Div. of Design, , USBR, Denver, Colo. 
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Hence the failure of the intake gates to close wonder unbalanced flow conditions 
q * 
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In the section “Other Exciting Forces” the author indicates that the: fixed 
cavity which forms on the trailing side of an obstruction in high velocity flow 
may produce an intermittent pressure pulse and hence a source of “ate 
example of this” type of vibration occurred at a high-head outlet works 
;  —-where | severe vibrationand noise were present at only one of three outlet pipes. — 
- occurred when the pipe was discharging in excess of 0.6 of its rated capa- 2 
city. The vibration of the pipe and the noise were considered serious enough 2 
to restrict the operation of the outlet pipe to not more than half of its full . 
! city. | When field conditions permitted, the outlet pipe was unwatered and the 
upstream ring follower gate frame and leaf inspected. With the gate latched in 
; its open position the ring follower leaf was found to be 11/ 32 of an inch a 
™ This irregularity i in the alinement of the circular opening in the leaf with | a 
respect to the conduit caused sufficient turbulence in the high velocity flow 
_ downstream from the gate leaf to produce the objectionable v vibr 


steel 136 ft long. "Vibration of the gate during flow regulation oc- 
curred from its initial operation. Various remedies were attempted in an ef- 
*< _fortto prevent the vibration. Initially the small guide | shoe clearances provided — 

de some friction damping and the vibration was not too objectionable. _ However, be 

7 as the wear of the shoes increased these clearances, there was a correspond- — 

oe ing increase in the vibration amplitude. . Spring loaded shoes were tried a 
= were found to be ineffective. _ The vibration was reduced somewhat by modify- 


ing the lower edge of the gate so as to discharge as a sharp-edged orifice. ea 
ss: Severe vibration under conditions of small gate openings with unbalanced 
‘heads was also observed ata 32-ft cylinder gate. ‘This 330,000- lb 
gate is suspended by three steel stems, 7-3/4-in. O. D. by 4-5/8-in. I. D., om 
approximately 330 ft long. In order to reduce the hydraulic downpull on the 
gate > the cylinder gate was modified bad part of tapered gate 
 Atanother installation a cylinder gate 20 ft 4i in. in diameter approx-— 
_ imately 85,000 lb is suspended on three 4 1-1/2- in. . diameter steel stems, 330 ft Rey 


long. During the ‘initial operation severe vibration of the gate ‘developed at 


small gate openings. The observed frequency was about 6 cycles per sec and 


oo amplitude approximately 1 in. At this installation aviscous damping 4 


__ together with a serrated gate bottom will be tried to ‘reduce the vibration Pee 
Me 


van 

— 
— flinder gate in- 
a _ gtallations, especially when the gate openings were small. At one installation _ 7 
a 
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VECTOR ASPECTS OF DYNAMIC SIMILARITY® 
Steven, 
CHARLES” E. BEHLKE,!2 ASCE.- —The state four difficulties 
~ which often arise in distorted open channel models. It has been the writer’s 
pega that each of the difficulties are also encountered ed in n non - distorted — 
steady flow is considered, assuming the Chezy resistance equation 


in which s is a displacement along a streamline, V describes the velocity, 


refers to the Chezy coefficient, his the depth, and z denotes the elevation of 


water surface above arbitrarily chosen datum plane. The writer has 
chosen the Chezy resistance for study purposes. 


‘The equations of motion must apply to the well as to the prototype 
‘The equations define the scales between model and prototype: 


(displacement along a streamline), 4, | 


™ 77 


| 

— 
il 
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which r is the radius of curvature of the streamiile, wie Ulstanc 
c &§ BG measured normal to the streamline, and V refers to the velocity along the &§ 

oe 7 ‘The first term of Eq. 45 is an acceleration and the other two terms repre- - 

—— - sent forces. In Eq. 46,the term on the left side of the equation is an accelera- gi 
w=) Gg  « tion and the term on the right side of the equation represents a force. Except _ 
: a ; for freely fallin rticles, fluid particles in open channel models do not ex- i 

ia 

| 

1 

: March 1961, by R. Koll aha Ww. L. neitmeyer oc. Paper 2763). 
= Assoc. Prof. of Civ. Engrg., Oregon State Univ., Corvallis, Oreg. 
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in which y, s, z, and C are e still of Eq. 45. Normal 


similarity is satisfied if in Eq. 48 


the ale i , the flow is 
similar “model prototype if, in addition of Eq. 50 


‘Thus, from Eas. 50, 51, (52 et 


in which a, is the acceleration scale | between model _ prototype. This sca le 


ratio need not be unity as indicated by the authors. 


Anker 
_ Eq. 53, which is the result of Eqs. 50,51, and 52, indicates that the acceler- 


ation scale is the same for both tangential and normal accelerations. The equa - 
u tions of motion, then, clearly indicate that the acceleration r ratio is not ae 


| 
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(57) 


_ There is nothing which limits any of the foregoing equations to equal hori- 

; 24 ‘zontal and vertical scales. — Irregardless of the distortion of the model, if the — | 
preceding scale factor equations are in the acts” 

w@ the model is too ‘rough or too © smooth, ho and the streamlines in 

_ the model are not similar to those in the prototype. If the model is too rough, — 
the model streamlines curves too much; and if it is too smooth, ‘streamlines | 

_ do not curve enough. This would be true for both distorted and non-distorted a 
models and can also help explain many of the which may 

A plot of transport versus velocity for : a specific depth can be made for the a 


¥, = 


a prototype and for the model. This usually indicates that the transport scale i 
changes with velocity. These ‘changes in transport scale could be nce an 
if a suitable bed material could be found for the model. _However, since an in- | 


— 


ot must usually be used. . The results obtained from this bed material must be - 
= carefully evaluated because it will react differently from the ideal material. | 
: _ The various things presented herein contribute more to the difficulties en- 
countered in distorted models than the distortion itself, and the equations given 


indicate that velocities must be scaled to the square root of the depth scale and — 


‘ 
x 


these equations, other useful equations cap be 
— 
q 


R. STEVEN,’ 13 M. ASCE. — -The authors present concepts 
which should be examined ee before they are accepted as — 
7 In Eqs. 6 , 7, and ¢ 8, the itd forces have be been taken as ene caused : 
by viscosity, pressure "and gravity. The authors have chosen, , for this analysis, 
to neglect the forces due to elasticity and surface tension. The inertial force 

hae been takenas the resultant of all forces and the usual dimensionless ratios 

F, R, and E have been obtained. _ The fact has not been recognized, however, 

a that the Euler number, which is ‘more accurately described physically as ~<a 

_ pressure coefficient, must be the same at all corresponding points in two _* 

-namically ‘similar systems regardless. of the e particular basic forces which ie 

_ tually shape the motion. Therefore, pressure or more significantly pressure 

a difference is not a basic or motivating factor that creates the fluid motion, but — 
rather, ‘it is a result of the motion itself. In general, in any solution in which © ee 
3 viscosity and gravity are the only basic factors of practical importance gid in 


than the boundary a generalized solution can be 


The classical this is is flow past a cylinder. When the fluid is ideal 
and the cylinder» is entirely immersed, | the _pressure coefficient is a function — 


identical at all points, This ‘clearly sho shows the dependent rather 

than primary role played by pressure in fluid motion. 
_ The authors have clearly emphasized the vector aspects of the equations. z 

Ba 7 conveys the sense e of a cs gps of vectors, in n three dimensions, the scalar _ 


‘remain The scalar of the unit vector diagram rectors 
are the dimensionless ratios R, F, and so forth. Obviously, for dynamic sim- r 
am between two systems, the scalar magnitudes of the corresponding parts | 
of the unit vector diagrams must be equal. Derived in this sense, Eqs. 13, 14 
a ont 1 have far greater physical mguitiennces | than when viewed merely as s the 
oe oe ‘The development of the Kolf number implies that this is a dimensionless — 
oe ratio, similar in nature to the Reynold’s, Froude, Weber, and Mach numbers in 
that it, like the others represents a ratio. between “basic forces at work ina 
system. This is not correct since the Reynold’s number and the others aie 


primary agencies or motivating forces to the resultant or initial for rce. The 


| 
ig 
coefficient 
factors. Consequently, division of each scalar magnitude by the scalar mag 
4 nitude Fj represents, physically, reduction of the actual spacial vectordiagram 
j 
4 
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Kolf number does not do this because centrifugal foree is not the: result ofa 
primary agency. The magnitude of the centrifugal force is determined entire- 
ly bythe streamline geometry and the velocity. It is thereforea resultant ai 
er than a motivating or basic element. Consequently, it is improper to suggest _ 
- Bh that the ratio V/u, which is simply an expression of kinematic similarity, has a ; 
_ important ratio since the direction cosines of the resultant velocity vector : 
~ neglecting radial components, are proportional to it. In fact this ratio in slight-— 
- different form has long been used in propeller analysis and hydraulic ma-— 


_ chinery under the name “slip function.” Take the forward velocity of the ship 
V and the tangential speed ¢ of the propeller tip as 7 D DN. ' The ratio 


“more usually expressed as: —~ is the slip function, This is simply the 
_ cal of the Kolf number. 
os 


vend 


= a common base with the Reynold’s and Froude numbers. Actually V/uis an a 


>, 
= 


— 


: 


— Th. 
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- the past. About all that can be hoped foris a greater understanding of the pre- 
& sent or initial state of things to serve as a guide in the predicting of the future — q - 
‘state. . To this end, all analyses are aimed. The forecasting technique pre- 
ee sented by “Rockwood and Jencks uses some unique forms of analytical logic | a 
= _ which seem worthy of additional consideration and study. ‘They are not entire- 


ly separable but fall into the following general categories: 


of abstract indexes as s forecast parameters; 

2. 2. Qualified acceptance of “Goodness of Fit” tests; 


4. the relative significance of 


rea isa tendency by most forecasters to select pastaien’ for analysis which | 
have a commonly- accepted high degree of causal relationship to the physical oat 
vr id factor to be forecasted. The less obvious and direct the causal relationship is, 
a a the more abstract it becomes. 5 But there are various degrees of abstraction. be 
_ The use of winter runoff of coastal streams to forecast the spring runoff of — 
. Columbia River is probably a far more abstract relationship, on the basis of 
common acceptance, than the use of precipitation and snow course measure- 
"ments from the area where the runoff of Columbia River originates. 
authors have rationalized the advantages of their uniquely defined Coast- 
al Flow Index (CFI) over the use of er eee by noting the deficiencies of 
Se and snow course measurements. They seem to be pointing toward | 
a conclusion that the runoff basins of the coastal streams constitute more ad- 
a "equate precipitation gages for an area several hundred miles distant than i 
perie _U. S. Weather Bureau gages located in the area. The big argument a 
: in favor of this conclusion is that the» coastal runoff basins provide a much — 
per The deficiencies of. measuring ‘streamflow from the coastal streams is ra- 
tionalized as acceptable because deficiencies also in the measurement 
a other indexes such as precipitation and snow courses. _ Whenever rationaliza- . 
used to | justify abstract considerations, the question | is immediately 
raised regarding how much abstraction should be justified in this way. It has © 


‘been emphasized that Columbia | River seasonal runoff is partly a function of 
westerly air flow and associated weather phenomena ‘during the 


March 1961, David M. Rockwood and Carlton E. Jencks Paper 2780). 


a 
1; 
— 
| 
7 | 


1961 

"previous winter season. 1. This observed climatic peraiotency is the medium 

which the causal relationship is rationalized. Some of the other 

nomena are partial functions of the same climatic conditions. — _ However , they 

- would not necessarily be considered as causing the variations in Columbia 

the number of people using golf courses during 

_ winter months at Portland, Ore., would not be considered as causing the Colum-— 7 

bia River runoff to be greater or ‘smaller, although the same climatic = 

tency medium could be used to rationalize their relationship, = | +r 

_-‘ The influences affecting precipitation from the time it falls to the ground ~ 

until it reaches the gaging station are many and difficult to evaluate individual-_ 

ly. The use of CFI may assume a direct proportion between these many in- 

fluences on the coastal streams and similar influences farther inland. If this ¥ 

assumption is valid, it would seem that use of the ratio of coastal precipita- 
: tion to CFI as a correction index to Columbia basin precipitation might be a a 

less abstract causal relationship. _ This combination may conform more di- 

rectly to the water balance of the hydrologic cycle and therefore more aco ; 


~ Frequent references in the paper to values obtained from applying Eq. 1 a 
the data from which it was derived as “Forecasts” may lead to the conclusion — 
_ that the authors consider the “Goodness of Fit” of 2 an equation as a measure of _ 
its predictability, or how reliably an equation will forecast in the realm of fu-— 
_ture events. It is stated, however, that the parameters selected for use in the 
equation do not necessarily represent the best fit of historic data of all para- . 
meters investigated. The purpose of regression analysis is to to iis oul eal , 
of the data in accordance with the “least squares” theory. 
___ If the goodness of fit tests leading to the least squares solution is not entire- al 
; ly acceptable, there would appear to be little justification for use of regression i ; 
a analysis at all, except if the limitations inherent in the least squares theory ec 
were recognized as resulting in an equation which is incompatible with pre- — <= 
conceived notions of the physical interacting forces. 
s ‘It appears that such recognition may be the basis for selection of an equa-_ i 
tion which does not represent the best fit of historic data and also explain the » © 
use of product functions as an adjustment of data prior to introduction into the ‘ 
regression analysis. The authors have attempted to justify their use of — af 


relationship, it in n accounting for runoff variance. “ It is assumed that the term, 


elation relationship, is intended to mean the oceaaead of fit of the equation to ay 
the data from which it was derived. 


ae _ The physical p process relating the CFI to the . April runoff of Chehalis River 
near Grand Mound, Washington, is not apparent to the writer and has not been 
_ fully explained. If the use of produce functions, as shown, is valid for inde - a = 
-_- pendent variables, it would seem that similar product functions involving» the 
va dependent variable would also be valid because all the independent ‘oi ri 
_ were selected because of their known physical relationship to the dependent on qi 
variable. When the CFI was multiplied by the other measured physical factors, 
the resulting variables became unidentifiable in relation to measurable physi- is 
ae cal forces. The abstractness previously ascribed to CFI made an — 
ae at this point leaving a v a vast void d of causal explanation. rain 


&§ 
a 
~ 
a 
a) 
alee eal process relating two independent variables as a product. The use of prod- im 
ae 
i 
— 
— 
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"DISCUSSION DISCUSSION ON 
an mene for the and terms in Eq. 1 and 
s follows, results in Eq. 2 for use with data available on April 1: 
0.7428 - 0. 2836 X 0 


0, 2836 X> - 0.678 X 

Eq. 2 wi will provide identical (except for decimal rounding) forecasts in the 

realm of future events as Eq. 1. qt 7 also show equal: goodness of fit to oll 


3 cients ‘would. have been significantly different from those in n Eq. 2. "Such a an 
a, analysis would also have provided an improvement in fit of historic data. — 
__ Rockwood and Jencks have considered such improvements in fi fit as falla-— 
cious. . Evidently, this conclusion is based on the questionable assumption that 
& 1 contains allof the known significant factors affecting runoff and is there- *. 


fore considered to be hydrologically sound. When the results of one regression | 
7 ia can be considered fallacious, the question arises as to what results 
can be considered as valid. Certainly, the amount of serious effort contribut- 
_ to the selection of parameters and analytical technique should qualify the — 
authors to judge the adequacy of their particular analysis. - But there still re 
mains the need for a ‘commonly- acceptable and universally-applicable 


of judging the validity of regression analysis results. 


Obviously, all of the physical ; factors affecting a a dependent variable can nev-— 7 


er be directly included in a. regression analysis except in the realm of ideal- ih 
ized theory. Few, if any of the factors are measurable in an entirely 


ent state. With limited measuring techniques as presently exist, it is often de- 
he p ly 

sirable to index many of the more abstract factors such as "Rockwood and 
Jencks may have done by use of their product functions. 


_ Further algebraic of Eq. 2 is possible for evaluating the relative 


values for the partial indexes as 


Y =0. 6240 x1 +20, 
: 
39 (107. 6) - 0 0. 2836 


— 
nd factoring 7 

— 

that if the four variables in Eq. 2 had been introduced into the regres- 

4 
— 

% 
— i 

| 
g 
jae 


107.24 - 0.3052 


9839. 


_ in which X9, denotes the winter runoff of Columbia River near The Dalles, 

4 _ Oregon(mean v value of 34 -07 from Table 3) ) and, 2 is the four-month average 
winter temperature (mean value of 38.7 from Table 4). 


From Eqs. 3, 4, and 5 the po sie variances have been assigned the ~ 


Inspection on of the mean values of these four that come 
: from a different field of measurement and numerical classification. Presum- 
= “ably, the variance of any of the variables about its mean value would be pro- 
_ portionate to the recorded variance even if that variable had been measured by 
4 a standard resulting in a meanvalue equivalent to the mean value observed for fp 
Y. Some part of the apparent relative variances previously indicated is — 
_ ably due to the relative magnitudes of the differing numerical classifications 
of the observed data. By standardizing | the numerical classifications to the 


4 _ value of Y (My), 4 a further definition of the relative weights assigned by 
regression analysis to the ee measured factors may be o obtained an 


‘in which 107.6 and iM, = = 96.8 


107.24 = - 0.1074 


125.07 = - 0.7295 


> 


- 
— 
— 
x — 
| 
in which Mgg 34.07 and M | 


HY5 DISCUSSION 


portionate varia variance significances a as 


ine algebraic factoring by which Eqs. 2 through 8 were derived is incom- 
_ patible w: with regression analysis because of differing philosophical bases. They 


to various conclusions about the relative significance of a contributing vari- i 
* able. Actually, equations derived from regression analyses are never com- — 
e he plete in themselves when viewed in the algebraic sense. The complete state - 


- ment contains an evaluation of how well the equation fits the universes of past 


_ and future. Such an evaluation in the complete sense can only be qualitative. 
’ Otherwise, algebraic equations are never complete in the probability — 
a the philosophy in which regression 1 analysis has its being. . Algebraic ¢ equation 
an 5 imply absolute equality, or perfection ‘The philosophy of probability denies 
absolute perfection. 


7 Victor A. Koelzer, M. ASCE, and Perry M. Ford havel5 somewhat arbitrar- 


rs ily ily considered a regression coefficient, b, as significant : if its standard error, 
_ Sp, was equal to or less than b/2. On this basis, the ratio 5,/b yoctetinn nu- =, 


eri ‘merical values inverse to the significance of the b coefficients and 1 - 
provides numerical values proportional to the significance. The latter form 
used | here for comparison t to the previous lists of proportionate variance sig- 
nificance. 4 Applying this concept to Eq. 1, the relative of the re- 


r+ the bg and bg coefficients becomes doubtful. Unqualified acceptance ir re- a 
SS significances of the independent variables as shown in Figs. 2 (a) and 2 s 

‘a may result in conclusions as misleading as would result from acceptance of 

- etther of the three tabular listings presented herein. The true relative causa- 

_ tion significance of the variables is not defined by regression analysis. a cao 

; ‘The author’s reluctance to completely rely on best-fit tests, and their use 

q  ¢ of product functions have strong intuitive appeal. Their bases can be simply 

demonstrated as follows: Given the measured values of Xj in Table 10 to fore- 

. = the measured values of Y, the first step may be to construct a table of Xy 


15 “Effect of Various ‘Factors on Snowmelt Runoff,” by Victor A. 
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— 125.07 = - 0. 

— — 
4 

a j 
| 

\ 
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times Y. In the future, whenever a value of is given, a forecast of Y could 
be obtained by finding the corresponding value of YX, and divide by the x ~" - 


ae In effect, this technique providesas many coefficients as there are samples, 4 = 2 


va) 


each of which is a product-function of the variables. Observed X; values in 


the future which are f fractional parts of the X, values in the sample, canbe 
- used by interpolation. Evaluation of “Goodness of Fit” for the sample in Table © 

io D provides perfection. If, however, the values of Y are always a function ot 


the =X1+X , future values of ‘may be more affected the 


ine 


value which was not considered in 10. A future value = 3would pro- 
a vide a Y value of 103 if Xp = = 10, whereas the procedure in Table 10 would pro- a 
vide a forecasted Y value of 147, _ A procedure which does not evaluate all the - 

ie causal effects cannot be ) completely eval evaluated ted for forecast reliability by good- 

The techniques. presented by Rockwood and Jencks constitute a bold step 

il forward in the use of abstract parameters for runoff forecasting. _ The ade- — 

a quacy of Eq. 1 will become more conclusive, however, as the application vand 


verification of samples are 


vor? 
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ROUGHNESS | IN RIGID sD OPEN CHANNELS 
Discussion by P. Biery and w. Delleur 


F. BIERY 28 a, M. ASCE and J. W. DELLEUR,29 M. ASCE.—The writers 
: _ Wish to extend the paper by showing the result of applying Sayre andAlbertson’s 
, analysis to a different type of roughness element consisting of round bars,and ~ 
aa tests were performed in a steel tilting flume 5 ft wide, 2 ft deep and 64 — 


ft long » rm me tests were run with two different boundary roughness 


longitudinal have Pesan 12 in. on ‘center and a 
in. on coutne, (b) along the side walls one layer of vertical bars 6 in, on center 
placed + in. from the wall. The bottom layer of bars were tied together with - 
wire. "The vertical bars were tied at the bottom to the transverse bars and 
clamped to the walls above the free surface. _ Fig. 13 shows the artificial — 
gt! ‘Uniform flow tests were run for smooth and rough | boundaries. The Darcy 


Weisbach friction factor, f, was computed from the oqention : 


to the "slope. In Fig. 14, the factor, ‘f, is plotted ver versus the 


v is the kinematic viscosity of t the flui 
‘The roughness elements used here are different from those used by eons: 
and Albertson. In particular, there is a definite amount of flow under the a= 
_ roughness elements, _ Fig. 15 shows a qualitative sketch of the flow around | E 


_ transverse bars. Centerline velocity profiles measured close to a transverse pe 
bar and ata point midway bars are shown | in . 16, for 
which h Q= 3.714 cfs and § 8= = 0.( 0125 _ * 
oe May 1961, by William W. Sayre and Maurice L. Albertson (Proc. Paper 2823). 
28 Research Engr. , Johns Manville Corp., Manville, N. J.; formerly Research Asst., 
School of Civ. Engrg., Purdue Univ. Lafayette, Ind. 
Assoc. ‘Prof. of Hydr. -Engrg., School of | Civ. v. Engrg., F Purdue ‘Univ. Lafayette, Ind 


— 
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| 
ae boundary, consisted of the steel flume walls finished with an epoxy resin paint. _ 
i| The second roughness pattern, which will be referred to as rough 
— 
— 
— 
— 
CO 
— 
= 
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September, 1961 
ee as) Six tests we were run to > determine the roughness parameter, ; In order to 
th have fully rough turbulent flow, the flume was set to its maximum slope i 7 


ra 


0.0125. . The test data are given in Table 2. 


plot of log n/a similar to Fig. 5 was prepared. Taking the 


with a slope of 6.06 confirming constant in Eq. 1%. The extrap-— 
_ olated value of C2 was 3.15. With these values of thea and a, x was determined A 


Centerline profiles were taken at a slope of 0. 0125 and a 


— 
| 

— 
— 
— 
is: compared to the veloc i 1 Eq. 20. 
) the velocity profile presented in Fig. 9 and Eq. 20. 
0.0126 Q= 3.714 cls, 8 = 0.0125 ant 
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_ ‘FIG. 14.-f-RRELATION FOR NORMALDEPTH TESTS 


September, 


was 


log * + 4.6 


With the change pattern the first of the empirical 6. 06, 
_ checked; but the second constant changed from 2.6 to ‘difference is 


0. 03 ft downstream of a cement 


~ 
= 
(=) 


Velocity, in feet per second 


FIG. 16. —EFFECT OF BARS ON VELOCITY 
that the roughness baffles used by Sayre and Albertson were ‘placed in such a 


g that there was no flowbeneath the roughness elements, whereas there was a - 


, certain amount of flow underneath the transverse bars” used in experiments 
. _ If Eq. 20 is accepted for the bar roughness, it would be possible to find the 


“value an equivalent roughness height, a, for the round bar roughness. Equat- 


— 
aig — 


P, Eqs. 20 and 34, the equ equivalent roughness p parameter, xX for the round — 


= for a, one obtains a =0. 0195 ft =0. 234 in. |, which is close to the a Ls 
ter of the bar of 0.25 in. Eq. 20 for the velocity distribution may also be used 
2 round bar roughness with a reasonable degree of accuracy by consifering 


equal to the diame her of of transverse bars. 


“TABLE 2, _TESTS FOR. THE ‘PARAMETER _ 


| 666 829 22.548 
44 650 | 21, 
960 


574 
7,283) (4,779 


(b) Velocity Profile Data (y measured from the bottom)@ 


| 6.556 
490 
3. 
8,48 
3,65 
9 \ 


| 


Fig. 18 shows a portion of the general diagram of Fig. 10, with» 
4 test data for the bar roughness added, in which the values of y/x indicated Je 


: + correspond to a value of x of 0.0126 ft. ” There isa ‘aad good a agreement 
substituted for y in the smooth boundary data). 
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@ Tests at maximum flows and slopes 


18.—GENERAL RESISTANCE DIAGRAM | FC 


y 


34, and Ve is the friction velocity ve Vey, 5 
4° at 0.2 and 0.8 time the depth 
> 


It is probable that the roughness parameter, x, may also be use 
and Sy, which could be determined from velocity measurements at 0.2y,, 
measurements, Eqs. 20 and 34 can 
— 
in which 6.06 log 1/¢ is equal 


06 Vy log — . 


Taking the ratio of Ses. 36 and 31, ating “Vo 2/%o. 8 and solving for x, 


The determination of x 
depend on the of € which is used. For the 
of. this analysis, using Eq. 34, 4.6 = 6.06 log 1/e and | € =0. 174, _ With this val- — 

ue of €, , Eq. 38 gives (with w = 0.7) for the velocity profile of Fig. 16 a value of ia 

x of 0. 0123 ft which compares favorably with the value of 0.0126 previously — 

obtained. instead Eq. is used, =0. and x is found to be 0. 0058 ft 


coefficient x, parameter €, 


G (U) = (log us 3 (log U)? - 6 
wry on the velocity profile of Fig. 13, the value of a computed by Eq. 4 
The authors have shown that Eq. 17 is more accurate than Manning? s for- 
a - mula over the range of conditions tested. _ The writers have shown that Eq. 17 
is also applicable to a different type of. roughness, and that the x 
_ may be used for field conditions where it can be obtained from velocity mea- 


‘ surements at two and eight tenths of the e depth. . The writers hope that suffi- 7 


cient information on the roughness parameter, x, may be collected in the near = 


Be i future so that designing engineers can use it reliably for field channels and a 
: natural streams, perhaps even including channels in alluvial terrains. Ree 
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38 gives the roughness parameter, x. in terms of the ratio of the velocity 
— 
, @, may also be given in terms of 
4 and the ‘ratio of the velocities at 0.2 Yp and 
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